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PROPAGATION OF SINGULARITIES IN MANY-BODY 

SCATTERING 



0\ ! ANDRAS VASY 



Abstract. In this paper we describe the propagation of singularities of tem- 
pered distributional solutions u G S' oi {H — X)u = 0, A > 0, where if is a 
many-body Hamiltonian _ff = A-fV, A>0, V = X]^ ^a, under the assump- 
^^ , tion that no subsystem has a bound state and that the two-body interactions 

Cn ' Va are real- valued polyhomogeneous symbols of order —1 (e.g. Coulomb-type 

with the singularity at the origin removed). Here the term 'singularity' pro- 
vides a microlocal description of the lack of decay at infinity. We use this 
. result to prove that the wave front relation of the free-to-free S-matrix (which, 

■^r ' under our assumptions, is all of the S-matrix) is given by the broken geodesic 

flow, broken at the 'singular directions', on §"~^ at time tt. We also present 
a natural geometric generalization to asymptotically Euclidean spaces. 



Propagation des singularite:s dans la probleme de diffusion a A'" 

CORPS 



Resume. Dans cet article on decrit la propagation des singularites des solu- 
tions temperees n S <S' de {H — X)u = 0, A > 0, oil ii est un Hamiltonien a N 
corps H = A + y, A > 0, V = ^^Va, ea supposant que les Hamiltoniens des 
sous-systemes n'ont pas de vecteurs propres (dans L'^), et que les potentiels a 
deux corps Va sont des symboles polyhomogenes reels d'ordre —1 (par exemple, 
de type Coulomb, mais sans la singularite a I'origine). Ici le terme "singularite" 
fournit une description microlocale de la croissance des fonctions a I'infini. On 
emploie ce resultat pour montrer que la relation de front d'onde de la matrice 
de diffusion, A^-amas A^-amas (qui est la seule partie de la matrice de diffusion 
r^ ' sous nos hypotheses), est donnee par le flot geodesique brise dans les "direc- 



geometrique naturelle au cas des varietes asymptotiquement euclidiennes. 



X 

H 1- Introduction 

In this paper we describe the propagation of singularities of generalized eigen- 
functions of a many-body Hamiltonian H = l^ + V under the assumption that 
no subsystem has a bound state. We use this result to prove that the wave front 
relation of the free-to-free S-matrix (which is the only part of the S-matrix under 
our assumptions) is given by the broken geodesic flow, broken at the 'singular di- 
rections', on §"^^ at distance tt. We remark that these results have been proved 
in three-body scattering, without the assumption on the absence of bound states, 
in [Q. Also, Bommier |^ and Skibsted |^ have shown that the kernels of the 
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2 ANDRAS VASY 

2-cluster to free cluster and 2-cluster to 2-cluster S-matrices are smooth, and pre- 
viously Isozaki had showed this in the three-body setting Q. However, as is clear 
from the smoothness statement, the microlocal propagation picture that is crucial, 
for instance, in the discussion of free-to-free scattering, does not emerge in the 
previous examples when the initial state is a 2-cluster. 

In this section we discuss the setup in the Euclidean setting, but in the following 
ones we move to a natural geometric generalization introduced by Melrose in pl[ . 
Namely, suppose that X is a manifold with boundary equipped with a scattering 
metric g and a cleanly intersecting family C of closed embedded submanifolds of dX 
with Co = dX G C. Thus, g is a, Riemannian metric in int(X) which is of the form 
g = x~'^ dx^ + x~^h near dX; here x G C°°(X) is a boundary defining function. 
We also assume that near every p e dX, C is locally linearizable (i.e. in suitable 
coordinates near p, every element of C is linear) ; this holds if every element of C is 
totally geodesic with respect to some metric (not necessarily h) on dX. Let A be 
the Laplacian of g and suppose that V € C°°([X;C];IR) vanishes at dX \ [J{C € 
C : C y^ dX}, and H — A + V - we refer to Sections |] and |6| for a more detailed 
discussion of the geometric and analytic aspects of the setup. We prove under the 
assumption that there are no bound states for each of the subsystems (we describe 
the assumption more precisely in Section ^, but it holds for example ii V > 0) 
that singularities of solutions u e C~°°{X) of {H — X)u e C°°{X) propagate along 
generalized broken bicharacteristics of A which are broken at C . We also show that 
this implies a bound on the singularities of the kernel of the free-to- free S- matrix. In 
effect, we show that many-body scattering is in many respects a hyperbolic problem, 
much like the wave equation in domains with corners, for which the propagation of 
analytic singularities was proved by Lebeau |0| . The geometrically simpler setting, 
where the elements of C (except Cq = dX) are disjoint, corresponds to three-body 
scattering in the Euclidean setting, and then the analogy is with the wave equation 
in smoothly bounded domains, where the results for C°° singularities were proved 
by Melrose and Sjostrand [p2|, p3J and Taylor [Esl, and for analytic singularities by 
Sjostrand ||ll. 

Here however we caution that another important aspect of typical many-body 
systems is the presence of bound states of subsystems. While propagation theorems 
indicate that geometry plays a central role in scattering, bound states afford a 
similar role to spectral theory. Thus, in general, the two interact, even changing 
the characteristic set of the Hamiltonian. The generalized broken bicharacteristics 
are also more complicated in this setting, and, as a quick argument shows, the 
'time tt' part of our result will not hold if bound states are present. In addition, the 
Hamiltonian must possess additional structure (as the Euclidean ones do) so that 
propagation in bound states can be analyzed. Hence, in this paper, it is natural to 
impose our assumption that there are no bound states in the subsystems. 

We now return to the Euclidean setting. Before we can state the precise defini- 
tions, we need to introduce some basic (and mostly standard) notation. We consider 
the Euclidean space M", and we assume that we are given a (finite) family X of 
linear subspaces Xa, a E I, of R" which is closed under intersections and includes 
the subspace Xi — {0} consisting of the origin, and the whole space Xq = M". Let 
X"" be the orthocomplement of Xa, and let vr" be the orthogonal projection to X"', 
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■Ka to Xa- A many-body Hamiltonian is an operator of the form 
(1.1) iJ = A + ^(7r'^)*K; 

here A is the positive Laplacian, Vo — 0, and the Va are real-valued functions in an 
appropriate class which we take here to be polyhomogeneous symbols of order —1 
on the vector space Xa to simplify the problem: 

(1-2) Va e ^phglX'^). 

In particular, smooth potentials Va which behave at infinity like the Coulomb po- 
tential are allowed. Since {■K°')*Va is bounded and self-adjoint and A is self-adjoint 
with domain _ff^(R") on L^ = L^(R"), H is also a self-adjoint operator on L^ with 
domain ij2(R"). We let R{\) = {H - X)'^ for A € C \ R be the resolvent of H. 

There is a natural partial ordering on / induced by the ordering of X'^ by inclu- 
sion. (Though the ordering based on inclusion of the Xa would be sometimes more 
natural, and we use that for the geometric generalization of many-body scattering 
starting from the next section, here we use the conventional ordering.) Let /i = {1} 
(recall that Xi = {0}); 1 is the maximal element of /. A maximal element of /\ /i 
is called a 2-cluster; I2 denotes the set of 2-clusters. In general, once I^ has been 
defined for fc = 1, . . . , ttt, — 1 , we let /,„ (the set of m-clusters) be the set of maximal 
elements of /^ = ^\Uj!!7/-^fc, if -^m is not empty. If /^ = {0} (so I'^^+i i^ empty), we 
call H an m-body Hamiltonian. For example, if / 7^ {0, 1}, and for all a, 6 ^ {0, 1} 
with a y^ bwe have XaOXi, = {0}, then _ff is a 3-body Hamiltonian. The A^-cluster 
of an A^-body Hamiltonian is also called the free cluster, since it corresponds to the 
particles which are asymptotically free. 

It is convenient to compactify these spaces as in [Q. Thus, we let §" to be the 
radial compactification of R " to a closed hemisphere, i.e. a ball, (using the standard 
map RC given here in ( ^.3[ )), and §"^^ = 9S". We write w = ruj, uj £ S"~^, for 
polar coordinates on R", and we let x e C°°(S") be such that x = (RC^^)*(r^^) 
for r > 1. Hence, a; is a smoothed version of r^^ (smoothed at the origin of R"), 
and it is a boundary defining function of S!f:. We usually identify (the interior of) 
§^ with R". Thus, we write S;^g(§+) and S'^g(R") interchangeably and we drop 
the explicit pull-back notation in the future and simply write x = r'^^ (for r > 1). 
We also remark that we have 

(1.3) 5p™hg(§'^) = x"C°°(§'|). 

We recall that under RC, ^^(S"), the space of smooth functions on §" vanishing 
to infinite order at the boundary corresponds to the space of Schwartz functions 
5(R"), and its dual, C-°°{E:l), to tempered distributions 5'(R"). We also have the 
following correspondence of weighted Sobolev spaces 

(1.4) H^c^iK) = H^'^ = H^'\^'') = {wy^H^{W) 

where {w) = (1 + liyp)-*^/^. Thus, for A G C \ R the resolvent extends to a map 

(1.5) R{\) : H^;i\^l) ^ Hi+^\^-+). 
Similarly, we let 

(1.6) Xa = C1(RC(A„)), Ca=Xan d^^. 



4 ANDRAS VASY 

Hence, Ca is a sphere of dimension ria — 1 where n^ — dvcaXa. We also let 

(1.7) C^{Ca: ae/}. 

Again, we write the polar coordinates on Xa (with respect to the induced metric) 
as Wa = TaUJa, ^a S Ca, and let Xa = r~^ (for Tq > 1). We note that if a is a 
2-cluster then Cq n Cfc = unless b < a. We also define the 'singular part' of Ca as 
the set 

(1-8) Ca, sing = Ub^a(Ch n Ca), 

and its 'regular part' as the set 

(1-9) Ca = Ca\ Ub^aCb = Ca\ Ca, sing- 

For example, if a is a 2-cluster then Causing = and C'^ — Ca- We sometimes write 
the coordinates on Xa ® X'^ as {wa,w°-). 

Corresponding to each cluster a we introduce the cluster Hamiltonian ha as an 
operator on L'^{X°-) given by 

(1.10) /ia-A + ^14, 

b<a 

A being the Laplacian of the induced metric on X". Thus, if i/ is a iV-body 
Hamiltonian and a is a /c-cluster, then ha is a {N + 1 — fc)-body Hamiltonian. The 
L^ eigenfunctions of ha play an important role in many-body scattering; we remark 
that by Froese's and Herbst's result, Q], speCpp(/ia) C (— oo, 0] (there are no positive 
eigenvalues). Moreover, speCpp(/ia) is bounded below since ha differs from A by a 
bounded operator. Note that X^ = {0}, h^ = 0, so the unique eigenvalue of /iq is 
0. 

The eigenvalues of ha can be used to define the set of thresholds of hb- Namely, 
we let 

(1.11) K = Ufc<a SpeCpp(/l6) 

be the set of thresholds of ha , and we also let 

(1.12) A^ = Aa U SpeCpp(/la) = U6<a SpeCpp(/lfc). 

Thus, G Aa for a ^ and Aa C (— oo, 0]. It follows from the Mourre theory (see 
e.g. ||, |2^) that Aa is closed, countable, and speCpp(ft.a) can only accumulate at 
Aq. Moreover, R{X), considered as an operator on weighted Sobolev spaces, has a 
limit 

(1.13) i?(A ± zO) : Ht'i^l) - H^+'''\Sl) 

for I > 1/2, I' < —1/2, from either half of the complex plane away from 

(1.14) A = AiUspeCpp(i/). 

In addition, L^ eigenfunctions of ha with eigenvalues which are not thresholds are 
necessarily Schwartz functions on X"^ (see [Q). We also label the eigenvalues of ha, 
counted with multiplicities, by integers m, and we call the pairs a — (a, m) channels. 
We denote the eigenvalue of the channel a by e^, write ipa for a corresponding 
normalized eigenfunction, and let Ca be the orthogonal projection to ipa in L'^{X"'). 
The definition of the free-to-free S-matrix we consider comes from the stationary 
theory, more precisely from the asymptotic behavior of generalized eigenfunctions, 
see 1 34 , and cf. [£l| |9|. Apart from the difference in normalization, it is the same 



as the S-matrix given by the wave operators, see |B8|. For simplicity, we state the 
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asymptotic expansion under the assumption that Va is polyhomogeneous of order 
—2 (so it decays as jw^j^^). Namely, for A G (0,oo) and g G C^{Cq), there is a 
unique u £ C~°°(§") (i.e. u G 5'(R")) such that {H — X)u — 0, and u has the form 

(1.15) u = e-'^'-r-("-i)/2^;_ + i?(A + iO)f, 

where v^ e C°°(§'f:), u_|s„-i = 5, and / e C°°(§'|). In addition, this u is of the 
form 

(1.16) u = e-*^''r-("-i)/2t;_ + e''^W-'^"-^^/^v+, v+ e C°°{§1 \ Co,sing)- 
The Poisson operator with free initial data is the operator 

(1.17) Po,+ (A) : C^iC',) ^ C-°^{Sl), Po,+ (A),g = u. 
Following pJ, we define the free-to- free scattering matrix, S'oo(A) as the map 

(1.18) 5oo(A):Cr(C^)-C-(Q, 

(1.19) Soo{X)9 ^ v+\ci,, 

so it relates the incoming amplitude v- Ign-i to the outgoing one, w+lgn-i . We recall 
from [Q that the wave operator free-to-free S-matrix is then given by i"~^S'oo(A)i? 
(as maps C^(Co) -^ C^°°{C'q)) where R is pull back by the antipodal map on Cq. 
There are only minor chan ges if Va is poly homogeneous of order — 1 . Namely, 



the asymptotic expansions in ( 1.15 ) and ( 1.16 ) must be replaced by 

(1.20) 

g±»v^r^-»a±-(„-i)/2^^^ a± = a±,A = ±V^'|c'/2%/AeC°°(C^), V = xV', 



(1.21) v±^Y.T. a,,.,±(c.)r-^"(logr)^ a,, .._ G C^iC'o), a,, ..+ G C°°iC'o). 

j=0 s<2] 

Note that a± are not defined at Cp.sjn gi but that does not cause any problems 
even in the uniqueness statement, ( 1.15 ), since u_ vanishes at §"~^ near Co,sing to 



infinite order. 

Our main theorem describes the structure of S'oo(A). We first introduce the 
broken geodesic flow (of the standard Riemannian metric h) on §"^^, broken at 
C. We denote by S'S"^^ the sphere bundle of §"^^ identified as the unit-length 
subbundle of T§"^^ with respect to h. Let / = [a, /3] C M be an interval. We 
say that a curve 7 : / -^ §"^^ is a broken geodesic of ft,, broken at C, if two 
conditions are satisfied. First, there exists a finite set of points tj £ I, a — to < 
ti < ... < tk-i < tk — (3 such that for each j, 7|[(.^f .^^j is a geodesic of h, and 
for all t e {tj,tj+i), 7'(t) e S'S"-^ Second, for all j, if -/(tj) £ C"^ then the limits 
7'(tj — 0) and 7'(tj + 0) both exist and differ by a vector in T^^^jS"^^ which is 
orthogonal to T^i^t.^Ca (i.e. the usual law of reflection is satisfied; see Figure p. We 
say that p,q £ SS'^~'^ are related by the broken geodesic flow at time tt if there is a 
broken geodesic 7 defined on [0, tt], such that 7'(0) = p, ■y'in) — q. Using the metric 
h to identify SdX and S*dX, this defines the broken geodesic 'flow' at time tt on 
S*dX. We refer to Definition |6.6| and Section ^ for a more complete discussion. 
We then have the following result: 



6 ANDRAS VASY 



Theorem. Suppose that no subsystem of H has bound states, i.e. for a ^ 0, 
speCpp{ha) = 0. Then the free-to-free scattering matrix, S'oo(A), extends to a con- 
tinuous linear map C~°°{Cq) —^ C~°°{C'q). The wave front relation of Soq{X) is 
given by the broken geodesic flow at time n. 




Figure 1. Broken geodesies on S^ starting at p. Here Ca = Cbf] Cc. 

In the actual many-body problem, w G Xa means that several particles are 
close to each other, namely the ones corresponding to the cluster decomposition 
a. Thus, w e Ca is a statement that the particles corresponding to cluster a col- 
lide. Hence, the Theorem describes how many-body scattering can be understood, 
modulo smoothing (hence in the C°° sense trivial) terms, as a sequence of a finite 
number of collisions involving the particles. Namely, each 'break' tj in the broken 
geodesic describes a collision involving the cluster decomposition a. In the three- 
body setting with Schwartz potentials it was shown in [^ that the amplitude of 
the reflected wave is given, to top order, by the corresponding 2-body S-matrix; an 
analogous statement also holds for short-range potentials. In particular, this shows 
that the Theorem is sharp as far as the location of singularities is concerned. 

We also remark that in the Euclidean setting, unbroken geodesic flow to distance 
TT amounts to pull-back by the antipodal map on S"~^ = 9§", so it corresponds 
to free propagation: particles leave in the direction opposite to the one from which 
they entered. 

Our approach to proving this theorem is via the analysis of generalized eigen- 
functions of H, i.e. of u G C^°°(§") satisfying {H — X)u = 0. We prove that 
'singularities' of generalized eigenfunctions of H propagate along broken bicharac- 
teristics in the characteristic set of H, similarly to singularities of the solutions of 
the wave equation. Here 'singularities' are not understood as the lack of smooth- 
ness: indeed H is elliptic in the usual sense, so every generalized eigenfunction is 
C°° in the interior of E^, i.e. on R". Instead, in this situation singularity means the 
lack of rapid decay u. Correspondingly, we define a wave front set, WFsc(m), at in- 
finity, i.e. at 9S" , and we will prove its invariance under the broken bicharacteristic 
flow. 

The two notions of singularities are very closely related via the Fourier transform. 
Here for simplicity consider A — A in place oi H — X. If (A — X)u = 0, then the 
Fourier transform of u, Tu, satisfies (|^p — X)Tu = where ^ is the dual variable 
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of w. Now, the multiplication operator P — |^P — A can be regarded as a Oth order 
differential operator. Hence, by Hormander's theorem, see e.g. [|ll|, WF(JFm) is 
invariant under the bicharacteristic flow in the characteristic variety of P, i.e. in 
the set {(C,^*) : |^P — A = 0} where we have written ^* for the dual variable of 
£,, so S,* is in fact w. Moreover, in the two-body problem, i.e. if F is a symbol (of 
say order -1) on K", F = A + V, and if {H - X)u = 0, we still have PTu = 
where now P = |^p — A + TYT^^. Since V^ is a symbol of order —1, !FV!F^^ is 
a pseudo-differential operator of order —1, hence lower order than |^p — A. Thus, 
the principal symbol of P is still |^p — A (recall that ^* is the cotangent variable, 
so this is indeed homogeneous of order in ^* - it is independent of ^*). Hence, 
Hormander's theorem is applicable and we have the same propagation statement 
as before. 

In the two-body setting the relevant wave front set measuring lack of decay at 
infinity is the scattering one, WFgc- For u G 5'(M"), WFsc(u) is essentially given 
by the usual wave front set of the Fourier transform of u, i.e. by WF(J^m), after 
interchanging the role of the base and dual variables. Since the Fourier transform 
interchanges decay at infinity and smoothness, WF(J^u) indeed measures the decay 
of u at infinity in a microlocal sense. Hence, Hormander's propagation theorem 
translated directly into a propagation theorem for WFsc(u). This result was de- 
scribed by Melrose in |gl| where he introduced the notion of WFgc- 

In the many-body setting conjugation by the Fourier transform is much less 
convenient. Hence, we will design an appropriate microlocal way of measuring the 
lack of decay at infinity without resorting to the Fourier transform. Instead, we 
introduce an algebra of many-body pseudo-differential operators \['sc(§" ,C) which 
reflects the geometry, and use it to define the wave front set at infinity. We then 
prove a propagation of singularities theorem for generalized eigenfunctions of many- 
body Hamiltonians H; here 'singularities' are understood in the sense of the new 
wave front set at infinity. The proof of this theorem is via a microlocal positive 
commutator estimate, similarly to the proof of Hormander's theorem, or indeed to 
the proof of the propagation theorems for C°° singularities of solutions of the wave 
equation with domains with boundaries [E2l. Finally, we relate such a result to 
the structure of the S-matrix. This step is comparatively easy as indicated by our 
description of the S-matrix in terms of generalized eigenfunctions of H. 

Positive commutator estimates have also played a major role in many-body scat- 
tering starting with the work of Mourre |Q , Perry, Sigal and Simon ||22l , Froese 
and Herbst [g, Jensen |Q, Gerard, Isozaki and Skibsted [g[ 0| and Wang Q. In 
particular, the Mourre estimate is one of them; it estimates i[H, w ■ D^ -|- D^ ■ w]. 
This and some other global positive commutator results have been used to prove the 
global results mentioned in the first paragraph about some of the S-matrices with 
initial state in a two-cluster. They also give the basis for the existence, unique- 
ness and equivalence statements in our definition of the S-matrix by asymptotic 



expansions; these statements are discussed in [38 in more detail. Correspondingly, 
these global estimates will appear in Sections |ll-12 of this paper where we turn 
the propagation results for generalized eigenfunctions into statements about the 
S-matrix. 

We remark that the wave-operator approach defines the S-matrix as a bounded 
operator L?{Cq) -^ L^(Co). Since Co, sing has measure 0, L^(Co) and L?{Cq) can 
be identified. As C^{Cq) is dense in L^(Cq), the asymptotic expansion S-matrix 
S'oo(A) indeed determines the wave-operator one. 



8 ANDRAS VASY 

The propagation of singularities of generalized eigenfunctions of H is determined 
by the principal part of H; terms decaying at the boundary do not change the anal- 
ysis. As opposed to this, the precise structure of incoming and outgoing functions, 
R{X ± iO)/, / e C°°(§" ), depends on lower order terms; a relatively trivial exam- 
ple is given by the appearence of r^'"± in (1.20D for long-range potentials. Since 



we consider S'oo(A) and Po,-i-(-^) £^s operators on distributions supported away from 
C'o,sing, we do not need to analyze the precise structure of incoming/outgoing func- 
tions at Co, sing, which is not 'principal type', although we certainly analyze the 
propagation of singularities there. Thus, we do not discuss what happens when the 
support of the incoming scattering data increases to Cg, even if the data are L^. 
But the behavior of Po,-i-(A), as the support of the data increases to Cq, plays an 
important part in asymptotic completeness, which states that all possible outcomes 
of a scattering experiment are indeed described by a combination of bound states of 
the cluster Hamiltonians, with asymptotically free motion in the intercluster vari- 
ables. Thus, our results cannot be used directly to supply a proof of asymptotic 
completeness. This completeness property of many-body Hamiltonians was proved 
by Sigal and Soffer, Graf, Derezihski and Yafaev under different assumptions on 
the potentials and by different techniques [|^ g|, ^ |2^, |[ ||, 0. In particular, 
Yafaev's paper pd[ shows quite explicitly the importance of the special structure of 
the Euclidean Hamiltonian. This structure enables him to obtain a positive com- 
mutator estimate, which would not follow from our indicial operator arguments in 
Section || and which is then used to prove asymptotic completeness. 

Finally we comment on the requirement that the collection C be locally lin- 
earizable. We show in the next section that it is equivalent to the existence of 
a neighborhood of every point p G dX and a metric on it, in terms of which all 
elements of C are totally geodesic. The importance of this assumption is closely 
related to the existence of a sufficient number of smooth vector fields on dX which 
are tangent to every element of C. Such smooth vector fields always exists once 
we resolve the geometry of C, i.e. on the blown-up space [i9X;C], but in general, 
without our assumption, there are not enough such smooth vector fields on dX . 
In the first part of the paper, we discuss the pseudo-differential algebra associated 
to many-body scattering. For this purpose we need to blow up C, in part for an- 
alyzing the indicial operators (see the following paragraph). Thus, in this part of 
the paper, the issue of local linearizability is irrelevant, and we do not assume it. 
However, in the second part of the paper, both the discussion of generalized broken 
bicharacteristics and the construction of the positive commutators would be more 
complicated without it, so from Section |g on, we assume the local linearizability of 
C. 

This paper is organized as follows. In the next section we describe the geometric 
generalization of the many-body problem which was outlined above. This includes 
a discussion of many-body geometry and the definition of many-body differential 
operators. In Section ^ we proceed to define and analyze the corresponding algebra 
of pseudo-differential operators, "^sc{X,C), which reflects this geometry. It includes 
many-body Hamiltonians, as well as their resolvent away from the real axis. It 
extends the definition of the three-body calculus presented in |^ , though here we 
emphasize the definition of the calculus via localization and quantization as opposed 
to the conormal description of the kernels on an appropriate resolved space. In 
Section Q we construct the indicial operators in this calculus. They provide a non- 
commutative analog of the principal symbol in standard microlocal analysis. Our 
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proof of positivity in conimutator estimates is based on replacing the argument of 
Froese and Herbst [g| by indicial operator techniques. In Section ^ we define the 
wave front set at infinity, WFsc('«), corresponding to the many-body geometry and 
pseudo-differential operators. The definition given here differs from the one in [ p9[ ; 
it follows the fibred cusp definition of Mazzeo and Melrose [^SJ. These definitions, 
however, give the same result for approximate generalized eigenf unctions of H. 

In Section O we discuss many-body type Hamiltonians and their generalized 
broken bicharacteristics. This section is, to a significant degree, based on Lebeau's 
paper p^ . In Section ^ we give a much more detailed description of the generalized 
broken bicharacteristics in the case when all elements C & C are totally geodesic. 
Of course, this is true in the Euclidean setting. In Sections p|-0 we build the tech- 
nical tools for turning a symbolic positive commutator calculation into an operator 
estimate. In Section nO we prove that singularities of generalized eigenfunctions 
of many-body type Hamiltonians propagate along generalized broken bicharacter- 
istics. This is the main new result of the paper. In Sections |ll| - [l2| we use this 
and adaptations of the global estimates, in particular those of Gerard, Isozaki and 
Skibsted |^ , to analyze the structure of the resolvent and that of the scattering 
matrix. Finally, in the Appendix we prove some of the results quoted from Lebeau's 
paper, using slightly different methods. 

The propagation estimates of Section |l^ lie at the heart of this paper. The 
reader may want to skip some of the technical sections when reading the paper for 
the first time. It may be useful to keep Mourre-type estimates and especially their 
microlocalized versions as in p, R] in mind while reading Section HO. 
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grateful to Maciej Zworski for introducing me to the work of Gilles Lebeau [O, 
for many helpful discussions and for his encouragement. It was Lebeau's paper 
that convinced me that the results presented here were within reach, and it plays a 
particularly central role in Section where the generalized broken bicharacteristics 
are described. I would also like to thank Andrew Hassell, Rate Mazzeo, Erik Skib- 
sted and Jared Wunsch for helpful discussions, their encouragement and for their 
interest in this research. 

2. Many-body geometry and differential operators 

It is convenient to carry out the construction in the general geometric setting. 
We first describe the many-body geometry. 

Thus, let AT be a compact manifold with boundary, and let 

(2.1) C^{Ca: a e /} 

be a finite set of closed embedded submanifolds of dX such that dX = Cq G C 
and for all a,b £ I either Ca and C(, are disjoint, or they intersect cleanly and 
CaCiCb = Cc for some c £ I. We introduce a partial order on C given by inclusion 
on C, namely 

(2.2) Ca < Cb if and only if Ca C Cb- 
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This partial order is the opposite of the partial order used traditionally in many- 
body scattering, discussed in the introduction, but it will be more convenient for 
us since it simply corresponds to inclusion. A chain is defined as usual as a set on 
which < gives a linear order. 

Definition 2.1. Let X and C be as above. We say that {X,C) is a space with 
A^-body geometry (or an iV-body space), A > 2, if the maximal length of chains 
is A^ — f . Similarly, we say that C'a is a fc-cluster if the maximal length of chains 
whose maximal element is Ca is A; — f. We also say that (A, C) is a many-body 
space if we do not wish to specify A^. 

Thus, if Ca is minimal, it is a 2-cluster, and if (A", C) is a space with A^-body 
geometry then dX is an A-cluster. The numerology is chosen here so that we 
conform to the usual definitions in Euclidean many-body scattering, described in 
the Introduction. 

Before defining the algebra of many-body scattering differential operators on 
(A, C), we discuss the simultaneous local linearizability of the collection C. As we 
have mentioned in the Introduction, the analysis of generalized broken geodesies as 
well as the commutator constructions of this paper become simpler if C is locally 
linearizable. To make this notion precise, we make the following definition. 

Definition 2.2. We say that a many-body space {X,C) is locally linearizable (or 
is locally trivial) if for every p G dX there exists a diffeomorphism (p from a neigh- 
borhood C/ of p in dX to a neighborhood U' of the origin of a vector space V such 
that for each C G C, the image oi C O U under cj) is the intersection of a linear 
subspace of V with U' . 

Remark 2.3. In three-body type geometry, where the elements of C except Co are 
disjoint, (A, C) is automatically locally linearizable. The same holds, essentially by 
definition, if C is a normal collection, see ||l^, Chapter V]. 

Local triviality of C is closely related to the question whether every element of 
C is locally totally geodesic with respect to some metric. In fact, 

Lemma 2.4. A many-body space (A, C) is locally linearizable if and only if every 
p G dX has a neighborhood U in dX and a Riemannian metric hu on U such that 
for each element C ofC, C HU is totally geodesic with respect to hu. 

Proof. Suppose first that p E dX and f/, hu are as above. By shrinking U if 
necessary, we can make sure that p ^ C implies C n C/ = for every C E C . By 
shrinking U further if necessary, we can arrange that the exponential map of hu at 
p E dX identifies a neighborhood U' of the origin in ]/ = TpdX and U. Moreover, 
the elements C E C for which p E C, are identified with TpC n U', since these C 
are totally geodesic. This proves that {X,C) is locally linearizable. 

Conversely, if (A, C) is locally linearizable, then the choice of an inner product 
on V induces a metric on TV, hence on U via the diffeomorphism (jj, and as linear 
subspaces of V are totally geodesic with respect to this metric on TV, the same 
holds for C over U. D 

After this brief discussion on the local linearizability of C, we turn to the setting 
of most interest, namely to Euclidean many-body geometry. Suppose that A = S" 
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is the radial compactification of R" and A" is a family of linear subspaces of M" as 
discussed in the introduction. Recall from |gl| that RC : M" -^ §" is given by 

(2.3) RC{w)^il/{l + \w\y/\w/{l + \w\y^^)eElcW'+\ u.eR". 

Here we use the notation RC instead of SP, used in |£1[, to avoid confusion with 
the standard stereographic projection giving a one-point compactification of R". 
We write the coordinates on R" = Xa ® X" as {wa,'w°'). Let m — dimXa- We 
again let 

(2.4) Xa - cl(RC(X„)), Ca^XaD dS^. 

We next show that polyhomogeneous symbols on X"', pulled back to R" by tt'^, 
are smooth on the blown-up space [X;Ca]- Recall that the blow-up process is 
simply an invariant way of introducing polar coordinates about a submanifold. A 
full description appears in pS( | and a more concise one in [Ell Appendix A] , but we 
give a brief summary here. Thus, suppose that X is a manifold with corners and 
C is a p-submanifold (i.e. product submanifold) of dX. Thus, near any p G C we 
have local coordinates Xi (i — 1, . . . ,r), yj (j = 1, . . . ,n — r), n = dimX, such that 
the boundary hypersurfaces of X through p are defined by Xi = 0, and X is given 
by cci > 0, i = 1, . . . , r. A tangent vector V € TqX, q near p, is inward-pointing if 
Vxi{q) > for all i. The normal bundle of C is the quotient bundle 

(2.5) NC = TcX/TC. 

The inward pointing normal bundle of C, N~^C, is the image oiT^X, consisting of 
inward pointing tangent vectors, in NC. Thus, near p, X is diffeomorphic to the 
inward-pointing normal bundle of C. The blow-up of X along C is locally defined 
as the blow up of the section of N~^C, i.e. by introducing the new C°° structure 
in N~^C given by polar coordinates in the fibers of the bundle and by the base 
coordinates pulled back from C. While this construction depends on some choices, 
the resulting C°° structure does not. The blow-up of X along C is denoted by [X; C] . 
The blow-down map [A^; C] — > X is the smooth map corresponding to expressing 
standard coordinates on a vector space, N^C, in terms of polar coordinates. It 
is denoted by I3[X;C]. The front face of the blow-up is the inverse image of C 
(i.e. of the zero section of N~^C) under P[X;C]. Hence, it is a bundle over C 
whose fibers are the intersection of a sphere with a 'quadrant' corresponding to the 
inward-pointing condition, i.e. to Xi > 0. In fact, it is the inward pointing sphere 
bundle S^NC which is the quotient of N^C \ o, o denoting the zero section, by 
the natural M+ actions in its fibers. 

We again return to the Euclidean setting. In particular A = S" . We denote the 
blow-down map by (3[X; Co] ■ [A; Co] — *■ X. Now S^NCa is a hemisphere bundle 
over Ca, which can be identified with the radial compactification of the normal 
bundle of Ca in dX whose fibers can in turn be identified with A". To see this in 
more concrete terms, we proceed by finding local coordinates on [A; Ca] explicitly. 
It is convenient to do so by using projective coordinates rather than the standard 
polar coordinates. Near Ca in §" we have \wa\ > c|w"| for some c > 0. Hence, near 
any point p d Ca one of the coordinate functions {wa)j which we may take to be 
{wa)m, satisfies \{wa)m\ > c'\{wa)j\, \{wa)m\ > c'\w"-\ for somc c' > 0. Taking into 
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account the coordinate form of RC we see that 
(2.6) 

X^ \[Wa)m\ , Zj = — - J = 1,... ,TO-1), y, = — - {] = I, . . . , n - m) 

\(Wa),n\ \[Wa)7n\ 

give coordinates on §" near p. In these coordinates Ca is defined by a; = 0, y = 0. 
Correspondingly, we have coordinates 

(2.7) X, Zj (j = 1, . . . , m - 1), Yj = yj/x {j = 1, . . . ,n-m), 

i.e. 

(2.8) 

X = \{wa)m\~^, Zj = /"Y (j = 1, . . . , m - 1), y,- = (ii;")j (j = 1, . . . , n - to) 

I i,^a jm I 

near the interior of the front face fF of the blow-up [X; Ca], i.e. near the interior of 
ff = /3[X;Ca]*a; see Figure a 



/ 



/ 



/ 

A\y\< 



/ 



\ 



\ 



\ 



\ 



\ 



X 

^ = const 



Y\= m <c 




Y = ^ = const 



Figure 2. The blowup of Ca — {x — 0, y = 0}; the z coordinates 
are normal to the page and are not shown. The thin lines are the 
coordinate curves Y — const and x ~ const in the region |y| < c 
(which is disjoint from /3[X;Ca]*dX), and their images under the 
blow-down map f3[X; Ca]- 

Near the corner d(3[X;Ca]*Ca = (3[X;Ca]*Ca n (3[X;Ca]*dX, in the lift of the 
region defined for some k by \yk\ > c\yj\ for some c > and all j ^ k, 

(2.9) X = x/yk, Yj = yj/yk {j ^ k), yk, z 

give coordinates. In terms of the original Euclidean variables these are 

X = \{w'^)kV\ Z, = ^^ (j- = 1, . . . , TO - 1) 

(2.10) 



Y, 



\{Wa)„ 



(j = 1,... ,n-m, 2^ k), yk 



(w'^)fe 



{w'^)k '■' ^'•■- '■- ■■-' -^ ' ■■^" ""■ \{wa)m\' 

Since in every region near the lift P[X;Ca]*Ca of Ca we can use one of these 
coordinate systems, and since away from there we can use coordinates as in ( ^.6[ ) 
but with Wa and w" interchanged, we have proved the following lemma. 
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Lemma 2.5. Suppose that X = §" and let j3 = I3[X] Ca] he the blow-down map. 
Then the pull-back /3*(RC"^)*7r" o/tt'' : M" ^ X" extends to a C°^ map, which we 
also denote by ■k"' , 

(2.11) T:'':[X-Ca\^X\ 

Moreover, ifx"" is a boundary defining function on X"" (e.g. x"" — \w"'\^^ for \w"'\ > 
1), then Pox = {'n°')*x^ is a defining function for the lift of dX to [X\ Ca], i.e. for 
(i*dX. 

Corollary 2.6. Suppose that X = §^, / e 5';^g(X''). Then 

(2.12) (7rT/ePaiC°°([^;a]). 

Here, following the previous lemma, we regard vr" as the map in ( |2.1l| ), and pgx is 
the defining function of (3[X;Ca\*dX , i.e. of the lift ofdX, and the subscript phg 
refers to classical (one-step polyhomogeneous) symbols (see Figure [^. 




(]*dS 




Figure 3. The blowup of Ca in §+; P = 0[^l;Ca] is the blow- 
down map and ff = 13* C a- X'^ and X'^ denote translates of Xa in 
W\ X'^ = cl(RC(X^)), etc. Note that the lifts of Xa, X'a and X'^ 



become disjoint on 



s^; Ca\ 



This corollary shows that for a Euclidean many-body Hamiltonian, H = A + 
J2a ^a, Va bccomcs & nicc function on the compact resolved space [S" ; Ca]. Thus, 
to understand H, we need to blow up all the Ca. In order to analyze this iterated 
blow-up procedure, it is convenient to generalize the clean intersection properties 
to manifolds with corners X . 

Let X be a manifold with corners. Given a finite family J- of closed p-submanifolds 
Fi of X we say that JF is a cleanly intersecting family if it is closed under intersec- 
tion (in the sense that any two members are either disjoint, or their intersection is 
in the family) and for any i and j, {Fi, Fj} form a normal collection in the sense of 
Melrose p^, Chapter V]. Thus, for any point p e F,; n Fj there are local coordinates 
on a neighborhood U oi p such that with some index sets /', /", 

(2.13) F,nU = {xr = 0, re I'„ y, = 0, se II'}, 

(2.14) FjHU^ {xr ^0, rei;, ys^O, se I';}; 

here the Xk are defining functions of the boundary hypersurfaces through p. This 
simply means that there is a common product decomposition for any pair of ele- 
ments of J-. In particular, if X is a manifold without boundary, then this simply 
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means that the Fi pairwise intersect cleanly. Hence, {X^C) is a many-body space 
if and only if C is a cleanly intersecting family in dX which includes dX. 

Just as in the case of a many-body space, inclusions gives a partial order on a 
cleanly intersecting family !F . Thus, F £ J- is minimal with respect to inclusion if 
there is no F' Q T such that F' y^ F, F' C F. Since ^ is closed under intersection, 
this means exactly that for all F' E T cither F' and F are disjoint, or F C F' . 

Lemma 2.7. Let T he a cleanly intersecting family of p-submanifolds of dX. Sup- 
pose that F Cz J-' is minimal with respect to inclusion. Then the lifted family, T' , 
consisting of the lifts of Fj, distinct from F, to [X;F], is also a cleanly intersecting 
family. 

Proof We claim that for any F„ Fj G .F the 4-tuple {F, Fk,F„Fj}, Fk = F, n Fj, 
is a normal collection in the sense of Melrose. Indeed, this is clear if Fk is disjoint 
from F; otherwise F C -Ffe by our assumption. 

So assume that F C Fk. By the normality of {F,Fk}, near any point p in F 
there are local coordinates Xr, j/s, on X such that 

(2.15) Fk = {x, == 0, r e 4, ys = 0, se 4'}, 

(2.16) F = {xr^O, re /', ys =0, se /"}, 

and /[, C /', I'k C /". Similarly, by the normality of {Fi,Fj} there are local 
coordinates a;^, y'^ near p on X such that 

(2.17) F, = {4 = 0, r e li y: = 0, se in, 

(2.18) F, = K = 0, re /;, y's^O, se /j'}. 
Thus, 

(2.19) Fk = {< - 0, r e /^ U /;- , yi = 0, s e /f U /j'}. 

Thus, the differentials of the coordinates x'^, r e I'i'Jlj, and y'^, s e /"U/", span the 
conormal bundle of F^. The same holds for the differentials oi Xr, r e I'k, ys, s e I'l- 
It follows that the differentials of a;^, r G ^ U Ij, Xr, r ^ I'^, y'^, s e /■' U /", ys, 
s ^ I'^ are independent at Fk in a coordinate neighborhood of p, so these functions 
give local coordinates on X near p in terms of which F, Fk, Fi and Fj have common 



product decomposition: F^, Fj and Fj. given by (2.17)-(2.19), and F by 
F = {xi = 0, r G /^ U /,', x^ = 0, r G /' \ I'k, 

(2.20) l r- , ^ J' »■ ' \ fc' 

y, = 0,seI,Ulj,ys^O,seI \ Ik }■ 

This proves that {F, F^, Fi, Fj} is indeed a normal collection. Hence, by [|l9|. Lemma 
V.11.2], it lifts to a normal collection of p-submanifolds on [X;F]. Writing /3 for 
the blow-down map, and P*Fk for the lift of Fk, etc., we see in particular that 
{/3*Fi, f3*Fj} is a normal collection whose intersection is P* Fk if Fk ^ F, and is 
empty otherwise. Putting together these facts we see that we have proved the 
lemma. D 

This lemma allows us to define [X; T\ if J^ is a cleanly intersecting family of p- 
submanifolds of dX. We do this by putting a total order on T which is compatible 
with the partial order given by inclusion. This can always be accomplished: pick a 
minimal element with respect to inclusion, and make it the minimal element of the 
total order. Proceeding inductively, if we already placed a total order on T' C T , 
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we choose any F Cz T\!F' which is minimal with respect to inclusion in JF \ JF' , and 
extend the total order to !F' U {F} by making F the maximal element with respect 
to it. Having imposed a total order on JF which is compatible with inclusion, we 
define [X; T] to be the blow up [X; Fi, F2, . . . , F„] where !F = {i^i, F2, . . . , F„} and 
Fi < F2 < ■ ■ ■ < Fn, < being the total order. Of course, a priori [X;^] depends 
on the total order. The following lemma shows that this is not the case. 

Lemma 2.8. If J-' is a cleanly intersecting family and <, <' are total orders on it 
which are compatible with inclusion, then the blow ups 

(2.21) [X;Fi,F2,... ,F„l Fi < F2 < . . . < F„, 

(2.22) [X;F[,F^,...,F:,], F[ <' F^ <'...<' FI, 

are canonically diffeomorphic. 

Proof. Since any total order compatible with inclusion can be obtained from any 
other one by repeatedly interchanging the order of adjacent elements, but keeping 
the order compatible with inclusion, it suffices to show that 

(2.23) [X-Fi,... ,Ffc,Ffe+i,... ,F„] and [X;Fi,... ,Fk+uFk,... , F„] 

are naturally isomorphic if both of these total orders respect inclusion. Now, cither 
Fk n Fk+i — 0, in which case the statement is clearly true, or Fk n F^+i = Fj for 
some j . Since inclusion is respected, we must have j < k. But upon the blow up of 
their intersection, any two closed p-submanifolds with normal intersection lift to be 
disjoint. Hence, on [X;Fi,... ,Fk-i] the lifts P*Fk and f3*Fk-i-i are disjoint, and 
thus they can be blown up in cither order. This proves the lemma. D 

Correspondingly, [X;.F] is defined independently of the total order used in the 
definition of the blown up space, assuming that it respects inclusion, so we can 
speak about [X;.F] without specifying such a total order. 

If Fi E !F, we can always specify the total order so that every Fj e !F with 
Fj < Fi satisfies Fj C Fi. Then the blow-up of Fi commutes with all the ones 
preceeding it. Hence, any function that is smooth on [X; Fi] pulls back to be 
smooth on [X;JF]. Applying this in the Euclidean many-body setting we conclude 
that 

Lemma 2.9. Suppose that X = S" and X is a linear family of subspaces of R" 
as in the introduction. Then V = J^a ^o-' ^^ ^ ^vha^'^'^^' ^^A^ ^^ ^^ ^^ element of 
Pg^C°°([X;C]) where pgx is the defining function of the lift of Cq = dX under the 
blow-down map 

(2.24) Psc^P[X-A--[X-C]^X. 

Our main interest is the study of differential operators, in particular the analysis 
of many-body Hamiltonians H . For this purpose we next investigate how vector 
fields lift under the blow up. First, we define Vh{X\ T) as the Lie algebra of smooth 
vector fields on X which are tangent to the boundary faces of X and to each element 
oiT. 

Lemma 2.10. Each element ofVb(X;!F) lifts to an element of Vb([X;!F]). 
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Proof. It suffices to show that V e Vh{X;J^) Hfts to be an element of Vb([^; F];^') 
where F is minimal with respect to inclusion and 

(2.25) T' ^{P*F' : F' eT\{F}}. 

Taking into account that for any F' ^ F, \F , F'} is a normal collection of p- 
submanifolds of X, this claim follows from [^ Proposition V.11.1], or it can be 
checked directly by using projective coordinates on [X; F]. D 

Remark 2.11. It is not the case in general that V\^{X;!F) lifts to span Vb([X;.F]) 
over C°°{[X;T]). This statement is true, however, if J^ is a normal collection (i.e. 
all elements of J^ have product decomposition in the same coordinate system, not 
just pairs of elements), see 11^, Proposition V.11.1]. 

We can now introduce the appropriate class of differential and pseudo-differential 
operators on many-body spaces {X, C). These will include many-body Hamiltonians 
in the Euclidean setting as well as their resolvents (in the resolvent set). 

First, we recall from |gl|] Melrose's definition of the Lie algebra of 'scattering 
vector fields' Vsc(-'f), defined for every manifold with boundary X. Thus, 

(2.26) Vsc{X) ^ xVb{X) 

where Vb(X) is the set of smooth vector fields on X which are tangent to dX. If 
{x,yi, . . . ,y„_i) are coordinates on X where x is a boundary defining function, 
then locally a basis of Vsc(^) is given by 

(2.27) x^d^, xdy^, j = l,... ,n-l. 

Correspondingly, there is a vector bundle ^'^TX over X, called the scattering tangent 
bundle of X, such that Vsc(-'^) is the set of all smooth sections of ^'^TX: 

(2.28) Vsc(X) = C°°(X; """TX). 

The dual bundle of ^^TX (called the scattering cotangent bundle) is denoted by 
^^T*X. Thus, covectors v G ^^T*X, p near dX, can be written a,sv = TJ^ + fi-^. 
Hence, we have local coordinates (x, y, t, fi) on ^'^T* X near dX. The scattering 
density bundle ^'^Q.X is the density bundle associated to ^'^T*X, so locally near 
dX it is spanned by a;^"^^ dxdy over C°°{X). Finally, Diffsc(^) is the algebra of 
differential operators generated by the vector fields in Vsc(^); E)ifT^(-'^) stands for 
scattering differential operators of order (at most) m. 

To establish the relationship between the scattering structure and the Euclidean 
scattering theory, we introduce local coordinates on X near p G dX as above, 
and use these to identify the coordinate neighborhood U oi p with a coordinate 
patch U' on the closed upper hemisphere §" (which is just a closed ball) near 
its boundary. Such an identification preserves the scattering structure since this 
structure is completely natural. We further identify S" with R" via the radial 
compactification RC as in ( |2.3| ). The constant coefficent vector fields d^^ on R" 
lift under RC to give a basis of '''^TS" . Thus, V E Vsc(S" ) can be expressed as 
(ignoring the lifting in the notation) 

TL 

(2.29) V^ = ^a,a„^., ajEC^iSl). 

As mentioned in the introduction, aj G C°°(§") is equivalent to requiring that 
RC* ttj is a classical (i.e. one-step polyhomogeneous) symbol of order on R". 
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This description also shows that the positive Euchdean Laplacian, A, is an element 
of DifFg(,(§" ), and that '"^fiS" is spanned by the pull-back of the standard Euclidean 
density \dw\. 

If X is a manifold with boundary then any element of Vsc(^) — xVh{X) is 
automatically tangent to any submanifold C of dX. Hence, due to Lemma ^.10| , we 
can define the algebra of many-body differential operators as shown by the following 
proposition. 

Proposition 2.12. If {X, C) is a many-body space, then Vsc(X) lifts to a suhalgebra 
ofVh{[X;C]). Correspondingly, 

(2.30) Diffs,(X,C) =C°°([X;C])(»c-(x)Difr,,(X) 

is an algebra. 



Proof By the first part of the statement, for any V E Vsc(-'^), / G C°°([X;C]), the 
commutator [ V^, /] = F/ is in C°° ( [X ; C] ) . D 

Since in Euclidean many-body scattering A G Diffgj,(S![:) and V = X^a ^a ^ 
C°°([§!;i;C]), it follows immediately that H = A + V e Diff^^(§!;:,C). 

3. Many-body pseudo-differential calculus 

Let {X;C) be a many-body space, and /3sc '■ i^'iC] -^ X the blow-down map. 
There are two equivalent way of defining many-body pseudo-differential operators. 
We can either specify their kernels as conormal distributions on an appropriately 
resolved space, or we can define them as the quantization of certain symbols. Here 
we give both definitions and show their equivalence. 

The b-double space, X^, has been defined by Melrose as [X'^; (dX)'^]. The front 
face of the blow-up is called the b-front face and is denoted by bf , while the lifts of 
the left and right boundary hypersurfaces oi X^, i.e. of dX x X and X x dX are 
denoted by If and rf respectively. The diagonal A of X^ lifts to a p-submanifold Ab 
of X^ which intersects dX^ in the interior of the b-front face, bf . (The definition 
of p-submanifolds and the blow-up process were discussed at the beginning of the 
previous section.) Moreover, Ab is naturally diffeomorphic to X. Hence, C can be 
regarded as a collection C of submanifolds of Ab, and, since Ab is a p-submanifold 
of X^, these submanifolds form a cleanly intersecting family in X^. Therefore, the 
blow up 

(3.1) Xl = [X^;C'] 

is well-defined by our previous results. Note that dX e C by our assumption, so 
the definition includes the blow up of the lift of 9Ab. It is easy to see that this 
space coincides with the Xg^ defined in |g^ if {X, C) is a 3-body space. 

Noting that even C U {Ab} is a cleanly intersecting family, we conclude that 
Ab lifts to a a p-submanifold, Asc, of X^^. Correspondingly, we define the set of 
many-body pseudo-differential operators by 

(3.2) %'/,{X,C) = {k e y^'"''(X|„ Ase; ""'ni^) : K = 0at/3* bf U/3* If U/3* rf}; 

here ^'^^r is the pull-back of the scattering density bundle from the right factor 
and P : Xg^ -^ X^ is the blow-down map. Similarly we define the corresponding 
polyhomogeneous operators 

(3.3) *™/(X,C) = {Kep'/™(X|„Asc;''f7fl,) : k = at /3* bf U/3* If U/3* rf } 
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where p is the total boundary defining function of X"^^. In particular, conornial 
distributions of order — cx) are smooth functions, so 

(3.4) ^g;^^'(X,C) = {«:ep'C°°(X2^,Asc;''f^K) : k ee at /?* bf U/?* If U/?* rf}, 

i.e. the kernels of operators in \E'gJ^' {X,C) are smooth up to all boundary hyper- 
surfaces of Xgj, (at least if Z is a non- negative integer), and vanish to infinite order 
at the lift of evey boundary hypersurface of X"^. Tensoring with vector bundles 
defines *";^'(X,C;£;,F) and ml';\x,C\E,F) for vector bundles E and F over X 
as usual. 

Since for all F E C we have F C 9Ab, we can do the blow up of 9Ab G C 
first, before blowing up other elements of C (normally we would do this blow up 
last by our total order construction). It follows that Xg^ is a blow up of the space 
X^^ = [X^;dA\,]. Hence, conormal distributions on X^^ pull back to be conormal 
on Xg^. Since the kernels of scattering pseudo-differential operators are conormal 
to Age and to the boundary of X^^ with infinite order vanishing at every boundary 
face except the scattering front face, we conclude that these kernels pull back to 



X|^ to be elements of the kernel space defined in (3.2), so ^^cc (-''^) C *"^,:(X,C) 



Suppose now that X = §" and C is a cleanly intersecting family of submanifolds 
of dX — (9§" — S"^^. Here we do not assume that C arises from a family X of 
linear subspaces of W\ An equivalent definition of ^scc(S+iC) is the following. 
Suppose that 

(3.5) aeA-"'-\[Sl;C]xSl). 

That is, identifying int(§") and int([S" ;C]) with R" as usual (via RC~^), suppose 
that a G C°°(R^ x M^) has the following property. For every P G Diff^(S!;:), acting 

on the second factor of S" (i.e. in the S, variable), and Q G Diff^ ([S";C]), acting 
on the first factor of §" (i.e. in the w variable), fc, k' G N, 

(3.6) PQaep^y'gL^^iSlxSl) 

where poo and pg are defining functions of the first and second factors of §" respec- 
tively. Let A = qL{a) denote the left quantization of a: 

(3.7) Auiw) = (27r)-" / e'^'"-'"'>^a{w, O^iw') dw' d^, 

understood as an oscillatory integral. Then A G ^g"i,(§",C). Indeed, the kernel of 
A is 

(3.8) K{w,w') — a{w,w — w') 

where 5, is the inverse Fourier transform of a in the ^ variable, i.e. a = 1F7 a. Thus, 
a(w, W) is smooth away from ly = 0, is conormal to T4^ = 0, and it is rapidly 
decreasing with all derivatives in W . More precisely, the rapid decay means that 
for aU /c and Q G Diffb([§'^;C]) and for all a, 

(3.9) sup {\w\'\W\''\Q^D'^~a(w,W)\)<^. 

\W\>\, luSR" 

Taking into account the geometry of ^1^., in particular that |w — w'j^^ vanishes at 



all faces of the blow-up (3.1) but the front faces (i.e. it vanishes at /?* If, /?* rf and 



/3* bf), we see that K vanishes to infinite order at these faces. Similar arguments 
describe the behavior of K near Age, proving that A G ^gJ;c(S+,C). 
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Conversely, ii A G ^n^'^(§",C) then there exists a satisfying (3.6) such that 



A — qhio). Namely, we let a(w, TV) = K{w,w — W) and let a be the Fourier 
transform of a in W . The conormal estimates for K (hence for a) give the symbolic 
estimates ( |3.6| ) for a. 

Similar conclusions hold for the right quantization B = qu{h) of a symbol h: 

(3.10) Bu{w) = (27r)-" / e*("'-™')-«6(u;', Ou{w') dw' d^ 

In addition, the polyhomogeneous class ^gj,' (§" ,C) is given by the quantization of 
symbols 

(3.11) aep^y^oC^{[§l;C]x^l). 

Since differential operators J2 o,a{w)D°' are just the left quantization of the symbols 
a{w,S,) = J2 '^aiw)£,°' , it follows immediately that 

(3.12) Dirs",(X,C)c*g^,(X,C). 

This conclusion also follows directly from the description of the kernels since the 
kernel of a differential operator is a differentiated delta-distribution associated to 
the diagonal. 

Note that, as usual, one can allow symbols a depending on w, w' and ^, so e.g. if 
a •= P^Pd,LPd,RC°°(i^+'^C] X [S+;C] X §"), pQ^L and pg^R denoting total boundary 
defining functions of the first and second factor of [§";C] respectively (i.e. they are 
pull-backs of a boundary defining function of S" ), then 

(3.13) Au{w) = (27r)-" / e*("'-^"'>«a(w, w',^ w(w') dw' d^ 



defines an operator A € ^^' ~^ (§",C). 

This characterization allows the application of the standard tools of the theory of 
pseudo-differential operators. In particular, if A g ^^' {X, C) is written as the left 

quantization of a symbol a and B G \E'^ ' (X, C) is written as the right quantization 
of a symbol 6, so 

(3.14) a e p^"^p'gC°^i[^l;C] x S'^), b e p^"^' p^^C^^ {[^l; C] x S'^), 
then the operator AB is given by 

(3.15) ABu{w) = {2tt)-" f e^^'"-'"">-^a{w, C)6(w', "("'') dw' d^. 



Here c{w,w',0 = a{w,Ob{w',0 is in P^'''-"'>^,iP^,flC-([§!^;C] x [S;;C] x S^), 
so we conclude that AB e ^gc ' (-^tC). In addition, the adjoint A* of A is 
the right quantization of a, so A* e ^g"' {X,C). Analogously, ^scc(S+,C) is also 
closed under composition and adjoints. These statements can be seen also from 
the standard more explicit formulae. For example, if B is the left quantization of 
a symbol b' , the composition formula, including the remainder terms, only involves 
derivatives of the form D'^DZb', and D^ € Diffl"l(S!^) C Diff{"'([§'^;C]), so we see 
that \l/scc(§+,C) is closed under composition. 

This discussion can be carried over to arbitrary manifolds with boundary X by 
locally identifying X with §" and using that our arguments are local in §" . More 
precisely, suppose that {C/i, . . . ,Uk} is an open cover of X by coordinate patches, 
and identify each Ui with a coordinate patch [// of §" . We write (pi : Ui -^ U[ for the 
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identification. Let C'^ denote the family given by the image of elements of C in U^. 
Then A e ^^/{X,C) if and only if there exists operators A'^ G ^^/{E>l;C') with 
kernel supported in the inverse image of U- x U^ in (§" )sc f^nd R € C°° {XxX; ^'^^b.) 
such that 

(3.16) A = ^(0:A:;(<^r')*)+^- 

i 

Note that the support condition on A^ ensures that this expression makes sense. To 
see this, just introduce a partition of unity pi G C°°{X) subordinate to the cover, 
and let ipi G C°°{X) be identically 1 in a neighborhood of supppi. Then 

(3.17) A^J2^P^^J2 ^'^^P' + 11(1 - ^i)Ap,. 

i i i 

It is straightforward to check directly from the definition of ^11^' {X,C) that the 
last terms is given by a kernel in C°°{X x X]^'^Ur), while A[ = {(j)Y^)*i^iApi(j)* g 
*Sc (^+i^i) ^itti the claimed support properties. Thus, our results for \1/^' (S" ,C) 
immediately show the following theorem. 

Theorem 3.1. Both "^sciXjC) and "^ scc{X , C) are ^-algebras (with respect to com- 
position and taking adjoints). 

Since ^^^^ (§!}:, C) C *S(I*"), where *™(K") is the class of pseudo-differential 
operators defined by Hormander |0, Section 18.1], arising by a quantization of 
symbols a G C°^\R" x R") satisfying 

(3.18) \DZDlaiw,0\<CMO"'-^^K 
and 

(3.19) *S(1^") ■■ {wy'H''{W") -^ (u))-''i7'^-™(M"), 

we immediately deduce the boundedness of elements of \l/™^ {X, C) between the 
appropriate weighted Sobolev spaces. 

Theorem 3.2. If A e ^^ci^^C) then A : iJJ'"(X) -^ i?J-™'"+'(X) is bounded. 

There is another way of characterizing the calculus ^^i,(§!Ji,C) via Hormander's 
Weyl calculus (see [|ll|. Section 18.5]). We describe it briefly here, only considering 
the Euclidean setting where the Ca arise from linear subspaces Xa', it is straight- 
forward to check that it agrees with the definition we have given above in terms of 



quantization of symbols as in (3.5). Namely, 'i!'^^ °°(S" ) is just the calculus on R" 
arising from the metric 

{wy {Q^ 

Similarly, if we take C to consist of a single element Ca, a 7^ 0, and if {wa, w") is 
the usual splitting of the coordinates, then $^^~°°(§" ,C') arises from the metric 

,,2 M„,,a^2 Je2 



(3-21) 9' - ,^.,, - ,^..a,, . ,,,,■ 



(,) ^ du^ (dw-y de 
{wy {w'^y (0^ 

In the three-body problem, CaHCb = if a, 6 7^ 0, we define the metric by localizing 
the g(°\ i.e. we consider a partition of unity 0a G C°°(S" ), a G /, supp^a n Ct = 
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unless 6 = 0, and define the metric 

(3.22) a = Y.^-9^"^- 

a 

(Here the (f>a are pulled back to the cotangent bundle by the bundle projection.) 
Since the g^°'' are equivalent near C'q, it follows that g is indeed slowly varying. 
Note that if (f>a is supported close to Ca, which we can arrange by enlarging the 
support of (/>o, dw'^/{w)'^ above can be replaced by dwi/ {wa)"^ ■ 

In general we simply repeat this procedure. Thus, to define the appropriate 
metric on T*X'^ if it has been defined on T*X° for every a with X" C X^, we 
define a partition of unity (j)a G C°°{X^) with supp^a n C^ = unless C^ C C^. 
Here X^ = X'^CBXl and C^ = ^X''r^c\{X''^)). We extend the metric g" on T*X'' to a 
symmetric 2-cotensor on T*X'^ using the orthogonal decomposition X'^ = X°'Q)X^, 
and let 



(3.23) g^'^^^g 



{dwlf , [d^lf 



c\2 l(:c\2 



(K)' (e 



Then 

(3.24) g-= Y. ^-9 



(a) 



a-.X-^CX" 



gives the desired metric on T*X'^. 

After this brief discussion of the relationship of ^g^*(S",C) with Hormander's 
Weyl calculus, we return to the general setting to describe the principal symbol 
map and its analog at dX. 

4. The principal symbol and the indicial operators 

Since the inclusion of Hl^'^ {X) to Hl^^{X) is compact for r' > r, s' > s, 
it suffices to understand A € \l/^'j.(X,C) modulo *^^ ' ^ iX,C) to analyze its 
spectral properties. Now, Hormander's principal symbol map on ^^"(M") restricts 
to a principal symbol map 

(4.1) asc,™ : ^sfi^+^C) ^ S^ ^^T* [SI; C]) , 

Sh{^'^T*[S^^;C]) denoting the space of smooth symbols which are homogeneous of 
degree m. Due to its invariance and its local nature, it immediately extends to a 
map 

(4.2) ascm : *™/(^,C) ^ SJ^e-T*[X;C]). 

We radially compactify the fibers of ^''T*[X;C] (i.e. replace the vector spaces by 
balls) and let ^"^S* [X; C] be the new boundary face (i.e. the boundary of ^^T* [X;C] 
at fiber- infinity). This allows us to write ffscm as a map 

(4.3) asc.„ : ^s"'°(^, C) ^ C^{Scg* [^. ^j . (^^Sc^* [^. ^^ymy 

The hue bundle N*^''S*[X;C] is locally spanned by the pull-back of d(|f|-^) from 
^'^T*[X;C], so (|4.3|) is obtained from (4.1) by writing homogeneous functions a{w,£,) 



of degree m as ao{w, ^)|^|'", ^ = C/|^|, considering ao as a function on the cosphere 
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bundle, and using N*^'^S*[X;C] to take care of the factor |^|'" invariantly. We then 
have a short exact sequence 

(4.4) 

^ *™-^'°(X,C) ^ *sc°(^.C) ^ C°°i^''S*[X;C]; {N*^'' S* [X ; C])-"") -^ 

as usual. 

Now, an operator A e vj/™' {X,C) is certainly determined modulo V?^' i^iC) 
by the restriction of its kernel to the front faces of the blow up ( |3.l| ). Our next 
task is to construct a multiplicative indicial operator from this restriction. We use 
oscillatory testing for this purpose as was done in [g9|. We start by discussing the 
effect of conjugation of A by oscillatory functions. 

Lemma 4.1. Suppose that A e *"/(X,C) and f G C°°{X;R). Then 

(4.5) A = e-'f'/^'Ae'f'/'' e *™/(X,C). 

Proof. It is convenient to use the explicit description of ^sc{X,C) in terms of 
localization and quantization (3.7). Thus, we may assume that X = §". Note 



that the pull-back of f/x to M" is a polyhomogeneous symbol of order 1 which we 
denote by F. Then the kernel of A is K{w, w') = e*(^("'')--^("'))X(w, w') where K 
is the kernel of A. But by the fundamental theorem of calculus 

(4.6) F{w') - F{w) = V« -Wj) / djF{w + t{w' - w)) dt, 



3 = 1 '" 

and djF is a polyhomogeneous symbol of order 0. Taking into account the rapid de- 
cay of ii" in PF = w-w' we immediately conclude that K e A'^-^iS'Dl^, Asc; KD|J 
vanishing with all derivatives at /3* bf U/3* If U/3* rf, so, returning to the global set- 
ting, ievl/-/(X,C). D 

We next discuss mapping properties on C°°{[X;C]). 

Lemma 4.2. If A e *"j'(X,C), u £ x'' C°° {[X ; C]) , then Au £ x''+^C°°{[X;C]). 

Proof. This result essentially reduces to the fact that ^sc(X, C) is an algebra. In- 
deed, write u — u ■ 1, and note that Au — {AU)1 where B = AU denotes the 
composite of A with the multiplication operator U by u. Since the latter is in 
x'^Diffg^(X,C), hence in ■^l'^{X,C), we conclude that B £ ^g'/+''{X,C). Thus, 
we only have to analyze Bl. Again, we can reduce the discussion to a local one. 



But writing B as the left quantization of a symbol b(w,^) as in (3.7), b satisfying 



(3.11) with I replaced hy I + r, and writing the oscillatory integral explicitly as a 
convergent integral, we see that 

(4.7) 

Bl{w) = (27r)-" f £■'("'-"'')•« (w - w')-^''(0"^'(l + ^iYHw, 0(1 + ^w'Tl dw' d^ 

for 2r > n, 2s > n + m. Changing the variables: 

(4.8) Bl{w) = {2Try f e'^-^iWy^^O'^'C^ + ^iyKw,OdW dC 



This is a convergent integral with w dependence only in b. Since 

(4.9) b£p^y'+^C^i[El;C]xSl), 
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we conclude that Bl e a;'+''C°°([§" ;C]). Hence, returning to the global setting, 
Auex^+''C°"{[X;C]) as claimed. D 

The previous two lemmas show that if u = e'f^''v, v e C'^{[X;C]), A e *"/ {X, C) 
then Au = e^^'^v' with v' e C°°([X;C]). Moreover, v' restricted to the boundary 
of [^;C] only depends on the restriction of v to i9[X;C]. It also follows from the 
above proof that if p € dX and v S C°°([X;C]) vanishes at P^^{p) then v' also 
vanishes there, i.e. composition is local in X (though not in the resolved space 
[X;C]). Similarly, if f{p) = f'{p) and dyf{p) = dyf'{p) (which really just mean 
that the scattering covectors d{f/x) and d{f' /x) agree at p) then e~^^^^Ae^^^'^v 
and e~^-^ ^^Ae^^ ^^v agree at p. This allows us to define the indicial operators of A 
at the boundary hypersurfaces of [X;C]. These depend on certain choices in gen- 
eral (though the dependence is via unitary equivalence), but if we have a scattering 
metric on X they can be constructed canonically, so we assume this in what follows. 

Recall first that a scattering metric 5 on X is a metric in the interior of X 
(smooth symmetric positive definite 2-cotensor) which is of the form 

r. .r.^ dx'^ h' 

(4.10) g=^ + — 

X^ X'' 

near dX, where a; is a boundary defining function oi X and h' is a smooth symmetric 
2-cotensor on X whose restriction to the boundary, h, is positive definite. Thus, 
g gives a positive definite pairing on ^'^TX, so it is (a somewhat special) smooth 
section of ^'^T*X (g) ^'^T*X. We remark that the choice of such a g fixes x up to the 
addition of functions in x'^C°°{X). 

Next, we recall the definition of the relative scattering tangent bundle ^'^T{C; X) 
of a closed embedded submanifold C of dX from |^ . 

Definition 4.3. For a closed embedded submanifold C of dX, the relative scat- 
tering tangent bundle ^'^T{C; X) of C in X is the subbundle of ^'^TcX consisting of 
V G ^'^TpX, p G C, for which there exists 

(4.11) VeVsciX;C)(lVsc{X) 
with Vp — V. Here 

(4.12) Vsc(X; C) == xVb{X; C) = x{V G Vb{X) : V is tangent to C} 
and tangency is defined using the (non-injective) inclusion map ^TX — > TX. 

Thus, in local coordinates (x, y, z) near p G C such that C is defined by a; = 0, 
y = 0, ^'^T{C; X) is spanned by x^dx and xdz , j — 1, . . . ,m—l where n — to is the 
codimension on C in dX. In the case of Euclidean scattering, X = §" , C — dXa, 
g the Euclidean metric, ^'^T{C;X) is naturally isomorphic to ^"^TcXa, i.e. it should 
be regarded as the bundle of scattering tangent vectors of the collision plane at 
infinity, spanned by dt^u^ . , j = 1, . . . , to,, to, = dim Xa- 

For C — Ca £ C, the metric g defines the orthocomplement {^'^T{C;X))-^ of 
"'=r(C;X) in'^'TcX. 

Definition 4.4. Given g, a scattering metric on X, the subbundle of ^'^TqX con- 
sisting of covectors that annihilate {^'^T{C; X))-^, is denoted by ^'^T*{C; X); we say 
that it is the relative scattering cotangent bundle of C in X. 
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This bundle of course depends on g. In the case of Euclidean scattering, ^'^T*{C;X) 
is naturally isomorphic to ^'^T^,Xa and is spanned by d{wa)j, j = 1, . . . , to. 

We now choose local coordinates (x, y, z) near p G C such that C is defined by 

X — 0, y ~ 0, and such that xdy- give an orthonormal basis of {'^'^T{C; X))-^. Note 

that a basis of '^'^T{C;X) is given by x^dx and xdz , while a basis of ^'^T*{C;X) 

is given by x^^ dx, x^^dzj. A covector in ^'^T*X can be written in these local 

coordinates as 

dx dy dz 

(4.13) T— +M- — + t^-— ■ 

x^ a; a; 

We will write this as 

, . , , dx dya dza 

X^ X X 

to emphasize the element C — Ca oi C around which the local coordinates are 
centered. Thus, local coordinates on ^'^Tq-^X are given by {y,z,T, ^i,u), while on 
^'^T*{C\ X) by (z, r, v) ~ {za, Ta^i^a)- Note also that at C the metric function of h 
is of the form |^p + h{z, v) with \ii\ denoting the Euclidean length of /x and h is the 
metric function of the restriction of h to TC; the metric function of g (also denoted 
by g) is thus 

(4.15) g = T'' + h+\ii\'' 

there. 

Now if C = Ca^ C'l, Cz C with Cq C C^, we can further adjust our coordinates 
so that C\ is defined by x = 0, y' = 0, for some splitting y = (y',y"). With the 
corresponding splitting of the dual variable, /i — {fi',fi"), we obtain a well-defined 
projection 

(4.16) TTba ■■ ^^r4(Cb;X) ^ ''^T*(Ca;X), 

(4.17) Trba{0,z,T,^i",u)^{z,T,iy). 

In the Euclidean setting this is just the obvious projection 

(4.18) TTba ■ ^'^Tg-^^Xb -^ ^'^Tg-^^Xa 

under the inclusion Xa C Xi,. We write tt for the collection of these maps. 

Before we define the indicial operators, we need to analyze the structure of the 
lift of Ca to [X;C]. For Ca e C let 

(4.19) Ca = {CbeC: CbCCa}, 

(4.20) r^{CbeC: CaCCb}. 

We carry out the blow-up [X] C] by first blowing up Ca- Since all elements of Ca are 
p-submanifolds of Ca, the lift f3[X-, Ca]*Ca of Ca to [X; Ca] is naturally diffeomorphic 
to 

(4.21) Ca^[Ca:Ca]. 

Thus, over C^, the regular part of Ca, Ca can be identified with Ca- The front face 
of the new blow-up, i.e. of the blow up of (3[X; Ca]*Ca in [X; Ca] is thus a hemisphere 
(i.e. ball) bundle over Ca, namely S~^NCa- We write the bundle projection, which 
is just the restriction of the new blow-down map to the front face, S^ NCa as 

(4.22) Pa : S+NCa ^ Ca- 
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In the Euclidean setting, these fibers can be naturally identified with X'^ via the 
projection tt" (extended as in Lemma ^.E| ). Every remaining blow up in [X;C] 
concerns submanifolds that are either disjoint from this new front face or are the 
lift of elements of C. The former do not affect the structure near the new front 
face, S~^NCa = P[X]Ca;Ca]*Ca, while the latter, which are given by the lifts 
of elements of C^, correspond to blow ups that can be performed in the fibers 
of S^NCa- Note that the lift of Cf, G C, meets the new front face only at its 
boundary since all Cb are subsets of dX. In particular, the lift Pg^Ca of Ca to 
[X;C] fibers over Ca and the fibers are diffeomorphic to a hemisphere (i.e. ball) 
with certain boundary submanifolds blown up. More specifically, the intersection of 
l3[X;Ca; Ca]*Cb, Cb e C^, with the front face S+NCa is the image of TP[X;Ca]*Cb 
under the quotients; P'sc^a is obtained by blowing these up in S'^NCa- Hence, the 
fiber of (3^^Ca over p € Ca is given by [S+NgCa; TqC] where q = /3[X; Ca] (p) G Ca- 
in particular, in the Euclidean setting, the fibers of Pg^Ca over Ca can be naturally 
identified with [X'';C°] via tt'^. Thus, we have the following commutative diagrams: 



(4.23) 



f^Sc^a 



^Ca 



/3sc 




Ca 



Psc^a 



Ca 



S+NCa 




with (3a being the fibration to the base Ca- 

We now define ^<=T*{Ca;X) denote the pull-back of '"'T*{Ca;X) by the blow- 
down map /3[Ca;Ca]- 



(4.24) ''T*iCa;X) 
If Ca C Cb then TTba lifts to a map 

(4.25) nba 



(3[Ca;Car'T*iCa;X)- 



^/5[Cb;Cb]*C„(^b'^) 



^T*{Ca;X) 



We recall from [^9[ Section 4] that the interior of the fibers S^NpCa = Pa ^(p) 
of Pa ■ S'^NCa -^ Ca, P & Ca, posscss a natural transitive free afRne action 
by the quotient bundle {P[X]Ca]*/''TX)/'"=Tp{Ca; X). Thus, the tangent space 
of S^NpCa at every point q e \ni{S^NpCa) can be naturally identified with 
(/3[X; Ca];''=rX)/"'=Tp(C'a; X), hence with the tangent space at other q' e \nt{S+NpCa). 

For each operator A G \I'^' (X, C), the Ca-indicial operator of A, denoted by 
Aa.i, will be a collection of operators, one for each C, G ^'^T*{Ca',X), acting on 
functions on the fiber [3a^{p) of /3a- So suppose that u e C°°{fia^{p)); we need to 
define Aa{C,)u. For this purpose choose / G C°°(X;R) such that d{f/x), evaluated 
at f3[Ca;Ca]{p), is equal to C- Then let A = e-'^^^x' 

= u. Then 



Ae'f/'^ G *™'°(X,C), and 



choose u' eC°°{[X;C]) such that u'l^^i/^^^ 

(4.26) ia,i(C)« = (i^')l^-(p), 

which is independent of all the choices we made. This can be shown by an argument 
which is analogous to the proof of the preceeding lemmas, but it will also follow 
from the explicit calculation we make below leading to ( 4.4C| ). li I ^ 0, then Aa^i 
would a priori depend on the choice of x up to O(x^) terms, but the choice of the 
scattering metric g fixes x up to such terms. We often simplify (and thereby abuse) 
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the notation and drop the index I, i.e. we write Aa = Aa.i, when the value of I 
is understood. Before discussing the Ca-indicial operators oi A E ^^' i^^C) in 
detail, we discuss how we can combine them into a single object. 

In the case of Euclidean many-body scattering, Ca = dXa and Aa^i is a function 
on (31^'^Tq Xa with values in operators on S{X'^); here 

(4.27) (3a=(3[Ca;Ca]:Ca=[Ca;Ca]'^Ca 

is the blow-down map. Note that (3a is simply the restriction of /3[Xa;Ca] to the 
lift Ca^P[Xa;Ca]*Ca. In fact, 

(4.28) Aa,i e C°°(/3r^^a*^^X„vI/™/(X^C'')) 

as we show shortly. Note that if Z is a (not necessarily compact) manifold with 
corners and {X,C) is a many-body space (in ( 4.28 ) we take Z — P*a^^T*^ Xa and 



{X°-,C°') for the many-body space), it makes perfectly good sense to talk about 
C°°(Z, ^P^' (X,C)), i.e. about smooth functions on Z with values in ^^' (X,C). 
The topology on $^' {X,C) is the standard one, namely that of conormal distri- 
butions on ^Ij,, conormal to Age, vanishing to infinite order at (3* bf U/3* If U/3* rf, 
/3 : Xg^ — !■ X^ the blow-down map. This is equivalent to the topology arising 



by localizing operators A e *gj.' {X,C) as in (3.17), and using the topology of 



the symbol spaces on the local pieces, i.e., with the notation of ( |3.11 ) and (3.5), 



of p^p'qC°^{[^\]C] X ^X) and ^-"'■'([§!;:;C] x ^\), in the polyhomogeneous and 
non-polyhomogeneous setting respectively (and that oiC°°{X x X-j^'^^r) for the 
remainder term). 

We need to generalize this example to accommodate the geometric setting. It 
should be kept in mind throughout following discussion that Z is simply a 'pa- 
rameter space'. So suppose first that </) : i? ^ Z is a fibration of manifolds with 
corners with fiber X, a manifold with boundary, Ce a cleanly intersecting family 
of p-submanifolds of E which is fibered over Z with fiber C, a cleanly intersect- 
ing family of p-submanifolds of dX that gives rise to a many-body space (X,C). 
That is, we suppose that there is an open cover {Uj : j G J} oi Z such that 
{(J)~^{Uj),Ce n (j)~^[Uj)) is diffeomorphic to Uj x {X,C); we denote the diffeomor- 
phism by tpj. Let d^E denote the fiber-boundary of E, i.e. locally it is given by 
Uj X dX (under the identification -0^). The algebra ^^'T°°(i<;,C£;) is then defined 
as the algebra of operators A acting on, say, functions u G C°°{E) which vanish 
to infinite order at d^E, with the following local characterization. For each Uj 
there is an operator A'^ £ C°°(C/j; *^'"°°(X,C)) such that for u e C°^{E) with 
suppu C (j)~^{Uj) and vanishing to infinite order at dcf,E, Au — ^*A'Atp~ )*u. 

This local description does not depend on any choices. Indeed, the local definition 
is equivalent to saying that the distribution kernel Ka of A on the fiber-product 
E Xz E (with values in scattering densities on the fiber X from the right factor, 
to be precise) is conormal on the appropriate blow-up E^^ ^ oi E Xz E. Here Ka 
gives rise to the operator A by fiber-integration 

(4.29) Au{w,z)^ KA{w,w',z)u{w',z)\dw'\, 



where z gives coordinates on Z, w and w' are variables in the left and right factor of 
the fiber X respectively, and we wrote Ka = Ka{w,w' ,z) \dw'\. Indeed, following 
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the discussion at the beginning of the previous section, we take E^ ^ to be the blow- 
up of d^E Xz dff,E in E Xz E, Ab,0 the Uft of the fiber-diagonal, S^Ab,^ its fiber- 
boundary which we identify with S^-E, C'^ the image of C^; under this identification, 
and -Egj. ^ the blow-up [E'^ z't^'e\- Then the definition of ^^'^°°(£^,Cb) is given 



by modifying (3.3) the natural way. Since all blowups can be done in the fibers 
over Z (i.e. Z can be regarded as a parameter), this description indeed agrees with 
local definition given above. 

This intrinsic definition of \1'^'T°°(£',Cb) given in the previous paragraph au- 
tomatically extends even to the setting where the fibration (j) is transversal to the 
collection C^, each fiber of being diffeomorphic to X. Note that in general there 
are no diffeomorphisms tpj even locally such that image of Ce takes a product form 
as above, though such diffeomorphisms exist, for example, if C^ is locally lineariz- 
able. In particular, we can take Z = ^'^T*{Ca;X), E to be the pull-back of Z to 
S^NCa by Pa, (f> : E ^ Z the map p^ induced by the pull-back, 

(4.30) E = p:^-T*{Ca;X), pi : E ^ ^-T*{Ca;X). 

Thus, E is a. vector bundle over S'^NCa with projection tt. Finally, we let Ce consist 
of the inverse images under tt of the lifts of Ct G C" to [X;Ca; Ca] intersected with 
the new front face, S'^NCa', in fact, we also add d^E to Ce to play the role of Co 
in C. We are then in the setting discussed above, so we have defined 

(4.31) 

*^jr(Pr7^*(^a;X),Ca), Ca = n-\S+NCan p[X;Ca;Ca]*C^)U {d4.E}. 

Recall that for CbEC, 

(4.32) S+NCanp[X;Ca;Ca]*Cb^Tp[X;Ca]*Cb, 

the right hand side understood as the image of the tangent space under the quotient 
map. We are now ready to prove the following proposition. 

Proposition 4.5. Suppose that A e ^g^(X,C). Then the indicial operators of A 
satisfy 

(4.33) Aa, e ^f^^^,yTT*{Ca-.X),Ca). 

Proof. We prove this statement by finding Aa (C) explicitly in terms of local coor- 
dinates. To simplify the notation we assume that A g ^I^g^, {X,C). We identify X 
with §" locally so that Ca is given by a; = 0, y = 0. In the interior of P^^Ca we 
can use the same coordinates as at the front face of [X]Ca], i.e. the ones given 



in (2/7)- (2^). So suppose that u' is supported in the region of validity of these 



coordinates. Then 

(4.34) Au'{w) = I K{w, w')w'(w') dw' = / d{w, W)u'{w - W) dW 

with a as in (p]q). Here the integral is understood as a distributional pairing in 
general, but it actually converges if r < —n. We now consider the coordinates 
( p?7| ) on the both factors, i.e. we take {x',Y',z') corresponding to w' ^ w — W, 
and (x, y, z) corresponding to w. Expressing (a;', Y' , z') in terms of (x, y, z) and W 
(using w' — w — W) gives 



(4.35) x' ^ x{i - x{WaU)-\ 4 - ? ,:;r^ n = Yj - {w^),. 



x{Wa)j 
l-x{Wa) ' '' '' 
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where we wrote W = {Wa, W'^) a nd {Wa)j, {W°')j denote the components of Wa 
and W^ respectively. Thus, (4.34) yields 



(4.36) 

Au'{x,Y,z)= f a{x,Y,z,W)u' 

Evaluating at a; = gives 



-,¥-]¥" 



x{Wa) 



l~x{Wa)m' 'l-x(Wa) 

Au'{0, Y,z)^ I 5(0, Y, z, W)u\0, Y - W, z) dW 



dW. 



a(0, Y, z, W) dWa w'(0, Y ~ W, z) dW . 



(4.37) 



Since d is the inverse Fourier transform in the ^ variable of the symbol a whose 
left quantization is A, and since the Wa integral above can be understood as the 
Fourier transform in Wa evaluated at the origin, we deduce that 

(4.38) 



ylw'(0,y,z) = (27r) 



— (n— 7n) / iW 



a(0, y, z, 0, r)w'(0, Y - VK", z) d^ dW". 



Thus, the indicial operator Aa{(p, 0)) where {p, 0) G ^"^T* {Ca\ X) is the zero covector 
above p = (0, 0, z) e Ca is given by 

(4.39) 

Aa{{j), Q))u{Y) = (27r)-("-'") / e*^°-«°a(0, Y, z, 0, i'')u{Y - W) d^ dW", 



i.e. by the left quantization in (Y, C) = (W", C") of a(0, F, z, 0, C)- Similar resuhs 
hold for Aa{C) in general, namely 

(4.40) 



Aa{z,^aHY) = i27ry 



[71— m) I ^iW°'-^'^ 



a(0, Y, z, ^„, CHY - W^"^) dr dVF^ 



Though the local coordinates are only valid in the interior of Pt^J^a, hence not at 
P*dCa, the continuity of Au up to P*dCa shows that ( 4.40| ) also holds with p € Ca- 
The explicit expression, ( 4.40| ) shows, in particular, that Aa(C)w is indeed inde- 
pendent of the extension u' of u that we chose, and also of the choice of / with 



d{f/x) prescribed at Psc{p)- Moreover, also from ( 4.40 ), for each (^ e ^'^T*{Ca;X), 

pe Ca, 

(4.41) Aa{0&%''(.Pa\p),T,n; 



here we wrote TpC"- for TpP[X;Ca]*C°- for simphcity. In fact, (4.40) shows the more 
precise statement which encodes the smooth dependence of Aa{C) on (, namely 
that 



(4.42) 



Aa,ie^''/l^APrT*iCa;X),Ca). 



In the Euclidean setting the many-body space {pa ^(p), TpC^) can be identified with 
(X^jC"), and we can write 



(4.43) 



^(C)e*^f(x^r) 
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and correspondingly 

(4.44) Aa,i e C^iPT^ryXa, vI;^'0(X^C'^)) 

as we have claimed. D 

If A e ^gj,(X, C), then the vanishing of AafiiC) for every a and every C, G 
scT*{Ca;X) impHes, by our explicit formula, that a S C°^{[X]C] x S") vanishes at 
{d[X\C]) X W%, so A e *g'^^(X,C). Thus, the vanishing of crsc,r(^) and aU indicial 
operators together, for A € Vfg^ {X,C), say, implies that A e ^gJT ' {X,C). 

An advantage of the oscillatory testing definition of the indicial operators is that 
it makes their multiplicative property clear. 

Proposition 4.6. If Ae ^'^^\X,C), B G ^'g/ {X,C) then 

(4.45) ABa.i+v{Ou = Aaj{C)Ba,i'{Ou- 

The indicial operators are related via the projections iiba ■ Thus, if C G ^'^T* {Ca\X), 
then the indicial operators of Aa,i{C,) are Ah^i{Q where Ga <Z Cb, Ca ^ Ch, and 
C e ^^T"b\c c I'Ca^^''''-^^ ^^ ^^''^ *^^* ^baiO — C- This follows easily from the 
explicit coordinate form of the indicial operators. In particular, if A e ^^' {X,C) 
with TO < 0, and if Ab.iiC) vanishes for all such b and (, then Aa^i{C) is compact. 
We thus have the following proposition. 

Proposition 4.7. //A e ^^^^(X, C) is such that a sc,r{A) never vanishes and Aa{C) 
is invertible in ^^'^ {Pa^ip)^ TpC"") (i.e. in ^P^'^ (X", C) in the Euclidean setting) for 
every a and for every ^ G scjr* j-^^^. j^-^^ i/ien t/iere exists a parametrix P for A such 
that PA - Id, AP - Id G *s,°°'°° (X, C) . 

Proof. Choose Pq G ^sT' (-'^jC) with principal symbol crsc,r(^)~^ and indicial op- 
erators AafliC)^^- Note that these match up under the projections nba since those of 
A match up. Thus, they indeed specify the restriction of some po & p'^C°°{[X; C] x 
§" ) to the boundary of [X; C] x S" . Quantizing pq gives an operator Pq with the re- 
quired indicial operators and principal symbol. Hence, E — Id— PqA e ^g^ (^tC) 
has vanishing principal symbol and indicial operators, so it is in ^Pg^ ' {X,C). Sum- 



ming the Neumann series J^iLi -^"' asymptotically to some F G ^g^. ' (X, C) and 
letting P = (Id -|-P)Po gives the required left parametrix. A right parametrix can 
be constructed similarly, and then the usual argument shows that they can be taken 
to be the same. D 

For A g VP^' {X,C) self-adjoint, to > 0, with asc,miA) never vanishing, we 
automatically have that {A - A)"^ e *g™'°(X,C) for A G C \ E. Moreover, the 
blow-up of {A — A)^^ in ^g^™' {X,C) can be analyzed uniformly as A approaches 
the real axis, see e.g. P, p9]. Therefore, the functional calculus for self-adjoint 
operators A and the Cauchy integral representation of </'(A) via almost analytic 
extensions, as in the work of Helffer and Sjostrand iQ, Derezihski and Gerard 0, 
see also |9|, gives immediately 

Proposition 4.8. Suppose that A e ^™J {X,C) self-adjoint, m > 0, and asc,miA) 
never vanishes. Suppose also that (j) g C^(M). Then (j){A) G ^5^' (^jC) and its 
indicial operators are 0(Aq(C)). If instead we assume G S~L{M.) then 4'{A) G 

*^r'°(x,c). 
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Ifm = 0, that isA e *gf(X,C), then 0(A) e *g;°(X,C) without any assumption 
on the invertibihty of crsc.o(^)- We thus have: 

Proposition 4.9. Suppose that A G *gf(X,C) self-adjoint. If (j) e C°°(R) then 

Proof. Since A is bounded, we can replace ^ by a function tp G C^(R) such that 
= f/; on the spectrum of A. Now agcfliA — A) = (Tsc.o(^) — A is invertible for 
A e C \ R, so (A - A)-i e *sc°(^' ^) fo'' ^ ^ K- Again, (A - A)-i can be analyzed 
uniformly up to the real axis, and then the Cauchy integral representation of ^p{A) 
now proves the proposition. D 

5. The wave front set 



front set WFg(,(A) of A S ^^' {X,C), at infinity will be defined as subsets of the 



The Sc-wave front set WFsc(u) of a distribution u, and the Sc-operator wave 
nt set WFg^(A) of A e ^™''(X,C), at infinity will 
npressed scattering cotangent bundle 

(5.1) ^^T*X = U/^T5,(Ca;X); 



we have defined ^'^T*{Ca;X) in Definition 4.4. This is very similar to the image of 
the cotangent bundle in the compressed cotangent bundle (the b-cotangent bundle) 
that Melrose and Sjostrand used to describe the propagation of singularities for the 
wave equation in domains with smooth boundaries [E2| and also to the correspond- 
ing phase space for domains with corners O, T^Jl, which was the setting for Lebeau's 
analysis of the singularities of solutions to the wave equation on f2. We topologize 
^''f*X using the projection tt : """T^x^ "* '"'f*X. Thus, WFsc(m) will contain less 
detailed information than WFsc(m), the latter being a subset of ^"^Tq-^X. However, 
neither of these wave front sets can be used to describe the other; they are simply 
different. The fact that WFsc('u) contains fewer details about u simply corresponds 
to the singular behavior of elements of ^Pg"' {X,C), compared to those of 4"™' (X). 

The definition of WFsc(u) and WFg(.(A) will be local in X. Thus, we can 
always work on S" instead. Just like when we defined \E'^' {X,C), we will be able 
to proceed either by giving an explicit description in S" via the Fourier transform, 
or by giving invariant definitions. 

We start with the operator wave front sets. The invariant definition proceeds by 
oscillatory testing. 

Definition 5.1. Suppose that A e *s"''(X,C) and C G "'=T;(Ca;X), p e C^. We 
say that C ^ WFgj,(yl) if and only if there exist a neighborhood f/ of ( in ^'^T*X 
and a neighborhood V oi p in X such that Au G C°°{X) for every oscillatory 
function u = e^^/^w, v G C°°{[X;C]) with 7r(graph(d(//a;))) n """f^^gx^ C U and 
su-ppv d (3s^\V). 

This definition implies immediately that WFgj.(A) is closed in ^'^T*X, 

(5.2) WF'g,(A + i?)cWF'g,(A)UWF^,(i?), A,B e ^sc{X,C), 

(5.3) WF'g,(Ai3)cWF'g,(A)nWF^,(B), A,B e ^sc{X,C). 

We can also formulate the definition explicitly. We thus locally identify X with 
S^ and consider A £ *™/(S!^, C). We also identify ^'^T*§1 with §!^ xR". So suppose 
that A is the left quantization of a symbol a £ p^"'p^gC°°{[Sl;C] x §1). Then 
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C i WFg^(A), C e ''''T*{Ca\ X), p e C^, if and only if there exists a neighborhood 
[/ of C in ^'^t*§l such that a vanishes at U' C (a[§!;:;C]) x M" to infinite order 
where U' is the inverse image of U under the composite map 

(5.4) (a[S!^;C]) X M" I:t^ (5§n) x M« = ='=Tg*„_i§'^ — ^^ «'=T*S!^. 

Note that this condition impUes automatically that a vanishes to infinite order on 
the closure of U' in (9[§";C]) x §", in particular, it is rapidly decreasing in some 
directions as |^| — > oo. It follows immediately from the usual formulae relating 
quantizations and the effect of diffeomorphisms that this definition is independent 
of such choices. For example, we could have equally well written A as the right 
quantization of a symbol with similar properties. 

In general, the structure of tt can be somewhat complicated in explicit coordi- 
nates. However, in the actual Euclidean setting, i.e. where C arises from a family 
of linear subspaces X, or indeed, if we merely assume that {X, C) is locally lineariz- 
able, it is particularly easy to give a sufficient condition for ( ^ WFg^. {A) . Namely, 
if there is a neighborhood V oi C = (0, z°, ^°) G "'=T;(Ca; S!^), p G C^, in dS^ x R™ 
such that a vanishes to infinite order at every point {q',£,) G (9[S";C]) x R" with 
{Psciq'),^a) € V, then ^ ^ WF'g^{A). Note that as p G C^, we can always assume, 
by reducing the size of V if necessary, that (g, C°) G V implies q & C'^ for some b 
with Ca C Cb- We can see that this condition is sufficient for C ^ WFg(,(A) since 
for nearby q G §"~^, assuming as we may that q E C'^, Ca C Cb, the restriction of 
TT to «=T;S!^ takes the form (g, 6, C^) ^ {q, ^b) and ^ splits as (C^, C") with C^ = ^a- 
Thus, the condition of the previous paragraph holds if we take 

(5.5) U - Ub{{q, 6) : g e C'b, 3^a, Cb s.t. (g, ^a) e V and 6 = (^;, Cb)}- 

The definition of the wave front set of a distribution u G C^°°{X) at dX is more 
complicated. To determine whether ( G '^'^Tp (Cq;X), p G C^ is in WFsc(u), we 
would like to cut off u to be supported near p, i.e. consider t/ju, ip G C°°{X), ip = 1 
near p, identify a neighborhood of p with an open set in S" near 9§" , and consider 
smoothness of the Fourier transform of u, Tipu. Indeed, in the two-body setting, 
hence in the many-body setting if we consider C, G ^'^T^, (Co; X) — ^^T^,X, written 
as a covector ^ • dw over p £ C'q, we have 

(5.6) C ^ WFsc(u) iff 3ip as above, s.t. Tipu is smooth near ^. 

In the general many-body setting, ( G ^'^T*{Ca;X), p G C^, C takes the form 
^a • dwa, and correspondingly we would like to say that J^Tpu is Schwartz in a 
region including the subspace S consisting of all points of the form (^a,C) where 
f is arbitrary. Here Schwartz takes the place of smooth since the region is not 
compact in R". However, as shown by the example of ordinary wave front set, 
we cannot expect that this wave front set behaves reasonably unless the region 
U is conic near infinity, i.e. unless it is a neighborhood of the closure of 5* in the 
radial compactification S" of R". This however introduces the complication that 
all parallel translates of S intersect U, and we are exactly interested in separating 
from each other the singularities on the various translates of S. This problem is not 
too serious, especially for generalized eigenfunctions of many-body Hamiltonians H , 
but it introduces additional terms into the following definition which is modelled 



on that of the fibred cusp wave front set by Mazzeo and Melrose |18 
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Definition 5.2. We say that 

(5.7) 

C i WFsc(u) n ^^T*c,^((Ja\X) iff 3A e *°f (X,C), ia,o(C) invertible in *°f (X^C") 

Here we used the Euclidean notation *gf (X°,C") instead of *sc°(Pa ^(p)> ^pC") for 
the sake of simphcity. Similarly, the filtered version of the Sc-wave front set is given 

by 

(5.8) 



C ^ WF^/(u) n ^^ra, (Ca;X) iff 3A e *^;y(X,C), ^a,o(C) invertible in vI/^-'^(X",C"), 

35, e *g-°°'°(X,C), C^WF^,(i?,), 
3u, e C-°-\X), J = 1, . . . , s, / g iI™''(X), 

s 

Thus, if p G C^, then the part of WFgc over p lives in '"'T*{Ca;X). If we 
define the scattering wave front set, WFgc (u), in terms of operators instead of the 
description of WFgc (w) given in ( |5.6| ) then the extra terms BjUj can be dropped. 



In fact, (5^) is equivalent to requiring that Au e C°°(§" ) where A — T ^(pTt/j G 
*rc°°'°(§+), ip as above, and (j) S C;?°(K) is identically 1 near ^. The additional 
terms BjUj for WFsc(u) thus arise because the invertibility of Aa{C) implies that 
o'Sc,o(^a(C)) cannot vanish which in turn means that crsc,o(^a(C')) i^ non-zero for 
every C' e ''''T*{Ca;X) since crsc,o(^a(C)) = CTSc,o(^a(C'))- This simply corresponds 
to the conic cutoff requirement discussed before the definition. 

With the topology we put on ^'^T*X, WFsc(w) is closed due to the relationship 
between the indicial operators mentioned above. Namely, the invertibility of Aa^o(C) 
implies that of Ab^o(C) with TtbaiC) = C: hence of Ab.o(C') for nearby (' . As the 
complement of WFgj.(_Bj) is open, this implies that the complement of WFsc(u) is 
also open. In addition, 

(5.9) WFsc(mi + U2) C WFsc(ui) U WFsc(m2) 

and the corresponding result also holds for the filtered wave front set. Moreover, 

(5.10) A e ^l'^'iX,C), u G C-°^{X) ^ WFse(Aw) c WF's,(A) n WFse(w), 
and similarly 

(5.11) 

A e *'s"/(x,c), u e C-°°{X) ^ wf^J-™^'+''(Am) c WF's,(A) n wf^J'''(u). 

We refer to |18[ for detailed arguments; we only need simple modifications of their 



proofs. We also refer to the remarks after Proposition 5.3 for connecting this wave 



front set to the one discussed in |39| in three-body scattering. 
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This wave front set, WFgc, gives a complete microlocal description of distribu- 
tions at dX. To state it generally, we would need to define the extension of the 
standard wave front set of u to give a subset of ^'^S*[X\C], but for us the following 



extension of (5.10) suffices: 
(5.12) 

A £ *g^°°''(X,C), u e C'°°{X), WF^,(A) n WFsc(m) = ^ A^i e C°^{X). 

We remark that in [gl|, WFsc(w) (or rather its part over dX) is defined as a subset 
of ^"^Tg-^X, the radial compactification of ^'^Tg-^X in the fibers. The part at fiber- 
infinity, i.e. at the boundary arising from the radial compactification of the fibers, 
extends the usual wave front set from the interior. However, for us this extension is 
not important; the operator wave front set of nearly all operators we are interested 
in is contained in a compact region of ^'^T*X. 

The description of the wave front set becomes simpler for generalized eigenfunc- 
tions of many-body Hamiltonians H . Namely, we have the following result. 

Proposition 5.3. Suppose that u G C"°°(X), H € *™j"(X,C), m > is self- 
adjoint and <Jsc,m{H) never vanishes. Let A € M, and define W C '^'^T* X by 

(5.13) 

C ^ PFn ^=T5,(a;X) tfJ3^p e Cr(K), V^A) = 1, 

3A e ^-s:^-\x,C), AaiC) = m)a^ Au e C^{X). 
Then 
(5.14) WF5c(w) C WFsc((ff - A)w) U W. 

The same conclusion holds with WF^c replaced by WF^'^' and Au G C°°{X) by 
AueHl^^iX). 



Proof Suppose that C ^ WFsc((-ff - X)u) and ( ^W. With -0 as above, let -ipit) = 

Jp-4(R) as V(A) = 1. Then ^(H) G *"" 



(1 - ijit))/it - A), soij e S-^^m as V(A) = 1. Then ^{H) € *-'"'°(X,C) and 



ld = ^lj{H){H-X) + 'ilj{H). WithAasabove, let A' = A+(Id-^(i7)) e *gf(X,C). 
Theni;.(C) = Id and 

(5.15) A'u^Au + 'tP{H){H~X)u. 
But Au e C°°{X) by assumption, so by ( ^.10| ) 

(5.16) WFsM'u) = WFsc{i'{H){H - X)u) C WFsc((i7- - X)u). 



Hence, there exist A" (in place of A), Bj, etc., as in Definition 5.2, A"A'u — 
f + J2 Bj'^jj ^^d the indicial operator of A" A' at ( is just the composite of those 
of A" and A' , hence invertible, showing that C ^ WFsc(u). D 

Remark 5.4. Our definition of WFsc(u), which is in particular valid if {X,C) is 
a three-body space, is different from the wave front set WF3sc(u) used in 139] in 
the three-body setting. Indeed, in the definition of WF3sc(u), the terms BjUj 
appearing in Definition |5.2| were not allowed. Consequently, ( 5.10| ), more precisely 



WF^sciAu) C WFssciu), and its filtered analogue did not hold in general. However, 



for the positive commutator proofs of both |39 and the present paper, one only 
needs ( p3.12D , which was proved for WFs^c under the weak additional condition 
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that WF3g^(A) is compact; this condition holds for aU operators appearing in such 
estimates in both papers. 

Note that WFscIm) C WF^sciu) directly from the definition. Moreover, if 
(Id-F)u e C°°{X) for some P € *sr'°(X,C) (e.g. P == ?A(iJ) in the setting 
of the proposition) then WFsc(u) — WF^sciu). In fact, suppose that ( ^ WFsc(u), 



so Au — '^BjUj + / as in Definition 5.2. Since A is invertiblc near (^, we can 



arrange (by inverting A nearby, i.e. by constructing a 'microlocal parametrix') that 
A'u = ^ B'jUj + /' with A' and B', as in the definition, but in addition with 
C ^ WFg^(Id— A') (cf. Q). Using the methods of Section^, given any neigh- 
borhood U of C, it is easy to construct an operator C G ^3^°°' {X^C) such that 
WFg(,(C) C U and C, ^ WFgj,(P - C) (hence the same holds for a neighborhood 
of C). Since the indicial operator of Q = C + (Id— P) at C, is the identity, and 
since (Id— P)u G C°°{X), we only need to prove that Cu £ C°°{X) to conclude 
that C ^ WFssciu). But Cu = i^CB'^^j + Cf + C(Id-A')u, so if U is chosen 
sufficiently small, then CB'^, C(Id-A') e *g^°°'°°(X,C), so Cu e C°°(X) indeed. 

6. The Hamiltonian and generalized broken bicharacteristics 

We next analyze the operator H — X where H = A + V and A is the Laplacian 
of a scattering metric 

, , dx2 h' 

Recall that h' is a smooth symmetric 2-cotensor on X whose restriction to dX (i.e. 
its pull-back), h, is positive definite. We assume that 

(6.2) V e C°°([X;C];R) vanishes at P^^Co, 
i.e. V vanishes in the free region. This implies that 

(6.3) HeBiSl^iX,C). 

Such a situation arises, for example, in actual Euclidean scattering if the potentials 
Va (in the notation of the introduction) are classical symbols of order —1 on X"'. 
Hence, we make the following definition. 

Definition 6.1. A many-body Hamiltonian is an operator H — A + V where A is 
the Laplacian of a scattering metric g, and V satisfies (p^). 

As indicated in the Introduction, from this point on we also make the assumption 

(6.4) (^jC) is locally linearizable; 

this will simplify the analysis. We recall that this is equivalent to the local existence 
of Riemannian metrics on dX, possibly different from h, with respect to which all 
elements of C are totally geodesic. 

Since CTsc,2(A) never vanishes, the same holds for crsc,2{H) which is the pull-back 
of the former. A simple calculation, see [p3. Sections 4 and 11] for more details, 
shows that the indicial operators of H are given by 

(6.5) ff,,o(0=^a,o((p,0))+r2+/i(z,z.), ^^{z,T,iy)e''T*{Ca;X), 

(6.6) H,^o{p,0)^^Y + V{p,Y) 



MANY-BODY SCATTERING 35 

where Y are 'Euclidean coordinates' on the interior of p'^^{p), i.e. that of /3~^(p), 
and Ay is the Euclidean Laplacian. 

More precisely, we have seen in Section| that {P[X;Ca]*p''''TX)/''''Tp{Ca; X) nat- 
urally acts transitively and freely on the interior oi p~^{p) = S^NpCa, so it makes 
sense to talk about translation invariant vector fields and differential operators on 
the interior of S^NpCa- Indeed, the restriction to S~^NCa of the lift of elements of 
DiSsc{X) (under /3sc) are such. We can see this since Vsc(-^) is given by sections 
of '^'^TX; the restriction of the lift of P G Vsc(-'^) is then given by the identification 
of {P[X;Ca]*p^''TX)/'"'Tp{Ca]X) with the tangent space at each point of the fiber 
p~^(p). Using the metric g to identify the quotient bundle with the orthocomple- 
ment of ^^Tp{Ca\ X), S^ NpCa becomes an afHne space with a translation-invariant 
metric (i.e. 'Euclidean') with the metric induced by g] Ay is the Laplacian of this 
metric. 

Equations ( |6.5| )-( |6^ ) show that Hafl{p,0) is uniformly bounded below, so for 
any V- e C;?^ (R) the set 

(6.7) U,cl({ee^^r*(C„;X): ^(i?„(0) 7^ 0}) 

is compact. 

The bound states of the subsystems of H play an important role in Euclidean 
many-body scattering. The appropriate replacement in the general geometric set- 
ting is given via the indicial operators of H. Thus, in this paper the statement 'no 
subsystem of H has a bound state' means that 

(6.8) HaAO has no L^ eigenvalues for any a 7^ and ^ e ''^r*(C'a; X). 



Due to (3.f:)-(p^), this means simply that 

(6.9) ha{p) = Ha.oiip, 0)) has no L^ eigenvalues for any a 7^ and p & Ca- 
in Euclidean scattering ha{p) is just t he s ubsystem Hamiltonian ha (which is then 



independent of p), so in that setting (x8) indeed means that the (proper) subsys- 
tems of H have no bound states. 

If no subsystem of H has bound states it can be expected that A — A governs 
the propagation of singularities of distributions u with {H — X)u £ C°°{X), except 
that the flow will break at the places where V is singular (i.e. where locally V ^ 
C°°{X)), similarly to boundary and transmission problems for the wave equation 
jn], Chapter XXIV], (2|, 0. Now, the symbol of A - A at 9X (i.e. its sc-indicial 
operator) is g — X. Hence, its characteristic variety is 

(6.10) S = Ea-a - U e ^%V^ : .9(0 - A = 0}. 

The rcscalcd Hamilton vector field ^'^Hg = x^^Hg of g (or q — A), introduced in 

|2l),is 

(6.11) '■''Hg ^ 2T{xd^ +p-d^ + i^-d^)~ 2hdr +Hh + xW\ W e Vb('"r*X), 
so its restriction to dX , also denoted by ^'^Hg, is 

(6.12) "'^iJg == 2t{pl ■df. + i^-d^)- 2hdr + Hh- 

Here {y, z, r, /i, v) denote coordinates about some C = Ca as before, though notice 
that ji-d^ + v -di, is simply the radial vector field in T*dX, so the above expression 
is indeed invariant (as it must be). The bicharacteristics of A — A are just integral 
curves oi^'^Hg. 
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We divide the image S C ^'^T*X of E under n into a normal and a tangential 
part, 

(6.13) i] = S„(A)UEt(A), 

as follows. Let tt be the restriction of n to E. We let 
(6.14) 

E„(A) = Ua{C e *"'T'5/ (Cq; X) n i; : -n-^^iO consists of more than one point} 

and 

(6.15) Et(A) = Ua{C e "'^T^, {Ca;X) n E : 7r"i(0 consists of exactly one point}. 



In terms of our local coordinates around C^, in view of (4.15) and |/iaP > 0, this 
means that 

(6.16) E„(A)=Ua{(za,r„^,)G^^T5,(Ca;X): T^ + h{Za, l^a) < X} 
and 

(6.17) ^t{X) =Ua{{Za,Ta,iya) e'''T^,JCa;X) : T^ + h{Za, l^a) = X} ■ 

Notice that for ^ — {za, Ta, Va) G ^tiX) and the unique point ^ = (0, Zq, Ta, fia, ^a) G 
'^'^Tgj^X with 7r(^) = ^ we have ^a = 0. As the dy^ component of ^'^Hg is 2/^q • 9j,^ 
at 2/a = (i.e. at Cq), for such ^ and ^, ^'^Hg{S,) is tangent to '^'^T^, X. On the other 

hand, if C G E„(A), | G 7r^^(C), then ''''Hg{^) is normal to "'^T^/ X, hence the choice 
of our terminology. Notice also that on ^"^T^iX, tt is the identity map, so 

(6.18) Sn'^^T^.J^C Et(A). 
We also define the radial sets R±{X) as the sets 

(6.19) R± (A) = TT{{{y, z, T,fi,iy): r = ±VX, h{y, z, ^, ly) = 0}). 

Thus, i?+(A)Ui?_(A) is the image (under tt) of the set where ^"^Hg vanishes. Notice 
that 

(6.20) i?+(A)Ui?_(A) cEf(A). 

Following Lebeau, we define generalized broken bicharacteristics of A — A as 
follows. First, we say that a function / S C°°{^'^Tq-^X) is 7r-invariant if for ^,^' G 
^'^TqxX, 7r(f) = 7r(^') imphes /(f) = /(!')• ^ 7r-invariant function / naturally 
defines a function fn on ^^T*X by fniO = /(O where f G '''^Tgj^X is chosen so 
that 7r(f) = ^. 

Definition 6.2. Suppose that {X,C) is locally linearizable. A generalized broken 
bicharacteristic of A — A is a continuous map -f : I ^ ^'^T*X, where / C M is an 
interval, satisfying the following requirements: 

(i) If ^0 = 7(io) G Ei(A) then for all 7r-invariant functions / G C°°("'=Tg*^X), 

(6.21) |(/. o 7)(fo) - ^'Hgf i^o), Co - ^-\Co). 

(ii) If ^0 = 7(to) G E„(A) n """T^, iCa;X) then there exists e > such that 

(6.22) tel, 0<\t~t„\<e^^{t)^^'^Tl,,{Ca;X). 
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The success of this definition (so that it indeed describes what we wish to de- 
scribe) depends on a plentiful supply of 7r-invariant functions on ^'^Tg^X. Under 
our local linearizability hypothesis, (|6.4|), there are always many such functions. 



Indeed, by (|6.4|), if p G C" = C^, we can choose local coordinates {y,z) on dX in 
terms of which all the Cb satisfying p £ Cb are linear, i.e. they are are given by 
Aby = where Ab is a (constant) matrix, and Ca is given by y = 0. Let (r, /i, ly) de- 
note the sc-dual variables of {x, y, z) as above. Choosing such coordinates, y, z, r, v 
are 7r-invariant near ^"^T^X. In general, without the assumption (6.4), v would 



not be TT-invariant, and we would not be able to modify it to make it such, so the 
definition would be inadequate. 

We can also arrange that the metric function is of the form h = h(z, v) + 
hnn{z, /i) at C^ by a further change of coordinates Zj = Zj + J2jk Zjk{z)yk, y' = y, 
which preserves the linear structure of the 0%. In general we cannot arrange that 
hnn{z,n) = ImP everywhere along C'^ without destroying the product-Hnear struc- 
ture of the Cb- However, by a linear change in the y coordinates we can make 
sure that at a fixed p G C^, h = ft,(z, z^) -f |/ip. The continuity of a generalized 
broken bicharacteristic 7 means that if 7(^0) G ^"^T^'iC'jX), then for t near to, 
t H^ {y{l{t)), z{-f{t)), T{'j{t)), i/(7(t))) is continuous, but ^{■j{t)) may be discontinu- 
ous. In terms of Euclidean scattering this means that at Ca the external momentum 
is conserved, but not necessarily the internal one, while image(7) C S corresponds 
to the conservation of kinetic energy. The latter cannot be expected to hold if the 
subsystems of the Hamiltonian have bound states; the relevant broken bicharacter- 
istics in that case exhibit more complex behavior. Another example of a 7r-invariant 
function in this situation is y ■ fi; this will play a rather important role in the prop- 
agation estimates. In fact, ^'^Hg{y ■ /i)(^o) = 2|/ioP if Co G ^^T*X is of the form 
(0, zo, ro,/io, I'o), so if 7r(Co) G S„(A) and ^o G S then y • /x is a parameter along 
generalized broken bicharacteristics near ^o ~ see also the following proposition. 

A stronger characterization of generalized broken bicharacteristics at I]„(A) fol- 
lows as in Lebeau's paper. Notice that if 7 : / ^ S is continuous then the conclusion 
of the following proposition certainly implies (i) and (ii) ((ii) follows as yj — {ya)j 
are 7r-invariant), so the proposition indeed provides an alternative to our definition. 

Proposition 6.3. (Lebeau, |17| , Proposition 1]) Ifj is a generalized broken bichar- 
acteristic as above, to € /, ^o = 7(^0); then there exist unique ^+,^^ G I](A — A) sat- 
isfying 7r(C±) = Co o.'i^d, having the property that if f € C°°{'''^Tq-^X) is ■n-invariant 
then t I— > fjri'jit)) is differentiable both from the left and from the right at to and 

d_ 
di. 



(6.23) (-)(/-°7)l*o± = ^'ffJ(l±)- 



We refer to Lebeau's paper for the proof in the general setting, but in the Ap- 
pendix we give the proof under the assumption that the elements of C are totally 
geodesic. In fact, we prove slightly more by giving a Holder-type remainder esti- 
mate. We present the proof in the Appendix, but we emphasize that it is simply 
a minor modification of Lebeau's proof. We remark that the most delicate part 
of the conclusion (under the totally geodesic assumption) is the differentiability of 
the 'normal' coordinate functions yj along 7, i.e. that of yj o 7. Here we dropped 
the projection w from the notation (i.e. we did not write (yj)^ o 7) to simplify it; 
we will often do this in the future for the other 7r-invariant coordinate functions r, 
Zj, Vj. The proof proceeds by induction using the order on C. Thus, we have to 
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understand what happens near tg if 7(^0) = Co G S„(A) n ^'^T^, {Ca;X). The induc- 
tive hypothesis is that we have already proved the proposition for b with Ca Si Cb- 
Thus, by Definition |6.2| , part (ii), it is true for to replaced by i ^ to, assuming 
|i — io| < £• Hence, we need to analyze the behavior of the coordinate functions 
using the Hamilton equation, (6.23|) whic h is a little more delicate than the positive 



commutator construction in Proposition 10.3, but the two proofs are very closely 



related via the use of same function (j) to localize near (and along) the generalized 
broken bicharacteristics. A rather similar analo gy a rises in ou r tangential estimates 



in the totally geodesic setting; see Propositions [7T| and |10.4| respectively. 

We now describe some corollaries of this proposition. First, we remark that the 
role of the globally defined 7r-invariant function t is somewhat analogous to the role 
played by the time variable in the wave equation in Lebeau's paper. In particular, 
T gives a parameter along generalized broken bicharacteristics with the exception 
of some trivial ones (namely the constant ones in R+{X) U i?_(A)). To see this, we 
show the following corollary of the above proposition. 

Corollary 6.4. Suppose that 7 : / — > S is a generalized broken bicharacteristic. 
Then T — Tj^ o j : I —> R is a C°° function. In addition, T has one of the following 
forms. Either 

(i) T{t) = y/X for all tel, or 

(ii) T{t) ^ -Vt for all tel, or 

(iii) T'{t) < for all t and if I ^R then T{t) -^ tVA as t -^ ±cx). 

Proof. As A = T^ + /i in Sa-a, we have for all ^ G Sa-a that 

(6.24) """HgTii) = -2hiO - 2(t(|)2 - A). 

Thus, with T = Tt^ o J, the previous proposition implies that for any t £ I, T is 
differentiablc from both the left and the right at t, and both of these derivatives 
are equal to 2(T(i)^ — A). (We remark that this is proved directly in the Appendix 
as a first step to the proof of the proposition.) Thus, T is C^ and it satisfies the 
ODE dT/dt = 2(T2 - A). But, given say T(io) ^ tq, this ODE has a unique 
solution which is C°°. The last statement follows by writing down the solution 
of the ODE explicitly, which, if T{to) E (—A, A) for some to £ I, takes the form 
T{t) = — vAtanh(4\/A(t — c)), t £ I, for an appropriate constant c. D 

Since for ^ g S with r(C)^ = A we automatically have ^ E R+W U i?_(A), in 
(iii) we see that (if / = R) as t ^- ±00, ^{t) approaches i?ip(A). In addition, in the 
same case, as T' never vanishes, T E (— vA, vA) can be used to reparameterize 7 
(reversing its direction). 

We proceed to examine generalized broken bicharacteristics in more detail, start- 
ing with cases (i) and (ii). Namely, we prove that generalized broken bicharacter- 
istics through i?+(A) U -R-(A) are constant maps: 

Proposition 6.5. // 7 : / — > E is a generalized broken bicharacteristic, 7(^0) = 
Co E i?+(A) U i?_(A), then j{t) — Co for t E I. Hence, tt~^ o j is a bicharacteristic 

Proof. The previous corollary and the above remarks show that for all t E I, 7(t) E 
R+{X) U R-{X). Let |(i) = 7r-^{j{t)). Thus, ^''Hg vanishes at i{t) E n-'^{R+{X) U 
i?_(A)) for all t. Since the base variables y and z are 7r-invariant, we conclude that 
d{{yj)TrO^)/dt vanishes identically, hence y is constant, and similarly for z, proving 
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that j{t) — ^0 for all t. The last statement of the proposition follows since ^"^Hg 
vanishes at 7r-i(i?+(A)Ui?_(A)). D 



Now, we consider case (iii) of Corollary 6.4. Namely, we show that if we rescale 
and reparametrize 7 and project off its r component, we obtain a generalized broken 
geodesic (of h) in dX, broken at C. This is a notion completely analogous to that 
of our generalized broken bicharacteristics, and we proceed to define it. Again, we 
need to introduce a 'compressed' cotangent bundle. The metric h on dX naturally 
identifies the cotangent bundle T*C of C € C as a subset of T*dX. The compressed 
cotangent bundle of dX is then 

(6.25) f*dX = U„T5,a. 

It is topologized by the projection ttq : T*dX -^ T*dX . We also define the com- 
pressed cosphere bundle as the image of S*dX under Trg; here S*dX is the set of 
covcectors of unit length: 

(6.26) S*dX = TTa{S*dX). 

The restriction of ttq to S*dX is denoted by ttq. This plays a role analogous to 
that of S. We also define its tangential and normal parts: 

(6.27) 

S^dX = Ua{C € Tq, Ca n S*dX : ttq^{C,) consists of more than one point} 

and 

(6.28) SldX = UJC e T^'Ca n S*dX : ng\0 consists of exactly one point}. 

Generalized broken geodesies are then defined as follows. 

Definition 6.6. A generalized broken geodesic of /i is a continuous map 7a : / ^ 
S*dX , where / C M is an interval, satisfying the following requirements: 

(i) If Co = ld{to) G SfdX then for all 7rg-invariant functions / e C°°{T*dX), 

(6.29) ^(/., o 7a)(to) = i/i,J(Co), Co = TTg i(Co). 

(u) If Co = 79(^0) e S*^dX n T^, Ca then there exists e > such that 

(6.30) f e /, < |i - tol < e ^ 7a(t) i T^,^Ca. 

Remark 6.7. Sometimes, with an abuse of terminology, we also say that the pro- 
jection of a generalized broken geodesic to dX (via the projection S*dX — > dX 
inherited from TdX) is a generalized broken geodesic. Indeed, this was the termi- 
nology used in the introduction. 

The metric g gives rise to a product decomposition 

(6.31) """ToxX = K^ X T*dX. 

The compressed scattering cotangent bundle is thus also naturally a product: 

(6.32) scy*^ = M^ X t*dX. 

We sometimes write the product variables as ^ = (r, ^"). We write 

(6.33) p : ^''t*X -^ T*dX 
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for the projection to the second factor. Note that ^^T*X inherits a natural R-action 
from ''''T*^X, and if ^ e S, T{^f ^ X, then C = p((A - r(0^)"^/^0 e S*dX since 
h = X-T^ on±. 

We also reparametrize generalized broken bicharacteristics 7 satisfying (iii) of 



Corollary IM by letting s = S{t) where S satisfies dS/dt = 2(A - r(7(i))2)i/2^ ^ith 
S'(to) = So picked arbitrarily. We have the following result. 

Proposition 6.8. Suppose that j : I ^ T, is a generalized broken bicharacteristic 
which is disjoint from R+{X) U i?_(A). Then 7 o 5"^^ : J -^ Y,, S defined above, is 
given by 

(6.34) T = \/Acos(s — si), ^" = vAsin(s — si)7a(s) 

where si is an appropriate constant and ^g : J —^ S*dX is a generalized broken 
geodesic, broken at C. IfI = M., then J — {si,si+tt), in particular J has length tt , 
and correspondingly the projection of ^q to dX is a curve of length tt. 

Proof. Let 

(6.35) ,ois)=p( ^^^"(^» ^ 



,VA-r(7(5-i(s)))V 
Condition (ii) of Definition ^^ implies (ii) of Definition |6.6| immediately. Let / G 
C°°{T*dX) be a Trg-invariant function. Let 

(6.36) FiO=fipiiX-T{0')-'/'Oy. 

here we slightly abuse the notation and write p : ^'^Tq-^X -^ T*dX. Then F 



is TT-invariant, so (i) of Definition 6.2 applies and gives d{FT, o ^)/dt{t^). Since 



(Fjr o 7) o S* = /ttq o 7a, the chain rule a nd a short calculation of ^'^HgF gives 



(i) of Definition |6^. The first equation in ( 6.34 ) follows since along 7, ds/dr — 
{ds/dt){dT/dt)-^ ^ -(A - r2)-V2. As (A - r^f'^ = \/Asin(s - si), the second 
equation follows as well. Since r — > =FvA along 7 as t — > ±cx) and r is decreasing, 
we deduce the last statement. D 

It is useful to introduce a relation on S* x S(A) using the structure of the 
generalized broken bicharacteristics given in this proposition. 

Definition 6.9. Suppose ^ G i;(A)\(i?_(A)Ui?+(A)), C G S*dX. We say that ^ --- 
C, if there is a generalized broken bicharacteristic 7 : R -^ ^(A) with 7(^0) = £, such 
that 79 : (a, a + 7r) — > S*dX, as in the above Proposition, satisfies lims^Q+ 7a(s) = 
C- We define ^ ~+ C, similarly by replacing a+ in the limit by (a + tt) — . 

We also need to analyze the uniform behavior of generalized broken bichar- 
acteristics. Here we quote Lebeau's results; they can also be proved completely 



analogously to the proof of Proposition 6.3 given here in the Appendix 



Proposition 6.10. (Lebeau, p7| , Proposition 5] J Suppose that K is a compact sub- 
set ofE, "fn '■ [o-ib] -^ K is a sequence of generalized broken bicharacteristics which 
converge uniformly to 7. Then ^ is a generalized broken bicharacteristic. 

Proposition 6.11. (Lebeau, p^ . Proposition 6] J Suppose that K is a compact sub- 
set of S, [a, 6] C R and 

(6.37) TZ = {generalized broken bicharacteristics 7 : [a,b] — > K{. 

If TZ is not empty then it is compact in the topology of uniform convergence. 
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Corollary 6.12. (Lebeau, |1^, CoroUaire 7]) If ^ : {a,b) -^ R is a generalized 
broken bicharacteristic then 7 extends to [a, b] . 

7. Generalized broken bicharacteristics for totally geodesic C 

We next examine the generalized broken bicharacteristics if all elements of C 
are totally geodesic with respect to h. First we prove that generalized broken 
bicharacteristics 7 : / ^ E with 7(^0) = Co, Co G ^tW H ^'^T^, {Ca, X) are actually 
bicharacteristics of ^^Hg (and hence stay in ^"^T^, {Ca;X)) for t near to. 

Proposition 7.1. Suppose that all elements of C are totally geodesic with respect 
to h. Let 7 : / — > E be a generalized broken bicharacteristic, 

(7.1) 7(io) = Co e (S,(A) n '■-Ta, (a; X)) \ (i?+(A) U i?_(A)). 

Then for t £ J , J a neighborhood of to, we have j{t) £ Sf(A) n ""^T^, {Ca;X), and 
7|j is a bicharacteristic of'^Hg. 

Proof. Our strategy consists of constructing a 7r-invariant function (f) with ^'^Hg(j) > 
c > in a neighborhood of 7r~^(Co)- Thus, by Proposition |6.3| , d/dt{(t)Tr{'j))\t± > 
c > for t € J, J sufficiently small, so (/)7r o 7 is increasing there. This will allow 
us to draw the desired conclusion for the correct choice of 6. We remark that 



this (f) will reappear in the proof of the propagation estimate in Proposition 10.4. 



Moreover, it is essentially the same as the corresponding function in the three-body 



propagation estimate |39, Proposition 15.4], though we will use slightly different 



methods to estimate ^^Hgcj). 

In fact, first we find a 7r-invariant function uj such that ^^Hg will be appropriately 
small near 7t^^(^o)- So introduce coordinates centered at C^ as after Definition 6.2, 
Then the metric function takes the form 

(7.2) ^ = ^ KiniV-, z)lMfJ.j + 2 ^ ft.'i(y, z)fiii^j + Y^ h\\ {y, z)v,Vj 
with 

(7.3) /il(0,0) = l, /i^,(0,z) = 0, 
and, due to the assumption that Ca is totally geodesic, 

(7.4) 9,ftj^,(0,z) = 0. 
We write 

(7.5) h{z, v)=Y h\\ (0, z)viVj 

for the restriction of the tangential part of the metric function to Ca , so 

(7.6) h\y^o = h + Y.hli^{0,z)^i,^,,. 

Now, the Hamilton vector field of h is given by 

(7.7) 

i,j,k ijjrk 'i',j,k 



LU 
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with W ^J2 "j^M, • Hence, ii lj e C°°(R;"-i x M™^) then 

(7.8) H^cj\y=o - Hf^^ + Y.^d,^ e^ - ~h))d-^^- 

k 

Now, /I, hence ft, — /i, is small near 7r~^(^o)i so to model 

(7.9) '^Hg = 2T{^i -d^ + vd,)- 2hdr + Hh, 
we introduce the vector field 

(7.10) W ^ 2T{iy ■ d^) - 2hdr + Hj^ 
locally (near ^o) on ^'^T*{Ca;X). Thus, we have 

(7.11) ''HgUj\y=o ^Wuj- 2{h - h)drLO + Y.^d,,{h - h))d,, 

k 

which is small if Wuj is small. 

We define w as follows. First, Wt — —2ft,, and hzgii^o) ^ since ^o ^ R+W U 
i?_(A), so near £^o, Wt ^ 0, i.e. W is transversal to the hypersurface t = tq. Thus, 
near ^o in ^'^T*{Ca\ X) we can solve the Cauchy problem 

(7.12) WiO^Q, Uj\r=r, = [z - Zof + [u - Uof . 

Since lo and dio vanish at ^O) the same holds on the bicharacteristic of W through 
^0, but <jj > and the Hessian is still positive in directions transversal to the 
bicharacteristics as these hold at ^o- Moreover, by |I|, Lemma 7.7.2], 

(7.13) \duj\ < Cw^/2. 
Let 

(7.14) ro = T^ + k{iy)-X, 

so WrQ — 0. At T = Tq we have tq = hziv) — ftzo(i'o), so 

(7.15) |ro| < C"\duj\ < C"cji/2 

when T — Tq, and then Wuj = = Wtq implies that this inequality holds every- 
where. Therefore, 

(7.16) |ft - ft| < |A - r^ - /i| + |A - r^ - ft| < |A - r^ - /i| + Clo^^^. 
Now, 

(7.17) 

"""Hgu = ""'HgLu -Wuj^-2{h~ h)druj 



+ Yl ^^-^ ('^M (2/, z) - hll (0, z)y,i^jd^,uj. 

Thus, using (7.3)-(7.4), for some C,C' > we have 

rHgC. - Wu;\ < C'{\t^ + ft - A| + ..1/2 + |yp + l^p + \ii\\y\)\div\ 

<C(|T2+ft-A|+C.l/2+|y|2 + |^|2)^l/2^ 



.w 
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Next, note that 



(7.19) 



"Hg\y\^ = 4 ^ K'nM,2/z + 4 ^ /i^V,2/^, 



hj 



«,j 



ri?,|yn<qy|(|y| + H). 



To-T + e-^\y\^ + e-'^uj. 



so by ( [7.3D , 

(7.20) 

For e > let 

(7.21) 
Thus, 
(7.22) 
rHg^ - 2h\ < C{e~'\y\i\y\ + \f,\) + e-W^\\y\^ + \f,\' + \t^ + h - X\ + co'^')). 

We next estimate fi. First, as /i^„(0, 0) — Sij, /i„„ is positive definite in a small 
neighborhood of (0, 0) and 

(7.23) \^i\'<2Y,hli^{y,^)f^^H 



there. On the other hand, 

(7.24) 



1,3 



ij 



(7.25) \J2h^n(.y,z)^^^^^J\<\h-h\+C^\yM+C2\y\'. 

Moving Ci|y||/i| to the right hand side and completing the square gives 

(7.26) Q^\.Cs\y\)'<\h~h\+C^\y\^ 
so 

(7.27) \f,\ < C{\h - h\'/' + \y\), i.e. \f,f < C'{\h - h\ + \yf). 
We can finally estimate ^'^Hgcf), using (7.16) as well: 



(7.28) 



rHg<p-2h\ < C{e-'\y\{\y\+u;'/' + \r' + h-X\^/') 



Note that (f>Tr{^Q) = 0, so near tt ^(^0)7 is small. So now suppose that < 6 < I 
and 

(7.29) < 25 and r - tq < 2S. 
Then 

(7.30) e-^lyp+e-^tj <4(S, 

so \y\ < (46(5)-'^/^, Lj < 4e^5. Hence, under the additional assumption 

(7.31) \T^ + h-X\<e6, 

i.e. that ^ = (y, z, r, /i, ly) sufficiently close to Ea-a, we have 



(7.32) rHg(b - 2h\ < Cie~'{e6y^\e^5) 



V.A^l/2/,2eNl/4 , -2 2 



e €' 



■S) < c's^/\ 



44 ANDRAS VASY 



Since h{'if ^(Co)) > 0, we have /i(^) > 2c > in a neighborhood of tt ^{£,o). Now 
choose So > sufficiently smah, so that C'6q' < c. Note that this requirement 



is independent of e. We thus conclude that for 6 e (0,(5o), ^ satisfying (7.29) and 



( 7.31 ), we have 

(7.33) ''i?s0(O > c> 0. 



Now, using the result of Proposition 6.3, let ^±{t) e I](A — A) be the unique 



points such that Tr{^±{t)) ~ 7(i) and for all 7r-invariant / 
(7.34) (|)(/.o7)|,±=-i/,/(|±(t)). 

Choosing a sufficiently small open interval J around to, r(7(i)), hence T{^±(t)), 



automatically satisfies ( 7.29) for t £ J, while ( 7.31 ) holds automatically as i^±(i) S 



S(A — A). Thus, applying ( 7.34 ) with in place of /, we see that, with 

(7.35) 5(i) = 0,o7(i), 
we have 

(7.36) i G J and g{t) < 26 => ( -^ j |t± > c> 0. 

As g is continuous and .9(io) = 0, this shows that g is increasing on Jn (— oo, to]. To 
see this, first note that g{t) < 26 on Jn(— oo, to], for otherwise g^^({2(5})n(—oo, to]r\ 
J is not empty, g^^{{26}) n (— oo,io] is closed, so taking ti = snp{g^^{{2S}) n 
(— oo,io]) < ^0 and ti e J. Thus, for t G [ti,to], g is differentiablc from cither 
side at t and the derivatives are both positive, so g is increasing on [ti,to], hence 
5(^1) < 5(^0) = contradicting ^(ti) = 26. Thus, g < 26 on J (1 (— oo,0], so g is 
increasing here, so g{t) < for i € Jn (—00, tg). Taking into account the definition 
of 4> we immediately deduce that 

(7.37) Ivilim < Ce^/^ teJn (-(X3, to). 

Since e € (0,1) is arbitrary, we conclude that y{'^{t)) = for t G J n (— oo,io]j 
so 7(i) G '"'T*{Ca;X) for such i. Similarly, w(7(t)) = for such t, so by the 
construction of w, 7(i) is the integral curve of W through ^0 (for t £ J, t < to). 



Of course, a similar argument (with a change of sign in tq — r in (7.21)) works for 
J n [0, 00), so we conclude that 7I j C ^'^T^, {Ca;X) and 7|j is an integral curve of 
W. As W preserves t"^ + h (being essentially its rescaled Hamilton vector field), 
T^(7(^)) + Hlit)) = A, t G J, so 7|j C St (A), and hence at 7r-i(7|j), "=iJg and W 
agree and j\j is a bicharacteristic of ^'^iJg as claimed. D 

Next, we prove that if ^0 G S„(A) n ^'^T^, {Ca, X), 7(^0) = £,0, 7 is a generalized 
broken bicharacteristic, then for a sufficiently small 6 > 0, 7l[o,5] is a generalized 
broken bicharacteristic of A — A, broken at C C C, where C is cleanly intersecting 
and Ca ^ C . This will not use that C is totally geodesic. 

Proposition 7.2. Suppose that £0 € S„(A) n '"^T^/ (C^; X), j is a generalized 

broken geodesic with 7(^0) = ^0 o.iT'd C+ *'S 'JS in Proposition p.jj . Suppose that 
^+ G ^'^T*(Cf,; X)) and & is minimal with this property (i.e. Cc C Cb and £^ £ 
'"'T*{Cc] X) imply c ^ b). Let 

(7.38) C' = C\{C,: C, n C, c CJ. 
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Then for sufficiently small S > 0, 7|[o.d1 's a generalized broken bicharacteristic of 
A - A, broken at C , and 7((0, 5]) is disjoint from '"'T*{Cc] X) if C'c i C . 

Proof. Let b be as above and introduce local coordinates centered at C^. We may 
assume that Cb is given hy y' — for a suitable splitting y = {y' , y"). Thus, ^+ is of 



the form ^+ — (0, 0, tq, 0, /Iq, I'o), and as ^+ S S„(A), /Xq ^ 0. By Proposition 6.3, 
taking into account that y is tt- invariant, 

(7.39) d{y] o 7)/dtk„+ = 0, d{y'; o 7)/^*^+ = (Mo)j- 

Since /Iq 7^ 0, there exist c > 0, Jq > 0, such that |y"(7(i))| > c(i — ip) for 
t £ {t(),to + So), while for any e > there exists ^1 > such that \y'{^{t))\ < e{t — to) 
for t £ {to, to + Si). In particular, for any e > there exists 6 > such that for 
t £ {to, to + S) we have \y'{l{t))\/\y"{'y{t))\ < e. By choosing e > sufficiently 
small we can thus make sure that 7(t) ^ Cc for t £ {to,to + S] ii Cc ^ C . Hence, 
7l[to,to+i5] <^^^ t)e regarded as a curve in ^c'aec^'^TQ, {Cc, X), C'^ taken with respect 

to C, if we let 7(^0) = 7rofc(Co) € ^^T*{Cb] X). Of course, ^\(to.ta+s\ '■s ^ generalized 
broken bi chara cteristic , broken at C (since it has no points above C\C'). Thus, by 
Corollary |6.12| , "i\{to.ta+s] extends to a generalized broken bicharacteristic, broken 
at C , defined on [to, to + 8\; by continuity of 7 this must coincide with 7, so 7 is a 
generalized broken bicharacteristic, broken at C , as claimed. D 

We can combine the previous results to deduce the structure of the generalized 
broken bicharacteristics if C is totally geodesic. 

Proposition 7.3. Suppose that C is totally geodesic with respect to h and j is a 
generalized broken bicharacteristic, broken at C with ^0 ~ 7(^0) G ""^T^, {Ca',X). 
Then there exists 5 > Q such that both ')\[tQ to+S) o.nd ^\{tQ~s to] ^'"6 bicharacteristics 

Of'^Hg. 



Proof. If ^0 G i?+(A) U -R-(A) then ^{t) — S,o for t near to by Proposition 6.5, 
hence near to, 7 is a (7r-projected) bicharacteristic of ^^Hg (as ^'^ Hg vanishes at 
R+{\)\jR-{\)). If Co e Sf(A)\(i?+(A)Ui?_(A)) then Proposition O applies and 
proves the result. If ^o G S„(A), then with C as in Proposition 7^, ^\[o.s) is a 



generalized broken bicharacteristic, broken at C, with 7(^0) G ^^T%, {Cb; X) n St (A) 



(prime taken with respect to C). Thus, Proposition 7.1 applies again and proves 



the result. D 

A compactness argument gives at once 

Corollary 7.4. 7/7 : [a,b] —f'Eisa generalized broken bicharacteristic, broken at 
C, and C is totally geodesic, then there exist to — a < ti < t2 < . . . < tm — b such 
that 7|rij,tj_|_i] is a bicharacteristic of A — X (i.e. it is not broken). 

8. Positive operators 

In the following two sections we discuss technical points of the microlocal positive 
commutators constructions. In this section we show roughly speaking that the 
positivity of the indicial operators of A S ^^g^' {X,C) implies the positivity of A 
modulo compact operators. We prove this by constructing an approximate square 
root of A. In the next section we examine commutators [A, H] in more detail. 

Throughout this section we assume that H is a many-body Hamiltonian. We 
start with the basic square root construction. 
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Lemma 8.1. Suppose that H is a many-body Hamiltonian and A e M. Suppose 
also that A £ ^^c^' {X,C) is self-adjoint, and for some c > and ^j G C^(R) 
which is identically 1 near A, 

(8.1) ^{H)A^:{H) > c^{Hf. 
Then for any c' G (0, c) and (j) G C^(M) such that 

(8.2) supp0nsupp(l- V) =0, 
there exists B G ^^g^' iX,C) such that 

(8.3) 4>iH){A - c')(l){H) = (t){H)B*B(l){H). 
Proof. Let 

(8.4) P = ^{H)A^{H) + c{ld~i;{Hf) G *^;°(X,C). 

Thus, P > c, so P — c' > c — c' > 0. Since the spectrum of P — c' is a subset of 
[c - c', cx)) and c - c' > 0, we have (P - c'f/'^ = /(P - c') where / G C^(R) and 
f{t) = \ft if i is in the spectrum of P — c'. By Proposition [4.9^ 

(8.5) = (P - dfl^ = f{P - c') G vl/gf (X,C). 

Let -01 be identically 1 near suppc/) and vanish near supp(l — "ip). Then 

(8.6) MH)Q^MH) - MH){P - c)MH) = MH){A - c')MH). 
Now let (p G C^(R) be identically 1 near A and vanish near supp(l — ?Ai). Let 

(8.7) B^QMH)e^sc"'°{X,C)- 



Multiplying (|8.6[) from both sides by 4>{H) then proves (S.3). D 



We now show that under certain additional assumptions, the positivity of the 
indicial operators implies positivity of the operator modulo lower order (hence com- 
pact) terms in the calculus. We start by assuming strict positivity of the indicial 
operators when localized in the spectrum of H. 

Proposition 8.2. Suppose that H is a many-body Hamiltonian and A G M. Sup- 
pose also that A,C E '^s^' i^i^) '^'"6 self-adjoint and Ca.o{C) = Ca((^)'(/'o(^a(C))^ 
for every a and C G '^'^T* {Ca', X) where CaiC) 's o- function with Ca{C) > 0, f/'o = 1 
near A G K, 'i/'o G C^(R). Assume in addition that there exists ip G C^(R) which 
is identically 1 near \, supp V' H supp(l — ipo) = 0, such that 

(8.8) V(i?a(C))ia(C)^(i?a(C)) > V(i?a(C))Ca (OVX^a (0) 

for every a and ( e "'T*{Ca]X). Then for any e G (0, 1) and G Cf^(M) with 

(8.9) supp (j) n supp(l - V) = 0, 
there exists R G ^g^' {X,C) such that 

(8.10) (l){H)A(j){H) > (1 - e)(j){H)C<j){H) + R. 

Proof. We apply a parameter dependent version of the previous lemma to the in- 
dicial operators to conclude that for each Q there exists -Ba(C) with 

(8.11) 

0(7J,(C))(ia(C) - (1 - e)Ca{C))cP{Ha{C)) - </'(i^a (C))i?a (0* 4 (C)0(i?a (C))- 
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It follows from the Cauchy integral formula construction of the square root in the 



calculus and the explicit formulae (B.4), (Kj) and (p.7D that the indicial operators 



Ba{C) match up so that there exists B G ^sc°°' i-^^ ^) '^ith indicial operators Ba{C)- 
Here note that the set where ip{Ha{C)) does not vanish has compact closure, hence 
c is bounded below on it by a positive constant. Thus, we can take the same smooth 
function / in the expression ( |8.5| ) for the square root for every a and (. By ( ^.11 ), 

(8.12) <P{H){A - (1 - e)C)(p{H) = (j>{H)B* B<j>{H) + R 

with R e ^~^'^{X,C). Since 4>{H)B* B4>{H) > 0, rearranging this proves the 
proposition. D 

Similar conclusions hold if we assume a two-sided estimate on the indicial oper- 
ators of A. In essence, this forces the indicial operators, hence their square roots, 
to vanish to infinite order when c vanishes. 

Proposition 8.3. Suppose that H is a many-body Hamiltonian and A G R. Sup- 
pose also that A,C G ^5^' iX,C) are self-adjoint and Ca,oiC) = Ca(C)'0o(^a(C))^ 
for every a and (^ G "'^T*{Ca', X) where Ca{C) '* a function with Ca{C) ^ which van- 
ishes with all derivatives at each (^ with Ca{C) ^ 0, tpo = 1 near A G M, V'o G C|?°(R), 
^a(C) = */ca(C) = Oi and for any differential operator Q G Diff(^'^T*(Ca; ^)), 
all semmorms 0/ g(c,(C)-iia(C)) m *^^°°'''(p-i(p), TpC"), C e ^%*(C„;X), are 
uniformly hounded on the set of C,'s with Ca{C,) > 0. (This is almost, but not quite, a 
statement about the seminorms of Ca{C)^^Aa{C) "in ^1^°°'^ {p*g^'"^T*(Ca]X),Ca), be- 
cause we restrict our attention to the region where Ca(C) > 0, and do so uniformly.) 
Assume in addition that there exists ip G C^(R) which is identically 1 near X, 
swppip ^ supp(l — ?Ao) = 0; such that 

(8.13) ^(i?„(C))ia(C)^(^a(C)) > ^(i?a(C))Ca(C)V'(^a(C)) 

for every a and (. G ''"^T*{Ca', X). Then the conclusion of the previous proposition 
holds, i.e. for any e G (0, 1) and (j) G C^(R) with 

(8.14) supp0nsupp(l- V') =0, 

there exists R G ^5^' {X,C), with seminorms bounded by those of A and C in 
*s;f°'°(X,C), and with WF's^iR) C WF's^{A) U WF's^iC) such that 

(8.15) (l){H)A(j){H) > (1 - e)(j){H)C(j){H) + R. 

Proof. We define Ba{C) = if Ca{Q = 0, otherwise we define Ba{C) as in the 
previous proposition. The only additional ingredient is the analysis of i?a(C) near 
( with Ca(C) = 0. To do this analysis, we follow the construction of Ba{C) in detail. 
So let 

(8.16) Pa(C) = V'(i?a(C))ia(C)V'(i?a(C)) +C,(C)(Id-V(^a(C))'), 

and let 

(8.17) cl(C)-(l-6)c,(C). 
Thus, A(C)-<(C)>eCa(C). Let 

(8.18) QaiO = (PaiO - C^iOf'^ = Ca{Cf'^{c,{C)-'Pa{C) - (1 - e)f'\ 
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By our assumption, there exists M > such that the norm of Pa{C) in B{L^,L^) 
is bounded by Mca(C)- Now choose / e C^(R) such that f{t) = Vt on [1 - e, M]. 
Then M > Ca{C)-^Pa{C) - 1 + e > e, so 

(8.19) QaiO - CaiCy^^iCaiO-'PaiC) - (1 - e)) . 

By our assumptions, the seminorms of CaiO^^PaiC) in ^sc (Pa ^(-P)i ^p^")' C ^ 
^^T*{Ca] X)^ remain uniformly bounded as Ca{C,) —^ 0, so the Cauchy integral repre- 
sentation of/, via an almost analytic extension, shows that f {ca{C)^^ Pa{C) ^ (^ ^ ^)) 
remains uniformly bounded. Thus, Qa{C) is continuous as a function on ^'^T* (Ca, X) 
with values in \&g^ {p~^{p), TpC""). A similar argument also holds for the derivatives 
of QaiO- Lfit ipi be identically 1 near suppi/) and vanish near supp(l — ip), and let 

(8.20) BaiO^QaiOMH). 

Again, the -Ba(C) match up so there exists B G ^gc°°' iX,C) with these indicial 
operators. We can also make sure that the lower order terms also vanish where c 
does, i.e. that WFgp(_B) C suppc. Then the indicial operators of (p{H){A — (1 — 
e)C)4'{H) and (j){H)B* B(t>{H) are the same, so 

(8.21) (l){H){A - (1 - e)C)(j){H) = (j}{H)B* B4){H) + R 

-ex 

Sc 

9. Commutators 



with R e ^g(°°' {X,C), proving the proposition. D 



In this section we discuss the basic technical tool underlying the propagation esti- 
mates of the following sections. Thus, we show how an estimate of the commutator 
[A^ H] at Co, which is essentially obtained by a symbolic calculation in the scattering 
calculus, can give a positive commutator estimate under the additional assumption 
that HafiiO has no L^ eigenfunctions for any a ^ and ( G ^'^T*{Ca',X). In the 
Euclidean setting this means simply that the subsystems have no bound states. 

To do so, we extend the notion of a function being 7r-invariant to functions on 
scrp*-^ in a trivial way: q G (joo^scy*^^ jg 7r-invariant if ql^cx* x is 7r-invariant. 
Since the analysis of classical dynamics, i.e. of generalized broken bicharacteristics 
of A — A, broken at C, is based on the properties of 7r-invariant functions, we will be 
interested in quantizing 7r-invariant symbols. More specifically, we are essentially 
interested in operators of the form A = Qipo{H), ipo £ C^(R), where Q is obtained 
by quantizing a 7r-invariant function q € C°°{^'^T*X). Since such Q would not be 
in our calculus, we construct A directly. 

All considerations in what follows will be local, i.e. we will assume that the 
projection of the support of g to X lies near a fixed p € dX, so we can always work in 
local coordinates and identify X with §" . The problem with such q E C°°(S" x R") 
is that they are rarely in C°°(§" x §"), i.e. they are not symbols in ^, so Q will not 
be in ^g^°(§") or indeed in $g|, (§",C). This, however, is not a major difficulty. 
Fix ipQ e C^(K; [0, 1]) which is identically 1 in a neighborhood of a fixed A. Thus, 
ipo{H) G *sc°°' (^i'^)i so it is smoothing. At the symbol level, ipo{H) is locally the 
right quantization of some 

(9.1) peC°^i[Sl-X]xSl) 

which vanishes to infinite order at [S" ; C] x 9S" , which will enable us to write down 
A directly. 
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We are thus interested in the following class of symbols q. We assume that 
q e C°°(R'|, X R?) and that for every multiindex a, (3 E N" there exist constants 
Ca^is and ma,f3 such that 

(9.2) \{DZD^^q){w,0\ < C„,^(u')-l"l(0'"-^- 
This implies, in particular, that 

(9.3) ge^°(§'; xM"), 

i.e. that g is a 0th order symbol in w, though it may blow up polynomially in ^. 
Indeed, in the compactified notation, ( |9.2| ) becomes that for every P S Diffb(§"), 
acting in the base {w) variables, and for every /3 £ N" there exist Cp^p and mp^p 
such that 

(9.4) |(PZ)f)Q|<Cp,^(0"-■^ 

It is convenient to require that q be polyhomogeneous on §!{: x M": 

(9.5) geC°°(§!;: xM"); 

this stronger statement automatically holds for the vr-invariant symbols we are 
interested in. 

We next introduce the product symbol 

(9.6) a{w,w',O^Qiw,Opiw\0, 

where ipo{H) is given locally by the right quantization of p. The main point is 



Lemma 9.1. The symbol a defined by (9.6) is in C°°{Sl x [§l;C] x ED and it 



vanishes with all derivatives at [§" ; C] x 9§" . Hence, it defines an operator A G 
^Sc^' i-^^^) ^y ^^^ oscillatory integral ( p. 13 ). 



Proofi First, a G C°°(S!;. x [E.l;C] x M") follows from (|t]) and (jgj). Moreover, 
the infinite order vanishing of p at [S";C] x §"-~^ implies that for every P' G 

Diffb([S!^;C]), /? e N" and A^ e N, 

(9.7) \P'dIp\ < Cp,,0,^(e)-^. 

Thus, Leibniz' rule shows that for P £ Diffb(§!^) acting in w, P' E Diff b ( [§!^ ; C] ) 
acting in w' , 13 E N" and N 

(9.8) \PP'Dla\ < Cpp^^,jv(0"'^. 

But this means precisely that a E C°°(S" x [S";C] x §" ) and it vanishes to infinite 
order at the boundary in the last factor. D 

The indicial operators of A are just given by the quantization of the appropriate 
restriction of a similarly to ( 4.40[ ) (except that now a depends on the base variables 



from both the left and the right factors of S!f:). This takes a particularly simple 



form if q is 7r-invariant, for then, in the notation of (4.40), q is independent of both 
Y and ^". Thus, we can take q outside the integral in ( 4.40 ), i.e. it simply multiplies 
the indicial operator of ■ipa{H) by a constant. 

Lemma 9.2. Suppose that q E C°°('"^T*S") is -K-invariant and it satisfies (p!4|). 
Let A E '^ sc°' i-^T^) ^6 as in the previous lemma. If ( E '^'^T*{C'a',X), then 

AaiO = q{c)MH),{C). 



Combining this lemma with Proposition 5.3 gives 
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Corollary 9.3. Suppose that ( e '"T^, {Ca;X) and u e C"°°(X). If A is as m 



Lemma 9.i, q{Q ^ 0, Au e C°^\X) and C i WFsdiH - X)u) then C i WFsc(w). 



Since the indicial operator of [A, ii\ ~ AH — HA in 'I'g(°°' {X,C) is just 
(9.9) [AHloiO = [AaAO,HaAO]=qiO[MHaAO),HaAO]=0 



for every a and C G ^'^T*{Ca'TX), we see that for every A as in Lemma 9.1, 
[A, H] E ^gc°°' (^)C). The additional order of decay corresponds to the one in 
the scattering calculus. Moreover, the indicial operator of [A, H] at Co, as an oper- 
ator in ^g^f"' {X,C) (so this indicial operator is just a function on '^'^T* {Cq; X)) , is 
given by the Poisson bracket formula from the scattering calculus. Since V vanishes 
at Co, this gives 

(9.10) i[AH],,o = -^'HgiqMg)) = -MgrHgq. 

If the indicial operators of H at the other faces have no L^ eigenfunctions, then this 
estimate combined with a compactness argument suffice to prove an estimate for 
[A, H] modulo lower operators (i.e. modulo ^gj^' {X,C)). However, to make the 

compactness argument work, we need to estimate the indicial operators, [A, H]^ ^, 
for all a. This is facilitated by the following lemma. 



Lemma 9.4. Let q and A be as in Lemma \9.^ . For every seminorm in 

and for every I G N there exist C > and m £ N such that for every a and every 
C e ''T;{Ca;X), p e Ca, the semmorm of [AJI]^AO ^n ^-^'"ip-^{p),TpC'') is 
bounded by 

(9.11) c(k(c)i+ Y. ^^p\in-\d^.dq){c,n\) 

\a\<m ^° 

where the differential dq is taken with respect to all variables, in '><^T*S^^, i.e. it is 
the differential of q C C°°C"=T*S![). 



-00,0 

which can be seen directly from our calculations in the following proof. 



Remark 9.5. Similar conclusions hold for every seminorm in ^°°|, {p*^^'^T*{Ca',X),Ca), 



Proof. This can be proved directly from the definition of the indicial operators, 
i.e. by computing x~^e~'f[A,H]e'fu' where / G C°°{X) and u' e C°°{[X;C]), 
similarly to |g^, Sections 7,13]. Since this is equal to x~'^[e~''^ Ae^^ ,e~^f He^^u' , 
and e~^^ Ae^f e ^sc°°' (^'^)j ^^ can assume that / = 0, the calculation being very 
similar in the general case. To compute the commutator, it suffices to commute 
both Av and Hv for every v G C°°([X;C]) modulo terms that vanish with their 
first derivatives in Pg^Ca- A straight-forward calculation can be performed just as 
in (4.34)-( Oc| ), where only the 0th order terms were kept. That shows with our 
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coordinates that 
(9.12) 

[Xi?]a,l(C)=[(&U(C),i?a,0(C)] 

+ {-{Drqm{[Y,HaAO]dY+YdYHaAC)) 

+ {D,q){C){dA..o){0 - {d,q){C){D,Ha,om 

+ {drqWil^ ■ D^Haflm - {y ■ D,q){0{drHafl))4'{HaA0)- 

Here dxA denotes the operator with kernel given by dx apphed to that of A. Since 
in our notation the kernel of A is 

(9.13) / e*("'-"'')-«g(u;, S)pW .0 d^, 

with the integral being convergent, rewriting this with the coordinates on the com- 



pactification [§";Ca], (2.7), so that q takes the form q{x,xY, z,S,) proves that all 
terms of (9.12) satisfy the stated estimate, completing the proof. 

Another approach to compute a-indicial operators is to use that near C^, A 
can be regarded as a (non-classical!) pseudo-differential operator in the free vari- 
ables (waj^a) with values in bounded operators on L'^{Xa) (in fact, with values in 
^g^'^iX'^X"))- More precisely, A e ^^^^''°{Xa;B{L^{X''),L'^{X''))). This allows 
us to use the scattering calculus for the computation of the commutators to give 
the stated result. D 

As an application of these estimates, we now show how, under the assumption 
that the subsystems have no bound states, a positive Poisson bracket with g can 
give rise of a positive operator estimate. We thus assume that 

(9.14) HaAO has no L^ eigenvalues for any a 7^ and C e '"'T*{Ca;X). 

To simplify the notation in the following proposition, we introduce the notation 
supple C """"T* {Ca, X) for 7r-invariant functions e e C°°{'"'Tgx^)- This is defined 
as the support of the function on ^'^T* {Ca, X) induced by e. Indeed, as e is tt- 
invariant, its restriction to ^"^Tq X can be regarded as a function on ^'^T*{Ca;X). 
Then supp^, e is the support of the pull-back of this function to ^'^T*{Ca\ X). 



Proposition 9.6. Suppose that H is a many-body Hamiltonian satisfying ( p.l4 ), 



and A G M. Suppose also that q,b,e € C°°('"^T*X; R) are ir-invariant, satisfy the 



bounds (9.4), q,b > 0, and that there exist S > 0, C > 0, Cq > 0, such that for all 



(9.15) |5(C) - A| < <5 ^ ^^H.qiO < -biCf + eiO 
and 

(9.16) 1.9(0 -X\<S and£_i suppe => q{0 < Cb{£,f and \{d'^dq){0\ < Ca,b{S,f. 

Let A G 4*^^°"' {^tC) be as in Lemma |ff. j| . For any e' > 0, a £ I and for any 
Ka C '^'^T* {C a', X) compact with supple f] K = ^ there exists 5' > such that if 
ip G C^(M) is supported in (A — S' , X + S') and ( e Ka then 

(9.17) t{ij{H)[A^H]i:{H))^AC) > (2 - e')b\^{HaAC)? ■ 
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Proof. Note that the estimate (9.17) is trivial if r^ + \i'a\z 
{za,Ta, J^a), S' < 1 arbitrary) since then both sides vanish as 

(9.18) m)aiO = i^{ha{Za) + rf + \Va\lj, 



> A + 1 (with C, 



ha denoting the subsystem Hamiltonian as in (S.E), and H^ > by the assumption 
on the absence of bound states of all subsystem Hamiltonians (including H^. with 

Ca C Ce). 

We prove (9.17) by induction on a. First, (9.17) is certainly satisfied for a — Q. 
In fact, as A € ^sc°°' (X,C), we can use the commutator formula in the scattering 
calculus, ( |9.10D , to find [A^H]^^ ^. Since V vanishes at the free face, /3g^Co, it does 
not contribute to [A, _//](, j^, so we indeed have, by ( |9.15 ), 

(9.19) zV(i/)[i^^A>]^(i/)o_i - -2q{--Hgq)ij{gf > 2b\^{gf = 2b\ij{Ho,of 
away from suppp e under the assumption that 

(9.20) suppV'C (A-(5,A + J). 



So suppose now that ( 9.17 ) has been proved for all c with Ca C Cc, Ca ^ Cc- This 

implies that all indicial operators of i[ip{H)A*A^{H),H]^{(), ( = {za,Ta,i^a) € Ka 
satisfy an inequality like (9.17). In fact, the indicial operators are of the form 

t[i,{H)A^{H),H]^{C) with /3sc(C) = i0,za) e Ca, ^ca(C) = C. Such a C is of the 
form C = {Y^',Za,Ta,fia,iya) whcrc Cc is given hy x = 0, y' = 0, so {Y^',Za) give 
coordinates along Cc- Note that as Ka is compact, so is 

(9.21) 

Kc = {C = (Y;:, Za, Ta,til, ^a) ■ {Za,Ta, ^a) G Ka, PsdO € C„, \^ll\ < A + 1} 

and as e is independent of /i" at Ca, i^'cH supp^ e = 0, so we can apply the inductive 
hypothesis. Taking into account that the estimate ( ^.17 ) is trivial at Cc for C, with 
\^.'a\ > A+ 1, we see that for all C = (0, Za, Ta, [tl'a, Va) with {za, Ta,Va) e Ka, we have 



(9.22) 



^{^{H){A*A,HM,{H))^^^{C) > {2 - e')b\ij{H)\^,{0. 



Since b^q is 7r-invariant on ^'^T*X, it is independent of C, for each fixed C, and if it 
vanishes at C, then so does [4>{H)A*A^{H), H]^^j^{C) by Lemmas |l[fj and ( |9.16D . 
Thus, by Proposition 8.2 , 

(9.23) imH)A^{H),Hl^,{C) > (2 - e' )b\^KH)l, oiO + ^(0 
where the seminorms of 

RiO e ^sr^\pa\p),Tpn, C e ^%*(Ca;X), 

are bounded by those of [ip{H) A* Ail}{H) , H]^ ^{() and by b{()'^q{(). By assumption 
( |9.16| ) and Lemma 9A the former are bounded by the latter, so i?(C) satisfies the 
estimate 

(9.24) \\RiO\\BiLUPl\p)),Hh\pl\p))) < C'qiCm? 
with C" independent of q and b. 
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We now use our hypothesis on the absence of bound states. So suppose that 
ipi,ip2 € C^ (M) , f/; = 1 near supp V'l , V'l = 1 near supp -02 ■ By assumption, 
X — T^ ~ Ivall is not an eigenvalue of the subsystem Hamihonian, ha{z). Thus, 

(9.25) MHaiO) = MKiz) + tI + \vJ^J ^ 
strongly as supp'i/'i — > {A}. Since Ka is compact, and the inclusion map 

(9.26) T:}ll-\p-\v))^Ll^{p-\v)) 
is compact, for ^/^i with sufficiently small support we have 

(9.27) m'^'^)a{0\WHh\,-\,)),LUp-.\m ^ ^'(^")"' 

for all i^Ka. Thus, 

(9.28) *(^i(iI)[i^iJ]7Ai(if))^ i(C) > (2 - ^)h\Mii\^{C,) - e'62g, ( G K^. 
Multiplying by ip2iH) from both left and right we finally conclude that 

(9.29) tiMH)[A^H]MH)),^, > (2 - 2e')b^qMH)lo- 

Relabelling ip2 and 2e' as ip and e' (thereby putting stronger restrictions on supp ■0) 
provides the inductive step and completes the proof of (9.17). D 

In the following corollary we add an extra term to the commutator that will 
enable us to deal with other terms arising in the propagation estimates. 



Corollary 9.7. Suppose that the assumptions of Proposition \9. (\ are satisfied and 
let C be as in ( p.lq ) . Suppose in addition that for any differential operator Q on 
^'^T*(Ca',X) and multiindex a there exist constant Cq and Ca,Q such that 

|5(0 - A| < (5, 5(0 7^ and S, ^ suppe 
^^■^°^ ^ \Q{h-\m\ < Cq and \Q{b-\d^dq)m\ < C^,Q. 

For any e' > 0, M > 0, and for any K C ^'^T* X compact with suppeflif = there 
exists S' >0, B,Ee ■^^^'"{X,C), F e ■^^^'\X,C) with 

(9.31) 

WF'g^iE) n i^ = 0, WF's,(F) C suppg, Ba,o(C) = &(C)g(C)'/'V'(^a,o(C)), C e K, 

such that if ip (z C^IR) is supported in (A — (5', A + S') then 

(9.32) 

i^P{H)x-^/^[A* A, H]x-^/^i:{H) - Mil){H)A* Ail){H) > (2 - e' - MC)B*B + E + F. 



Proof. Let p e qoo ^scrp* -^^^ ^^ 7r-invariant, p > 0, satisfy estimates ( |9.4[ ), and such 
that suppp n suppe ~ and supp(l — p) K = 9- (Here p can be regarded as a 
function on "'=T*X.) Let ipo G Cf (R; [0, 1]) identically 1 near [A - (5, A + (5], and let 
P e *s,°°'°(X,C) be such that Pa,o(C) = PiOMHaiO) and WF's,{MH) - P) n 
K = %. For example, P can be constructed as in Lemma |9.l| . 
The indicial operators of 

iiP{H)P*x-^^^[A* A, H]x-^/'^ Pil;{H) ~ M%Ij{H)A* A%Ij(H) 



54 ANDRAS VASY 

are 
(9.33) 

iijiH)P*x-^/'^[A*A^H]x-^/^P'iJj{H)^_giC) - M{il;{H)P*A*APtP{H))^^„iC) 

since V'oV' = "0- Thus, by Proposition |9.6| and as Mq < MCb^, we have 
(9.34) 

iiP{H)P* x-^l'^[A*A^H]x-^l'^P'ijj{H) ^^C) - M{il;{H)P*A*APil;{H))^AC) 

>{2-e'-MC)b\ij{Ha.o{C)f. 



Thu s, ta king into account ( 9.30 ) and the remark following Lemma |9.4| , Proposi- 
tion ^^ gives 

(9.35) 

ii}{H)P*x-^l'^[A* A,H]x-'^/'^Pil:{H) ~ Mil:{H)A* Ai:{H) > {2-e')B*B + F, 

with B e ^sr^'{X,C), F e ^sr-\x,C), 

(9.36) KoiC) = PiOKCMCV^', 



so the second statement of ( 9.31 ) holds. Moreover, writing 'iJj{H) — Ptp{H) + 
(ipoiH) — P)iMH), and expanding the left hand side of ( |9.32|) , every term but the 



one given in (9.35) has operator wave front set disjoint from K. Letting E be the 



sum of these terms proves the corollary. D 

10. Propagation of singularities 

In this section we prove that singularities of generalized eigenfunctions of the 
many-body operator H propagate along generalized broken bicharacteristics under 
the assumption that that no (proper) subsystems of H have a bound state. That 
is, due to our definition in Section |6|, we assume that 

(10.1) HbfiiO has no L^ eigenvalues for any & 7^ and ^ e ''''T*{Cb; X). 



The technical reason for this assumption lies in the argument of Proposition ^ 

in which a symbolic estimate is used to deduce positivity estimates for the indicial 

operators. However, it is clear that the generalized broken bicharacteristics of A — A 

cannot be expected to describe propagation if the subsystems have bound states 

since in this situation even the characteristic set of H (i.e. the set where i?6.o(C) is 

not invertible) changes. 

Suppose that p £ C'^^ ~ C (the regular part of C). As in Section ^ let {x, y, z) = 

{xa^ Ha, Za) be Coordinates on X near p with x defining dX as usual, C defined by 

X = 0, y = 0, chosen so that every C'b with p g C;, is a product-linear submanifold 

of dX in these local coordinates, i.e. it is of the form {{y,z) : Aby — 0} where 

A = Ab is a matrix. In addition, as in Section |6|, we arrange that at C, dy- — 9(y„) 

is perpendicular to TC for each j (with respect to h) and they are orthonormal with 

respect to each other at p. Let (r, ^, j/) = {Ta, fJ-a, Za) denote the sc-dual variables, 

so we write elements of ^'^T*X as 

dx dy dz 

10.2 T^+^i-^ + v—. 

X'^ X X 
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Thus, 'dtp (i.e. on^'^T*X) the metric function oih is of the form |/ip+/i(z, ly) with \^\ 
denoting the Euchdean length of /i and h is the metric function of the restriction of 
h to TC. When talking about C},, we sometimes write the corresponding orthogonal 
sphtting of 2/ as y = (y', y"), so Cj, is defined by At,y = y' = in dX. 

Recall that TTpa : """T^X -^ '''=r*(C;X) is the (orthogonal) projection given 
by the metric at C. Thus, in our local coordinates {y, z,t, fj,,i>) on ^'^Tg-^X, 
7roa(0, z, r, /i, J/) = {z,T,iy). We use composition with the projection ^'^Tg-^X to 
^^TqX given by our choice of local coordinates, (y, z, r, /x, v) i— s- (z, r, /i, v), to extend 
TToa to a map, denoted by 7rg^, from ^'^Tg^X to ^^T*{C; X). Thus, vrg^(y, z, r, /i, t^) = 
(z,r,i/). 

The propagation of singularities estimate in directions tangential to C proceeds 
much as in the 3-body case. In fact, essentially the same operator as there gives a 



positive commutator, see Propositions 10.4-10.5; the functional analysis part of the 



argument is much as in the normal case which we proceed to examine. Recall that 
the normal part of the characteristic set of iJ — A over C" is 

(10.3) S„(A)n"'=T5,(C;X) = {(z,r,i/) : r^ +h{z,v) < \}. 

Since the characteristic set Sa-a of A — A is given by r^ + \v\^^ + |^p = A at p, 
the condition 7r(|) € S„(A) n "'=Tp*(C;X), | e Ea-a implies that /^ ^ 0. Since the 
rescaled Hamilton vector field ^^Hg of A (restricted to ^'^Tg^X) is given by 

(10.4) ^^ffg = 2t{ii ■d^. + y d,) - 2hdr + H^, 

the dy component of ^'^Hg at pis 2fi-dy, meaning that bicharacteristics of A through 
^ are normal to ^'^TqX. In addition, with rj — y ■ fi, rj is 7r-invariant and can be 
used to parameterize bicharacteristic curves near ^ = 7r(^). In fact, at each Cb 
with p Cz Cb, T] = fi ■ y has the property that if we split y = {y' , y") so that a; = 0, 
y' = defines Cb then fj, ■ y — fi' ■ y' + fi" ■ y" is independent of fi' at y' = 0, so 77 is 
TT-invariant. Moreover, ^'^Hgri{^) = 2|/ip > 0, so 77 can be used to parametrize the 
generalized broken bicharacteristics near ^ as claimed. We remark that r is another 
possible variable to use for the parameterization, as usual. 

We now proceed to prove two normal propagation estimates. The first one will 
be less precise, but it works under our most general assumptions. On the other 
hand, the second estimate requires that all elements of C be totally geodesic, but 
it locates the incoming singularities more precisely. Although the consequences 
are the same, as far as propagation along generalized broken bicharacteristics is 
concerned (due to the geometry of these bicharacteristics), the finer estimate is 
worth proving since it is closer to the tangential estimates in spirit and it applies 
in the setting of most interest, Euclidean many-body scattering. 

We only state the following propagation result for propagation in the forward 
direction along the generalized broken bicharacteristics. A similar result holds in 
the backward direction, i.e. if we replace r]{S^) < by 77(^) > in ( 10. 5| ); the proof 



in this case only requires changes in some signs in the argument given below. 



Proposition 10.1. Suppose that H is a many-body Hamiltonian satisfying (10.1). 
LetueC^°°{X), A >0. Let^o = {zo,to,vo) S S„(A)n"'=T5,(C; X) and let r] = y ■ ft 
be the TT-invariant function defined in the local coordinates discussed above. If there 
exists a neighborhood U of ^0 *^ E such that 

(10.5) ^£U and f]{0 < ^ ^ ^ WFs,(m) 
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Remark 10.2. Note that r/(^) < implies y ^ 0, so (, (^ "'^T(^,(C;X). 

Proof. The main step in the proof is the construction of an operator which has a 
microlocally positive commutator with H near ^o- In fact, we construct the symbol 
of this operator. This symbol will not be a scattering symbol, i.e. it will not be in 
C°°(S" X §" ), only due to its behavior as ^ ^ oo corresponding to its 7r-invariance. 
This will be accommodated by composing its quantization with a cutoff in the 
spectrum of H, 4>{H)^ (j) G C^(R) supported near A, as discussed in Lemma |9.1[ 
This approach simply extends the one taken in the three-body scattering proof 
of [p9|, though the actual construction is different due to the more complicated 
geometry. 

Employing an iterative argument as usual, we may assume that ^q ^ WFgj, (u) 

and we need to show that ^o ^ WFg^ ' (u). 

First we define a distance function to ^o- Thus, we let 

(10.6) CO = |y|2 + 1^ - zoP + |t - Top + 1^. - j/oP, 

|.| denoting the Euclidean norm. Then uj vanishes quadratically at ^q, so \duj\ < 
C'lLO^/"^. In particular, 

(10.7) rHgUj\ < CiLU^/^. 

Next, we use the variable ri — y ■ fi to measure propagation. Let 

(10.8) co = A-r2-|^.o|^„>0. 

Since the dy component of ^'^Hg at (0, zq, t, fi, v) is 2/i, we see that 

(10.9) rHgf^ - 2\fi\^\ < C'^{\y\ + \z- zqD < C2UJ^^^. 
In addition, 

\>^-r^-Wo\l-\^^n<\x-9\ + \9~T^~Wo\l-\^^^\'\ 

<\X-9\ + C\\y\ + \z-z„\ + \T-To\ + \u-u,\)<\X-g\+C3u;^/^ 
so we conclude that 

(10.11) rHgTj - 2co| < C4(|A -g\+ uj^^^). 

For (3 > 0, 6 > 0, with other restrictions to be imposed later on, let 

(10.12) <j> = v+^uj, 

so is a TT-invariant function. Let xo G C°°(M) be equal to on (— oo,0] and 
Xa{t) = exp(-l/t) for t > 0. Thus, Xo(i) = t~^Xo{t). Let xi e C°°{^) be on 
(—00,0], 1 on [l,oo), with x'l > satisfying x'l G C?°((0, 1))- Furthermore, for 
Aq > large, to be determined, let 

(10.13) q = xo{Ao\2 ~ 4'/d))xi{y ■ ii/5 + 2). 

Thus, on suppq we have (j) < 26 and y ■ fJ. > —25. Since uj > 0, the first of these 
inequalities implies that y ■ fJ. < 26, so on suppq 

(10.14) \yf^\<2S. 
Hence, 

(10.15) uj<{S/P)(26-y-n)<4:6'^p-\ 
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We now proceed to estimate ^'^Hgcj). First, by ( 10.11 ) and ( |10.7| ), 

(10.16) rHgC^ - 2co| < C4(|A -g\+ u:^'^) + ^uj^'\ 

So let 



(10.20) 



' ' ' ' (8Ci)2 " 8C4 ■ 

Under the additional assumptions 

(10.18) J < (5o and |A ~ g\ < -^ 

we have w^/^ < co/(4C4), so we conclude that \^'^Hg(j) — 2co| < cq, hence 

(10.19) """Hgcj) > Co > 0. 

This at once gives a positivity estimate for ^'^Hgq near ^o- Namely, 
^^iJgg - -Ao-i,5-ixo(^o"'(2 - '/'/'5))xi(y • ^^|5 + 2)-i/g0 

+ r 1x0(^0 ^ (2 - 0/<5))x'i (y • m/-^ + ^rHgT^. 
Thus, 

(10.21) '"'Hgq=-b^ + e 
with 

(10.22) 52 ^ Ao-irVo(^o '(2 - 0/'5))xi(y • m/'5 + 2rHgCJ,. 
Hence, with 

(10.23) b^ = coA^'S-^'aiAo'C^ - ^/S))xiiy ■ 1^/5 + 2), 
we have 

(10.24) """Ugq < -b^ + e. 
Moreover, 

(10.25) 52 > {caAo/l6)q 
since (t> > y ■ IJ- > —2(5 on supp q, so 

26) ^o(^"'(2 - HS)) = Al{2 - cl^/S)-\o{Ao\2 - <t>/5)) 

>(A2/16)xo(V(2-0/^)). 
On the other hand, e is supported where 

(10.27) -2(5 < 2/ • ^ < -(5, ^1/2 < 2/5-1/2^^ 

so, for (5 > sufficiently small, in the region which we know is disjoint from WFsc(u). 
Moreover, on supp q, 

(10.28) -2(5 < y • ^ < 2(5, ^1/2 < 2/3^1/2 ^, 

so, for (5 > sufficiently small, we deduce from the inductive hypothesis that 
supp (7 (hence supp 6) is disjoint from WFgj, (u). In addition, by choosing (5 > 
sufficiently small, we can assume that the support of q, e and b are all disjoint from 

WFsc((i?-A)w). 



58 ANDRAS VASY 

Moreover, with d denoting a partial derivative with respect to one of (y, z, r, /u, v), 

dq = -A^'S-'xoiAo\^ - HS))xi{v/S + 2)90 

~-S-'xo{A^\2~4>/S))x[{v/S + 2)dr,. 

As 2/ = is outside the support of the second term, and as 9^0 vanishes at y — 0, 
we conclude that for any multiindex /3, 

(10.30) \d^dq\ < Cpb^ at y = 0. 

More generally, at any Cb with p g Cb, defined by x = 0, y' = 0, as above, 4> is 
independent of fj! at y' = so outside supp e 

(10.31) \dl,dq\<Cph'^ aty' ^0. 

In fact, outside suppe, but in the set where h is positive, 

(10.32) b-^dq^Co^d(l), 



so the uniform bounds of ( 9.30 ) also follow. 

Let tp G C^(M) be identically 1 near and supported sufficiently close to so 
that the product decomposition of X near dX is valid on suppi/;. We also define 

(10.33) q^rljix)q. 



Thus, q e Qoo^scrp*^^ jg ^ TT-invariant function satisfying (9.4). Let A be the 
operator given by Lemma |9.l| with q in place of q, so in particular its indicial 
operators are g(C)'0o(^h,o(C))- Note that (|l|) holds with C = IGc^^A^^. So 
suppose that M > and e' > 0. Choose Aq so large that MC < e' . By Corollary 3.7 
and the hypothesis ( 10.1 ), we deduce the following statement. For any K' C ^'^T*X 
compact with suppe n iiT' = there exists S' > 0, B,E £ *g^°°'°(X,C), F e 
*g;^'^(X,C) with 

(10.34) 

WF^,(£;) n if' - 0, WF^,(F) C suppg, S,,o(C) = 6(C)<z(C)'/'^(i?a,o(C)), C e K', 

such that ii ip E C^(M) is supported in (A — 5', A + S') then 

(10.35) 

i'ilj{H)x-'^/'^[A*A, H\x~^/'^il){H) ~ M%l){H)A*Ail){H) > (2 - 2e')B*B + E + F. 

Let 

(10.36) Ar = x-^^^/^{l + r/x)-\ r e (0, 1), 

soA^ e *s^"'+^/^(A,C)forr e (0, 1) audit is uniformly bounded in *g^"'"^/^ (A, C). 
The last statement follows from {l + r/x)~^ being uniformly bounded as a 0th order 
symbol, i.e. from {xdx)''{l + r/x)~^ < Ck uniformly {Ck independent of r). We 
also define 

(10.37) Ar = AArX-'^/'^ij{H), B,. = BAr, Er = ArEAr. 
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Then, with ipo G C^(K; [0, 1]) identically 1 near suppV", 
(10.38) 

= i{l + r/xy'^'ilj{H)x~^/'^[A*A, H]x-^l'^ijj{H){l + r/x)"^ 

+ iilj{H)A*x^+^/^ [krX~^/^,H] (1 + r/x)-^x-^^^MH)AiP{H) 

+ iiP{H)A*MH)x-'^^{l + r/x)-' [ArX-^/^,H]x^+^^^Aip{H) + Hr, 

where Hr is uniformly bounded in ^S(°°' {X, C). Note that Hr arises by commuting 
A, powers of a; and A^ through other operators, but as the indicial operators of A 
and X are a multiple of the identity. A, x and A^ commute with these operators to 
top order, and in case of A^, the commutator is uniformly bounded as an operator 
of one lower order. Then, multiplying ( 10.35| ) by {l + r/x)~^ from the left and right 
and rearranging the terms we obtain the following estimate of bounded self-adjoint 
operators on L'^^{X): 

ix'+^/^[A;Ar,H]x'+^/^-tPiH)A*{G; + Gr)A^p{H)-M^{H)A*A'tPiH) 

^^^'^^^ > x'+i/2((2 - e')B:Br + Er + Fr)x'+'/^ 

where 

(10.40) Gr^tMH)x-'^^{l + r/x)-'[ArX-^/^,H]x^+^^^, 

and Fr G *g^°° "^'+^(X,C) is uniformly bounded in ^g^'^^\x,C) as r ^ 0. Now, 
Gr G ^Sc°°' (-^^ ^) ^^ uniformly bounded in ^g^ {X, C), hence as a bounded operator 
on Lg(,(X). Thus, if M > is chosen sufficiently large, then Gr + G* > —M for all 
r e (0,1), so 

(10.41) iP{H)A* {Gr + G;+ M)AxP{H) > 0. 
Adding this to ( 10.39| ) shows that 

(10.42) ix'+^^^[A;Ar, H]x'+^/^ > a;'+i/2((2 - e')S,!5, + Er + Fr)x'+^/^. 
The point of the commutator calculation is that in Fl^{X) 

{u\AlAr,H]u) 

(10.43) = {u, AlAr{H - \)u) - {u, {H - A)^;^!^^) 
= 2ilm{u,A*AriH - X)u); 



the pairing makes sense for r > since Ar G 'I'sc^' (^i^)- Now apply ( 10. 3£ ) to 
^-/-i/2y and pair it with x^'^^/^u in Ll^{X). Then for r > 

(10.44) \\Bru\\'^ < \{u,Eru)\ + \{u,Fru)\+2\{u,A;Ar{H ~X)u)\. 



Letting r — > now keeps the right hand side of (10.44) bounded. In fact, 
Ar{H - X)u e C°°(A:) remains bounded in C°°{X) as r -> 0. Similarly, by ( |10.34D , 
ErU remains bounded in C^{X) as r ^ if we chose K' so large that WFgc(u) C 
K'. A lso, Fr is bounded in B{H^^\X),H-^"^^-\X)), so {u,Fru) stays bounded by 
( 10.34| ) as well. These estimates show that BrU is uniformly bounded in Ll^[X). 
Since (1 + r/x)-'^ -^ Id strongly on B{H^''^' {X), H^' '^' (X)), we conclude that 
Bx^'^^^^u e L'^^{X). By ( 10.34 ) and Proposition |5.3| this implies that for every m, 

(10.45) ^0 ^ WF^:''+^/'(zi). 
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This is exactly the iterative step we wanted to prove. In the next step we decrease 
5 > sHghtly to ensure that WFg(,(i^) C suppq is disjoint from WFg"' (u). D 

To state and prove the finer estimate under the assumption that all elements of 
C are totally geodesic, first note that in geodesic normal coordinates around p e C", 
h— (|a*P + W\'^) vanishes together with its first derivatives at p = (0, 0). Thus, ^'^Hg 
agrees with W^ on ^'^T*X where 

(10.46) W^ = 2T{fi ■df. + i^-d^)- 2(|/ip + \iy\^)dr + 2iy -8^ + 2(1 ■ dy. 

We will use W^ to model the bicharacteristic flow of H. Note that W^ is the 
(rescaled) Hamilton vector field of the metric function r^ + |i^p + |/ip, i.e. where 
we replace the actual metric hhy a flat one. 

We remark that it is the 9^ and di, components of ^'^Hg that differ from W^ on 
^'^T*X if we do not assume that the elements of C are totally geodesic. The former 
is inconsequential since we only consider 7r-invariant functions (in particular, the 
only /x-dependence is via rj = j/-/x), but the latter rules out the more precise location 
of the singularities given in the following proposition. 



Proposition 10.3. Suppose that H is a many-body Hamiltonian satisfying (10. 1) 
and that every element of C is totally geodesic with respect to h. Let u G C~'°°{X), 
A > 0. Given K d E„(A)n"'=r(;^,(C;X) compact with Kr]WFsciiH-X)u) = there 
exist constants Cq > 0, Sq > such that the following holds. If £,q = (0, tq, vq) G K 
and for some < S < So, CqS^'^ < e < 1 and for all a = (j/, z, t, /i, i^) G '^'^Tq^X n 
Sa-a 
(10.47) 

a G "'T*,^X and \TT^^{exp{6W^)(a)) - ^o\ < eS and |2/(exp((5VK^)(a))| < eS 
=» TTQbia) i WFsc(w) 
i/ienCo^ WFs,(u). 

Proof. We again employ an iterative argument, so we assume that ^o ^ WFg^ (u) 
and we need to show that ^o ^ WFg^ (u). 

We first construct a C°° function uj oi z, t, v, rj ~ fi ■ y and s = \y\^ which 
measures the distance of bicharacteristics of A in Sa-a from tt^q (^o) H Sa-a- 
Thus, T^ + \v\'^ + l/ip — A will be small along these bicharacteristics. We will take 
uj of the form 

{y ■ i^r 

^-To - Wo\ 
where luq only depends on z, t, v and rj — y ■ fi. Note that 

(10.49) \y\'~^-^ = \y-'rii^\' 

is the squared distance of the integral curves of H\i_i\^, which are just straight lines, 
from y = 0, so near Sa-a the second term in to gives the fourth power of this 
distance. 

Pushing forward W^ by the map F : {y, z, t, h, v) i~^ (z, t,i^, h ■ y) at some point 
a — {y, z, r, /i, ly), we obtain the vector 

(10.50) F^l^W^ = 2(Tr/ + \pi\^)d,^ + 2tv ■ 8^ - 2(|^p + \v\'^)dr + 2v ■ 8,. 



(10.48) ^-^g + (!#- . Z,^\,, ? 
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|2_ 



Since we are interested in what happens near Sa-a H ^'^T*X , where A = r^ 
l/xp, we are led to consider the constant coefficient vector field 

(10.51) Wo = 2(A - r2 - \i.o\^)d„ + 2toUo ■ d. - 2(A - t^)^, + z/q • 9, 
in the variables (z, t, v, 77), so 

(10.52) F,\^W'' = 1^0 + 2(A - tI - |^o|' - \M^)){~d^ + 5r). 
Note that the 9^ component of Wq is nonzero. Let 

10.53 zo W = l£-i. ro{t) = TO + -^f, i^o W = '^0 + TT^t, 

WqT] WqT] WqT] 

so 

(10.54) j:t^{zo{{Wov)t),To{{Wov)t),MiWov)t)AWov)t) 

gives a curve through (^Qi 0) with tangent vector Wq- Now we define loq by 

(10.55) wo = (^ - zo{v)f + ir- Mv)? + (^ - ^o(^))' 

so ujQ vanishes exactly quadratically along 7 and is positive elsewhere, and 

(10.56) Wqujo = 0. 
Note that by the triangle inequality 

(10.57) \z\ + |t - rol + 1^. - z^ol + M < C^^^ + \v\) 
for sufficiently large C. 

Since for a G 7r~^(^o) we have F^,\a'^'^Hg — F^\aW^ — Wq, we see that 

(10.58) ^'Hgiz - zoin)) = at ^-^(^0), 
i.e. when y = 0, z = 0, t = tq, u = i/q, g = X, so 

(10.59) rHg{z - zo{v))\ < Ci\y\ + c.^/' + |A - g\). 
Hence, 

(10.60) rHg{z - zoivjfl < 2Cu;l^\\y\ + w^' + |A - g\). 
Similar conclusions hold for t — TQ{ri) and j^ — i'o(7?), so 

(10.61) rHgUol < C\{\y\ + Lo'/' + \X - .gDc.^^ 

Next, we calculate ^''Hg{\y\'^ — (y-//)^/(A — Tq — |i^oP))- Since the function we are 
differentiating vanishes quadratically at j/ = 0, the same follows for its derivatives 
with respect to any vector field tangent to y — Q. Since the dy component of ^"^Hg 
(and of W^) is of the form 2^- dy + Y^ Pjdyj with (3j vanishing at y = 0, z = (i.e. 
at p), we conclude that 

(10.62) \rHg - 2f, ■ dy){\y\' - {y ■ f^f/iX r^ l^oH)! < C,\y\i\y\ + \z\). 
On the other hand, 

(10.63) (2m • dy)i\y\' ~ {y ■ m)V(A - r^ ~ Wof)) = 4(m • J/) ^ ~ """ 'J''"', ~ '^'' ■ 

But, as in ( [lO.lOD , 

|A - r^ - Wo\' - ImH < |A - g\ + C'{\y\ + \z\ + \r - to| + W- i^ol) 

(10.64) ^,2 
<C^i\X-g\ + \y\+4^'). 
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Thus, 

(10.65) rH,i\y\^ - [y ■ ^^f/{\ r^ \u,\^))\ < C,\y\i\X - g\ + \y\ + 4^'). 
Our results thus far imply that 

(10.66) rHgiul < C5Lo'/^{\y\ + \X - g\ + uj^^f. 

Now let 1 > e > 0, (5 > 0, with other restrictions to be imposed on these later, 
and let 

(10.67) = ro - T + -Igco. 

We use tq — T to measure propagation along the bicharacteristics; rj = y-fj, would also 
work. We again let xo G C°°(R) be equal to on (— oo, 0] and Xo{t) = exp(— 1/t) for 
t > and we let xi e C°°(M) be on (-oo,0], 1 on [l,oo), with Xi > satisfying 
xi G C^((0, 1)). Furthermore, for Aq > large, to be determined, t e (0, 1), let 

(10.68) qt^q = XoiA-\l + t - 0/<5))xi((to - t + 5)/{e5) + t). 

We usually simply write q in place of qt . We only use t to slightly shrink the support 
of q in our inductiv e pro of (i.e. as I is increasing), instead of adjusting 6 as in the 



proof of Proposition 10.1. Thus, on suppq we have (j) < 26 and tq — r > —26. Since 



w > 0, the first of these inequalities implies that tq — t < 26, so 

(10.69) |t - Tol < 25 and lo < e^6^{26 + r - ro) < 4e^6^. 

Hence, luq < 2e'^6^, which together with \t — tqI < 26 gives |ry| — \ij. ■ y\ < Cq6 since 
the dr component of W in non-zero. Since we also have 

(10.70) \\y\' - {y ■ /i)V(A - r^ - l^oHl < 2e'6', 

we conclude that \y\ < CyJ. Thus, under the additional assumption 

(10.71) \X-g\<6 
we deduce that |'^'^-ffgw| < Cae'^6^, so 

(10.72) \'''Hg(j)-2h\<Cs6/e'^. 

Hence, for cq > 0, Co > appropriately chosen and for e E (0,1), 6 > satisfying 
6/e^ < Co, we have 

(10.73) """Hgct) > Co > 0. 

Again, this directly gives a positivity estimate for ^'^Hgq near ^o- Now 
(10.74) 

^^Hgq = -Ao i,5-ixo(^o '(1 + i - HS))xi{{ro -r + 6)/{e6) + tTHg^ 

- {e6)-\o{A,\l + t - HS))x'i{{ro - r + 6)/{e5) + tf^HgT. 
Hence, with 

b' = coA^'6-^x'oiAoH^ + t - 0/<5))xi((ro - r + 6)/{e6) + t), 
e = -{e6)-'xo{Ao\l+t-(f>/6))x[{{ro-T + 6)/ie6)+tY^HgT 
we have 
(10.76) '"Hgq < -6^ + g. 
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In addition, similarly to ( |l0.25D -( [l0.26 ), we see that 



(10.77) 6^ > (c„Ao/16)(7. 

Moreover, with d denoting a partial derivative with respect to one of (y, z, r, ^, ^), 

dq - -Ao-M-ixo(^o"'(l + t - HS))xi{{ro ' r + 6)/{eS) + t)d<t> 

- (6<5)- Vo(Ao ^(1 + i - 0/<5))x'i((^o - r + <5)/(6<5) + t)ar. 



Thus, ( |l0.30|) -( [l0.32|) hold, and hence the uniform bounds of ( |9.30 ) also follow. 
Now e is supported where 

(10.79) -5 - te5 < tq - T < -5 + {I ~ t)eS, uj^/^ < ^/2 eJ, 

so near the backward direction along bicharacteristics through i^Oi in the region 
which we know is disjoint from WFgc(u). In addition, by choosing S > sufhciently 
small, we can assume that the support of g, e and b are all disjoint from WFgc((-ff — 
X)u). 



From this point we can simply follow the proof of Proposition 10.1. Thus, we 
conclude that for every m, 

(10.80) Co i WF^/+^/'(i.). 



This is exactly the iterative step we wanted to prove. In the next step we decrease 
t slightly to ensure that suppgt is disjoint from WFg"' (u). D 

Before proving the general tangential propagation estimate, we first do it in the 



totally geodesic case (C totally geodesic). Proposition [7. l| shows that for sufficiently 
short times there is a unique generalized broken bicharacteristic through any point 
in S((A), namely the integral curve of ^^Ug. The simplicity of this description 
may already give a hint that it is particularly easy to prove the corresponding 
propagation estimate for singularities. Indeed, in the proof of the aforementioned 
proposition, we have essentially already constructed the pseudo-differential operator 
A to commute through H by defining the 7r-invariant function (which will play 
an analogous role to that of in the proof of normal propagation). The following 
argument may also clarify the close relationship between proving results about the 
geometry of the generalized broken bicharacteristics and the positive commutator 
proof of propagation estimates. Again, we only state it for forward propagation. 



Proposition 10.4. Suppose that H is a many-body Hamiltonian satisfying (10.1). 
Suppose also that every element of C is totally geodesic with respect to h. Let 
u e C-^{X), A > 0. Lei Co e St(A)n-T5,(C;X), C = Ca, satisfy i^ ^ WYsdiH- 
\)u). Then there exists e' > such that if in addition for some s G (— e', 0) we have 

noaieMs'"Hg){TT-\Co))) i WFsc(w) 
thcnini^^scin). 

Proof. First note that there is nothing to prove if ^o G R+W U i?_(A), so from 
now on we assume that ^o ^ ^+(A) U i?_(A). The proof is very similar to the 
previous one and the positive commutator construction is exactly the same as in 
three-body scattering [p9| Proposition 15.4], based on the 7r-invariant function cj) 



used here in the proof of Proposition 7.1. Thus, we take local coordinates centered 



at C as above, i.e. of the form (y, z), and let 4> = 4''"'^^ be defined by (7.21), so in 
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particular is 7r-invariant. In the proof of Proposition |7.1| we showed that there 
exists (5o G (0, 1) such that for any S £ (0, So) and any e G (0, 1) 

(10.81) <j>{i) < 26, T{i) -tq<2S and \T^{i) + h{£,) - X\ < eS 

imply that 

(10.82) 



'H„ 



satisfies ( |7.33D , so 

'"Hgm > Co > 0. 

We define q as in ( |10.68| ). Then ( |l0.74| ), hence ( |l0.75D -( p79l ) also hold. Since 
e > can be taken arbitrarily small, we can choose it and 6 G (0, So) so that 
suppe is a small neighborhood of exp(s*''^iJt,)(7r~-'^(^o)); in particular, 7ro6(suppe) is 
disjoint from WFsc(m) for each b. We can then apply the compactness argument 
of Proposition 10.1 to prove (10.39) for the operators A, B, etc., defined in that 
proof, and conclude that ^o ^ WFsc(u). □ 

We now return to the general setting of not necessarily totally geodesic C. 



Proposition 10.5. Suppose that H is a many-body Hamiltonian satisfying (10.1). 
Let u G C-°°{X), A > 0. Given 

(10.83) K C (St(A) n ''T^,iC;X)) \ (i?+(A) U i?_(A) U WFsdiH ~ X)u) 

compact there exist constants Cq > 0, Jp > such that the following holds. If 
Co = (-^Oj To, t'o) G K and for some Q < 5 < 5q, CqS < e < 1 and for all a = 

(2/,z,r,Ai,^)G^^TaV^nl]A-A 
(10.84) 

a G '"'T*.^X and |7rgJexp((5W^^)(a)) - ^q\ < eS and |y(exp(5W^^)(a))| < eS 

=> TTQb{a) i WFsc(m) 
f/ienfo^ WFs,(u). 

Proof. The proof is very similar to the previous ones and now the positive commu- 
tator construction follows that of [^ Proposition 15.2] in three-body scattering. 
Thus, we take local coordinates as above, i.e. of the form (y, z) with Cb defined by 
linear equations in y. Then we construct ujq G C°° {^'^T^r{C; X)) (defined near ^o) 
to measure the squared distance from integral curves of 

(10.85) W^ = 2ti^ ■ d^ - 2hdr + Hj^; 



this is achieved by solving a Cauchy problem as in [p9| and in (7.12) here. (In- 
deed, an approximate construction, like that of loq in the normal case discussed 
above, would also work). Then we extend wq to a function on ^"^Tq^X (using the 
coordinates {y, z,t, ii,v) near dX), let 

1 



(10.86) 



UJ =LJ0 + \y\ , (/) = To - T + 



;2,5 



^, 



and define q as in ( 10.13| ). The difference in the powers of e and S in this defi- 
nition of (/) in the (general) tangential setting and that in the normal case (given 
in (10.67)) arises since in the normal setting u approximates the fourth power of 
the distance from the generalized bicharacteristics while here it approximates the 
squared distance. The estimates on '^'^HgCJ) are just as in ||3^, Proposition 15.2], see 



also the proof of Proposition 7.1 here in the similar totally geodesic setting (the 
estimates are simply better but no different in nature under the totally geodesic 
assumption since now we do not have (|7.4D), giving a slightly better result than in 
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the totally geodesic normal case: it is (5/e, not 5J€^^ that has to be bounded below 
by an appropriate positive constant. The difference arises as the model integral 
curves in the tangential setting are closer to the actual ones than in the normal 



setting. Thus, one obtains (10.25) here as well. The functional analysis part, under 



the assumption that there are no bound states, is exactly as in the normal case. D 

An argument of Melrose-Sjostrand pa, E3|, see also [O, Chapter XXIV] and 
|l7[ allows us to conclude our main result concerning the singularities of general- 
ized eigenf unctions of H. Here we concentrate on totally geodesic C (since that 
is the case in Euclidean scattering), in which case the more delicate tangential 
propagation argument of Melrose-Sjostrand is not necessary. The proof presented 
below essentially follows Lebeau's paper |l^. Proposition VII. 1]. We thus have the 
following theorem. 

Theorem 10.6. Let (X^C) he a locally locally linearizable many-body space, and 



ssupposc that H is a many-body Hamiltonian satisfying (10.1). Let u € C^°°{X), 
A > 0. Then WFsc(w) \ WF sc((^ — ^)'^) ^^ '^ union of maximally extended gener- 
alized broken bicharacteristics o/ A — A. 

Proof. We only need to prove that for every a, if ^o G WFsc(w) \ WFsc((ff — A)w) 
and ^0 G ^'^Tq, {Ca;X) then there exists a generalized broken bicharacteristic 7 : 

[-60, eo] ^ S, eo > 0, with 7(0) = ^o and such that 7(f) e WFsc(m) \ WFsc((i/ - 
X)u) for t E [— Eq, Eq]- III fact, if this statement holds for all a with Cc C Ca, let 

(10.87) 
TZ —{generalized broken bicharacteristics 

7 : (a,/3) ^ (WFsc(^) \ WFsc((-ff - X)u)) D (Jl^^T^, (C^; X) : Cc C Ca}, 

a < < /3, 7(0) = Co}, 

and put the natural partial order on 7?., so 7 < 7' if the domains satisfy (a, /?) C 
(a',/3') and 7 = 7'|(q,/3)- Then TZ is not empty and every non-empty totally 
ordered subset of TZ has an upper bound, so an application of Zorn's lemma 
gives a maximal generalized broken bicharacteristic of A — A in the intersection 
of WFsc(u) \ WFsc((i/ - A)u) with Uc^ccJ^'T^'^iCalX) which passes through ^o- 
A similar maximal statement holds if we replace C^ C Ca by Cc ^ Ca- 
Indeed, it suffices to show that for any a, if 

(10.88) Co e WFsc(u) \ WFsciiH - X)u) and Co e ^^T^, (Ca;X) 

then 
(10.89) 

there exists a generalized broken bicharacteristic 7 : [— eo, 0] ^ S, eo > 0, 

7(0) = Co, lit) e WFsc(w) \ WFg,((7J - A)u), t e [-eo, 0], 

for the existence of a generalized broken bicharacteristic on [0, eo] can be demon- 
strated similarly by replacing the forward propagation estimates by backward ones. 



and, directly from Definition S.2, piecing together the two generalized broken 
bicharacteristics gives one defined on [— eo,eo]. 

Note that if every element of C is totally geodesic and Co S 5]f(A) then, due 



to Proposition 10.4, (10.88)^(10.89). Indeed, to prove this statement for Co G 
^'^T^, (Ca'jX), we only need that Cb be totally geodesic for Cb with Ca C Cb (since 
the result is local); in particular, it always holds at Cg. 



66 ANDRAS VASY 



We proceed to prove that (10.88) implies (10.89) by induction on a. As remarked 
above, this is certainly true for a = 0. We only prove the implication here under 
the assumption that all elements of C are totally geodesic, the general case repeats 
the argument of Melrose-Sjostrand. Thus, we already know the implication if ^o G 



So suppose that ( |10. 88 )^( 10. 89 ) has b een proved for all b with Ca CI Cb and that 



^0 e S„(A)n'''^T^, (Ca] X) satisfies ( |10.88D . As noted in the first paragraph, we thus 
know that the intersection of WFsc(u) \ WFsc((i? - A)u) with lic^cc\'"'T^' {Cb;X) 
is a union of maximally extended generalized broken bicharacteristics of A — A. We 



use the notation of the proof of Proposition 10.1 below. Let U he a, neighborhood 
of ^0 = (0, 2:0, To, I'o) in S which is given by equations of the form |?;| < S' , |z — zqI < 
S', \t - Tol < S', \v - vo\ < S', 6' > 0, such that """"Hgr) > on TJ-'^iU) and 
U n WFsc((i/ - X)u) = 0. Such a neighborhood exists since fo i WFsc{{H - X)u) 
and^'^Hgr]{£_o) = \—TQ—h{zo,i^o) > for every ^0 G ^^^{^o)- Also le t U' be a subset 



of U defined by replacing d' by a smaller S" > 0. By Proposition 10.1 , there is a 
sequence of points ^n S S such that ^„ £ WFgc(u), ^„ -^ fo as n — + 00, and ?7(^n) < 
for all n, so we may assume that ^„ G U' for all n. By the inductive hypothesis, 
there exist generalized broken bicharcteristics 7„ : [— eoiO] -^ S with 7„(0) — ^„ 
and such that if for aU t e [-e',0], < e' < eo, we have 7„(i) ^ ^''^T^, (Ca; X), 
then jn{t) e WFsc(u) for t G [—e',0]. But i] is increasing on generalized broken 
bicharacteristics in U since ^^HgT] > there, so we conclude that y{'yn{t))-n{^n (t)) = 
Vilnit)) < vilniO)) < for i € [-60,0], hence y(7n(i)) ^ 0, so 7„(t) G WFsc(w) 
for all t G [— eo,0]. By Proposition |6.11 , apphed with K = WFsc(m), there is a 



subsequence of 7„ converging uniformly to a generalized broken bicharacteristic 
7 : [-60,0] -^ WFsc(u). In particular, 7(0) = ^o and -/(t) G WFsc(u) for all 
t G [— eo, 0], providing the inductive step. 

Note that this argument for normal points ^0 does not use that the elements of C 
are totally geodesic; it works equally well in the general case. Thus, for non-totally 
geodesic C now we only need to consider ^0 G St(A), and, as mentioned above, this 
can be done by the Melrose-Sjostrand argument. D 

We remark that the result is optimal as can be seen by considering the Euclidean 
setting, taking potentials singular at a specified Ca, thereby placing ourselves into 
the three-body framework. As p^ shows, singularities do reflect in all permissible 
directions in general, the reflection being governed to top order by the (two-body) 
S-matrix of the subsystem. 

11. The resolvent 



Before we can turn Theorem 10.6 into a result on the wave front relation of the 
S-matrix, we need to analyze the resolvent. More precisely, we need to understand 
the boundary values 

(11.1) RiX±iO) = {H-{X±iO))-^ 

of the resolvent at the real axis in a microlocal sense. To do so, we also need 
estimates at the radial sets R± (A) . Since the Hamilton vector field of the metric g 
vanishes at i?+(A)Ui?_(A), the estimates must utilize the weights x^'^^ themselves. 
In this sense they are delicate, but on the other hand they only involve x and its 
sc- microlocal dual variable r, so they do not need to reflect the geometry of C. 
The best known positive commutator estimate is the Mourre estimate, originally 



MANY-BODY SCATTERING 67 

proved by Perry, Sigal and Simon in Euclidean many-body scattering [E6l , in which 
one takes q — x~^t with the notation of Section 0. Since it is easy to analyze 
the commutator of powers of x with H (in particular, they commute with V)^ the 
functional calculus allows one to obtain microlocal estimates from these, as was 
done by Gerard, Isozaki and Skibsted |^, Q]. Thus, nearly all the technical results 
in this section can be proved, for example, by using the Mourre estimate and 



Theorem 10.6. In particular, apart from the propagation statements, they are well- 
known in Euclidean many-body scattering. The generalization of these Euclidean 
results to our geometric setting is straightforward; the arguments essentially follow 
those in three-body scattering that were used in [^ . 

We first state the weak form of the limiting absorption principle, namely that for 
/ e C°°(X), R{\ ± ii)/, t > 0, has a limit in H^'^{X), m arbitrary, I < -1/2, as 
t ^ 0. To simplify the asymptotic expansions of R{X±iO)f which we also describe, 
for A > we introduce the functions 

(11.2) a±-a±.A-±— ^GC°°(X\Co.,i„g), 

and the set of polyhomogeneous functions A^y^ {X \ Casing) on X \ Co, sing with 
index set 

(11.3) /C = {{m,p) : m,p € N, p < 2m}. 

Recall from [Q that v £ ^5;^ {X \ Co, sing) means that v is C°° in the interior of X 
and it has a full asymptotic expansion at Cq which in local coordinates {x,y) take 
the form 

oo 

(11.4) t;(x,j/)~^^x^(logx)Xr(y), a,,reC^{C'o). 

j=0 r<2j 

Thus, V € C"(X \ Co,sing) and \v{x,y) - ao,o(2/)| < Cx\ logxp. 



Theorem 11.1. Suppose that H is a many-body Hamiltonian satisfying (10.1), 
A > 0. Let f £ C°^'{X), uf = R{\ ± it)f, t > 0. Then uf has a limit u± = 
R(\ ± iO)/ m iJ™''(X), I < -1/2, as t^O. In addition, 

(11.5) WFsc{u±)ciR^iX). 
IfV is short-range, i.e. V £ x^C°°{X \ Co, sing), then 

(11.6) u± - e±'^/-x("-i)/\±, v± e C°^{X \ Co,™.g), 
while if V is long-range, i.e. V merely satisfies ( p.2| ), then 

(11.7) u± = e±*VA/.^(„-i)/2+»±^^^ ^^ g j^K^^(^x \ Co,«„g). 

Remark 11.2. The first statement in the theorem also holds if we merely assume 
/ e i/™'''(X) with I' > 1/2, but then WFsc(w±) has to be replaced by the filtered 
wave front set WFg|l' ^ (u±). Moreover, R{X ± iO) give continuous operators from 

H^-''{X)toH^+^''{X). 

Proof. This result is a weak form of the limiting absorption principle and can be 
proved by a Mourre-type estimate. In the Euclidean setting, it is a combination of 
the Mourre estimate, proved by Perry, Sigal and Simon |2q], and its microlocalized 
version obtained by Gerard, Isozaki and Skibsted Q. In the geometric setting, 
the Mourre estimate describes the commutator of H with a self-adjoint first order 
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differential operator A e x^^ r)iSl^{X) such that A - xD^ <E r)iSl^{X) (this is 
of course a restriction only at dX). Namely, it says that for cj) G C^(M; [0,1]) 
supported sufficiently close to A, we have 

(11.8) i(/)(i/)[A,i/]0(i/) >2(A-e)(/)(i/)2+i?, e>0, 

where R G ^sc°°' (^i^), hence compact on Ll^{X). It was proved in the geometric 
three-body setting (with an appropriate adjustment to allow bound states of subsys- 
tems) in [|9| , following the Euclidean argument of Froese and Herbst |^ . The proof 
given there goes through essentially unchanged for more than three bodies. Under 



our assumption (10.1), the symbolic commutator calculation in the scattering cal- 
culus, '''=iJg(a;-V) + 2g e xC°° {''''T* X) , and a shght modification of Corollary |7, 



prove the Mourre estimate. The argument of ||26| then proves the existence of the 
limits u± in H^i^{X), I < -1/2, and {H - X)u± = / € C°°{X) shows that the same 
holds in _ff™''(X) for every m and for every I < —1/2. 

To show the fiavor of the arguments, we prove here a version of the estimate of 
Gerard, Isozaki and Skibsted l^] . Such arguments as this can be combined to prove 
the limiting absorption principle without a direct use of the Mourre estimate as was 



done in the geometric two-body type setting by Melrose |21 and in the geometric 
three-body setting in [p9| . Here, however, we concentrate on proving the wave front 
set result. The major difference between the propagation estimates of the previous 
section and the ones near R±{X) is that ^"^Hg is radial at i?+(A) U i?_(A): it has 
the form 2Txdx- Thus, we need to use a weight x~''~^ to obtain a positive symbol 
estimate. So for ^ > —1, let 

(11.9) q = x-^-^xi'r)Tl'{x)>0 

where ■0 G C^ (M) is identically 1 near and is supported in a bigger neighborhood 
of (it is simply a cutoff near dX), x ^ C^(K; [0, 1]) vanishes on (-co, V^ — 2e), 
identically 1 on (\/A — e, oo), e > 0, x' > 0, and x vanishes with all derivatives at 
every t with x(i) = 0. Then for sufficiently small S > 0, \g — X\ = \t^ + h — X\ < S 
implies 

"""H^q = -2((/ + l)rx(r) + h^ {t))x-'-' < ~\?x-'-\ 

& = (2(? + l)Tx(r) + (A - rV(r)/2)i/2. 

Thus, both x^^^q and x'+^6 are 7r-invariant. Let A £ ^sc°° (^i^) be a quan- 



tization of q as in Lemma 9.1, except that now q is not supported in a single 
coordinate chart, so we need to define A as the sum of localized operators (of 
course, this is not necessary in the actual Euclidean setting). Thus, roughly speak- 
ing, A is the product of a quantization of q and ipo{H), ■i/'o € C^(IR). The fact that 
q £ x^^^^C°° {^^T*X) does not cause any trouble, and the argument of Corollary 3.7 
shows that for tp E C^(R; [0, 1]) supported sufficiently close to A we have 

(11.11) 

ix^+^/^iJj{H)[A*A, H]i;{H)x^+^/^ > x^+^^^{{2 - e')B* B + F)x^+^''^ , e' > 0, 

where 

(11.12) F e M'gf ■-2'(X,C), WF's,(F) C supp(x'+ig), 

(11.13) B e vf-f -'-i/'(x,C), B,,_;_i/2(C) = KC)9(C)'/V(ffa,o(C))- 
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H-{X + it))u+)-2t\\Anu+\\'^, 

(11.16) \\Bu+\\'+2t\\A„u+f<\{ut,Fut)\+2\{ut,A*,AoiH-i\ + tt))ut)\. 

Since i > 0, the second term on the left hand side can be dropped. Since u^ —^ u^ 
in H°/{X) for /' < -1/2, we conclude that for / e (-1,-1/2) the right hand 
side stays bounded as t ^ 0. Thus, Bu^ is uniformly bounded in L'^^{X), and 



Let 




(11.14) Ao = 


--A^iH)e^^:^^ 


Since 




(11.15) {u+,i[AlAo,H]u+) 


= -2Im(w+,A;.4 


we conclude that 





as u+ -^ u+ in ij£;' (X), we conclude that Bu+ G L'^d^)- But then (11.13) 



shows that for any ^ with q{() ^ 0, we have C ^ WF^' ' (u_|_) for every m 
This proves that the incoming radial set, i?+(A), is disjoint from WF^' ' (w+). 



1 + 1/2 G (—1/2, 0). Iterating the argument, as in the proof of Proposition 10.1, gives 
that WFsc('u+) n i?+(A) = 0. Since WFgc(u+) is closed, the same conclusion holds 
for a neighborhood of i?+(A). Finally, as all generalized broken bicharacteristics of 
A — A tend to R+{X) as t ^ — oo a nd ( H — A)u+ = / € C°°{X), the propagation 



of singularities theorem. Theorem 10.6 , implies that WFgc(M+) C i?_(A). The 
existence of the asymptotic expansions is a local question, so at Cq we can work in 
the scattering calculus to prove it, see Q] for details of the proof. D 

A pairing argument immediately shows R{X ± iO)v also exists for distributions 
V g C^°°{X) with wave front set disjoint from the incoming and outgoing radial set 



respectively. Combining it with the propagation theorem. Theorem 10.6, we can 



deduce the following result; as usual, we assume that {X,C) is locally linearizable. 



Theorem 11.3. Suppose that H is a many-body Hamiltonian satisfying ( 10.1 ), 



A > 0. Suppose also that v G C"°°(X) and WFsc(w) H i?+(A) = 0. Let uf - 
R{\ + it)v, t > {). Then uj' has a limit u_|_ — R{X + iO)v in C^°°{X) as i — > and 
WF5c('U-|-) ni?+(A) = 0. Moreover, if ^ € I]\i?_(A) and every maximally backward 
extended generalized broken bicharacteristic, 7 : (—00, ig] ~* S, with 7(^0) ^ ^ is 
disjoint from WFsc{v), then ^ ^ WFsdu^)- The result also holds with i?+(A) 
and R-(X) interchanged, R{X + it) replaced by R{X~it), (— oo,io] by [to, 00) and 
correspondingly 'backward extended' by forward extended'. 

Proof. As mentioned above, the first part follows from the self-adjointness of H, so 
that ioit>0,ve C-°°{X), f e C°°{X), we have v{R{X + it)f) = R{X + it)v{f)- 
recall that the distributional pairing is the real pairing, not the complex (i.e. i^) 



one. The wave front statement of Theorem 11.1 and the assumption on v show the 
existence of the limit U+ = _R(A + iO)w in C~°°{X) and that in addition WF^' (u +)C\ 



R+{X) = for every I < —1/2. The positive commutator argument of Theorem 11.1 
then applies and shows that WFgc(M+) n-R+(A) = 0. In the Euclidean setting these 
results follow from a microlocalized version of the Mourre estimate due to Gerard, 
Isozaki and Skibsted |0]; see ^ for a detailed argument. 

Finally, since WFsc(u+) is closed, a neighborhood of i?+(A) in S is disjoint from 
WFgc(u+). Since all generalized broken bicharacteristics approach R+{X) as t 



-00 by Propositions^, the last part follows from {H—X)u+ = v and Theorem 10.6. 



It can be also proved by modifying the argument of Propositions 10.1-10.5 along the 
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lines of our proof of Theorem 11.1. Namely, we consider the family u^ ^ C °^{X), 
t > 0, and note that for t > 0, R{X + it) e ^~^'°{X,C), so WFgcK^C WFsc(w). 
Let Aq, etc., be defined as Ar with r ~ where Ar is given by ( |10.37D (i.e. we do 
not need to use the approximating factor (1 + r/x)^^). Then 

(11.17) {u+,t[A*oAo, H]ut) = -2Im(u+, A^A^iH - (A + it))u+) - 2t\\Aoutf. 

Note that the pairings make sense since now WFgj.(Ao) is disjoint from WFsc(^tf^), 
t > 0. Thus, 

(1118) ll^o"^ll' + 2^ll^""^ll'^l("^'^o^^)l + l("^'^""^)l 

+ 2\{u+,A*„AoiH-iX + it))u+)\. 

Since t > 0, the second term can be dropped from the left hand side. Thus, knowing 
that uf — > M+ in C"°°(X) as i ^ 0, and assuming that ^q ^ WF^' (m_|_) has already 
been proved and ( 10.47 ) is satisfied by u+, we co nclude that ^o ^ WF^' ' (m+) 



The iteration of this argument of Proposition 10.1 and the similar arguments for tan- 



gential propagation allow us to conclude the forward propagation estimates which 



can then be turned into maximal statements as we did in Theorem 10.6. This 
argument also shows the influence of the sign of t: if i < 0, the inequality ( 11.18 ) 
cannot be used to derive results on m+. Instead, the signs are then correct in the 
backward estimate, just as expected. D 

We conclude this section with the following uniqueness theorem on solutions of 
{H — X)u = 0. It is essentially a geometric version of Isozaki's uniqueness theorem 
|l5| , Theorem 1.3], though we allow arbitrary growth of u away from one of the 
radial sets, say i?+(A). 



Theorem 11.4. Suppose that H is a many-body Hamiltonian satisfying (10.1), 
A > 0. Suppose also that u € C-°°(X), [H - A)u = and WF"j'(u) n R+{X) = 
for some m and some I > —1/2. Then u = 0. The same conclusion holds if we 
replace i?-|_(A) by R^{X). 



Proof. Just as in the proof of Theorem 11. 3|, the positive commutator estimate 



of Theorem |ll.l| (but now applied with a regularizing factor in a;) shows that 



WFsc(u) n i?+(A) = 0, and then Theorem ^Oj shows that 



(11.19) WFsciu) d R-{X). 

We remark that although we need a regularizing factor here which requires some 
changes in the proof, e.g. see the argument of the paragraph below, the regularizing 
factor (whether (1 + r/x)~^ or another one) commutes with V , so the additional 
arguments for dealing with it are essentially the same as the two-body ones. Thus, 
the regularization part of the proof of WFsc(u) H R+{X) = essentially follows [ pH 
Proposition 10]. 

We proceed to show that 

(11.20) meM., /< -1/2^ WFg"''(u)ni?_(A) =0. 

We give the details below since regularity arguments for distributions which are 
large at infinity seem to appear less often in the literature than the 'finer ones'; in 
particular, |jl3. Theorem 1.3] assumes u G if™''(X) with ^ > — 1. We essentially 
follow the proof of ]^, Proposition 9] below. 
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So suppose that (11.20) has been shown for some I < —1; we now show it with 
I replaced by / + 1/2. This time we consider 

(11.21) q = x-'-\{T)tl;{x), K-1, 

where ip £ C^ (K) is identicaUy 1 near and is supported in a bigger neighborhood 
of (it is simply a cutoff near dX), x & C^(R; [0, 1]) identically 1 on (— oo, — -s/A+e), 
vanishes on (— \/A + 2e,oo), e > 0, and x vanishes with all derivatives at every t 
with x(i) = 0. Then 

(11.22) ^^H,q = -2{{l + l)rx(T) + hx'{T))x-'-' = {-b' + e)x-'-\ 

(11.23) 62^2(/ + l)Tx(r). 

The first key point now is that on WFsc(w) we have r — — vA, so WFsc(w) n 
7r(suppe) — 0. Let A S ^sc°°'~^ (^iC) as in Lemma 9.1. Corollary 9.7 again 
shows that for ip £ C^(R; [0, 1]) supported sufficiently close to A we have 

(11.24) 

ix^+^/^ijiH)[A*A, H]iP{H)x^+^/^ > x'+i/2((2 - e')B*B + E + F)x^+^''\ e' > 0, 

where 

(11.25) E£^s:°^-''-'{X,C), WF's,(ii;)nWFsc(ii) =0, 

F £ ^-:^'-'\x,C), WF's^iF) C supp(x'+ig). 

Let 

(11.26) 

Ar^A{l + r/x)-^^P{H), Br = B{l + r/x)-\ Er = {1 + r/x)-^E{l + r/x)-\ 

so 
(11.27) 

Ar £ ^s^-^\x,C) for r > 0, Ar is uniformly bounded in *s^°^'"'"^(X,C); 
analogous statements also hold for B^ and Ej.. Thus, 

tx'+^^^[A;Ar,H]x'+^/^ 

(11.28) = 1(1 + r/x)-^x^+'^/'^^{H)[A*A, H]iIj(H)x^+^/^{1 + r/x)-^ 

+ iij{H)A*x^+^{Gr + G;.)x^+^ Aij{H) + Hr 
where Hr is uniformly bounded in ^sc°° (^i^) ^'^'^ 

(11.29) Gr - iMHf^~\^ + ^/2;)"M(1 + r/a;)-\ iJ], 

tpo £ C^iM,; [0, 1]), ^0 = 1 on suppf/", so Gr is uniformly bounded in ^g^f"' {X,C). 
Thus, we need to estimate the commutator [(1 + r lx)~^ ,H], and now we do not 
have a large M as in the proof of Proposition 10.1 to help us deal with it. 

The other key point is thus that we have z[(l + r/x)^^, H] — i[{\ + r/x)"^^ A] 
and 
(11.30) 

-H,(l + r/x)-' ^ 2r-^ = -cl + d,„ c, = Xi(r)(-r)i/2 (^\ , 
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Xi G C^(IR;[0, 1]) identically 1 on (— oo, — \/A + 3e), vanishes on (— a/A + 4e, oo), 
e > 0. Let Cr be the quantization of Cr multiplied by ipo{H) as in Lemma p3, and 
define Dr similarly but with V'o(^) replaced by ipo{H)'^. Thus, as (1 + r/x)~^ is 
uniformly bounded in the symbol class S'°(X), 

(11.31) iMH)x-^^^[{l + r/x)-\H]x-^/^MH) = C*,Cr + Dr + K 

with Cr and Dr uniformly bounded in ^~^'°{X,C), Cr & %^'^^'^iX,C) for r > 0, 
Dr e *g^°°'^(X,C) for r > 0, and ij; uniformly bounded in *7^°°'^(X,C). Moreover, 
-Dr^ G ^g^°°'°°(X, C) uniformly due to the disjoint operator wave front sets. Thus, 

(11.32) Cr + C; = 2(1 + r/x)-^/^{C*rCr + Dr){l + r/x)-^/^ + H'; 
with H'r uniformly bounded in '^q^' {X,C), so 

(11.33) 
il;{H)A*x^+^{Cr + C*r)x^+^A^H) 

= 2i;{H)A*x^+^{l + r/x)-^'^{C*rCr + Dr){l + r/x)-^'^x^+^Ai;{H) + H^ 

, — OQ.l 



Hi, Hi uniformly bounded in *s^°°'^(X, C). Combining ( |ll.24| ), ( |ll.2S| ) and ( |ll.33| ), 
we see that for e' > we have 



(11.34) ix^+^'^[A*rAr, H]x^+^/^ > a;'+i/2((2 - e')B*rBr + Er + Fr)x^+^ 



/2 



We deduce as at the end of the proof of Proposition 10.1 that WFg*' ' {u) n 
i?_(A) = for every m and for every 1 + 1/2 < —1/2, so ( |ll.2C ) holds. In particular, 
u G H^'^X) for every m and for every / < —1/2. 

In the Euclidean setting we can now simply refer to Isozaki's uniqueness theorem 
pSf Theorem 1.3] to conclude that u = 0. Here we give some details to indicate 
how this conclusion can be reached in general. The crucial step is improving the 

=^ — 1/2 

estimate past the critical regularity Hsi (X). In the Euclidean setting this was 
done by Isozaki Q, Lemma 4.5] and his argument was adapted to the geometric 
setting in |3|, Proposition 17.8]. We thus conclude that WF™/('u) n i?_(A) = for 



Sc 

I e (0,-1/2). This is the point where (7J-A)u = 0, and not just {H - X)u € C°° {X) 



is used. Now we can apply a commutator estimate like that of Theorem |1 1 . 1| but 
near i?_(A). Thus, we conclude that WFsc(u) n i?_(A) = 0, so w € C°°{X). The 
theorem of Froese and Herbst Q on the absence of bound states with positive energy 
adapted to the geometric setting, as discussed in [^9|, Appendix B], concludes that 
u = 0. D 

12. The Poisson operator and the scattering matrix 

Just as in |35|, |39|] where three-body scattering was analyzed, the propagation 
of singularities for generalized eigenfunctions of H implies the corresponding result 
for the (free-to-free) scattering matrix, S{X), of H. Note that this is the only S- 
matrix under our assumption of the absence of bound states of the subsystems. We 
start by discussing the Poisson operator, then we use it to relate the propagation 
of singularities for generalized eigenfunctions to the wave front relation of the S- 
matrix. 
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The result that allows us to define the Poisson operator is that if V is short- 
range, i.e. V € x^C°^{X \ Co, sing), then for A g (0,oo) and g G C^{C'q), there is a 
unique u G C~°°(X) such that {H — X)u — 0, and u has the form 

(12.1) u = e-'^^a;("-i)/2^_ + R{\ + iO)f, 

where f_ G C°°{X), V-\ox = 9, and / G C°°{X). For long-range V the same 
statement is valid with the asymptotic expansion replaced by one similar to that 
of Theorem 



11.1 



(12.2) u ^ e-*^"x("-i)/2+'"-f_ + R{X + iO)f, v^ G ylphg(X). 
The Poisson operator with initial state in the free-cluster is then the map 

(12.3) P+(A):Cr(C^)^C— (§';), P+{X)g^u. 

(Note that the subscript for the free cluster has been dropped here in contrast to 
the introduction and ||38| .) To see that such a u is unique, note that the difference 
V = u—u' of two distributions u and u' with the above properties satisfies {H—X)v — 



and WFg;r(w) n R+{X) = by Theorem |ll.l| , so v = due to Theorem |llj. To 



see the existence of such u, note that as suppg C C'q, we can construct 

(12.4) 

u_ = e-*^"a;("-i)/2^_, w_ G C°^{X), v^\gx = g, -f = {H - A)w_ G C°^{X), 

by a local calculation as in |gl[, i.e. essentially in a two-body type setting. (We need 
to make slight changes in the asymptotic expansion for long-range V as described 
above.) Thus, we construct the Taylor series of u_ at dX explicitly, so we can even 



arrange that supp?;_ n Co,sing = 0- Then m = m_ + R{X + iO)f is of the form (12.1) 
and satisfies {H — X)u = indeed. 

We need to understand the Poisson operator better before we can extend it to 



distributions. So first recall from |39, Section 19] that the Melrose- Zworski [ p4| 
construction of a parametrix for the Poisson operator in the two-body type setting 
(C is empty) gives 'the initial part' of a parametrix P+{X) for the Poisson operator 
with free initial state in the many-body setting (for three bodies in that paper, but 
this makes no difference). Although the construction is performed there for short 
range potentials, it can be easily adjusted to long range potentials decaying like x, 
see ^, Appendix A] and [^ Section 3]. In particular, the kernel of P±{X) is of 
the form 

(12.5) K'^ix,y,y') = eT^^™^'i-*(^'^')/^x^"^(«')a±(a:,y,2/')M/i|, 

where dist is the distance function of the boundary metric h, \dh\ is the Riemannian 



density associated with it, a± are given by (1.20), and a± G C°°{X x Cq) are cut 
off to be supported near y = y' . Here y' is the 'initial point' of the plane waves, so 
y' G Cq corresponds to considering free incoming particles. In Euclidean scattering 
Kj. takes the form e^™'^ ci±{w, y')\dh\, w = y/x is the Euclidean variable and \dh\ 
the standard measure on the sphere; and e.g. if the potentials Vb are Schwartz then 
a± are just cutoff functions supported near y = y' which are constant in a smaller 
neighborhood of y = y' . In general, a±(0, y, y) is determined by the condition that 

(12.6) P±Wg = e=F»^/^2;"T+("-l)/\^^ 

v± G Apy^g{X), v±\dx ~ 9, and then a±{0,y,y'), as well as the other terms of 
the Taylor series of a± at x = can be calculated from transport equations near 
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y = y' . Finally, we cut off the solutions to the transport equations close io y — y' 
before reaching Co, sing- Thus, for y' in a fixed compact subset in Cg, K{x, y, y') is 
supported away from Casing, so for g e C~°°{C'q), supp(P±(A)g) is disjoint from 

^0,sing- 

The most important properties of P± (A) are summarized in the following propo- 
sition. Although we state them for P+(A) only, they also hold for P_(A) with 
the appropriate sign changes. Here we use ~j_ as the relation on S*dX x Sa-a 
defined analogously to ~± (see Definition 3.£ ), but with 'generalized broken bichar- 
acteristics' replaced by 'bicharacteristics of A — A'. Note that generalized broken 
bicharacteristics are simply bicharacteristics in ^'^T^, X which is where we will apply 
to following result. 



Proposition 12.1. (^, Proposition A.ljj Kl e C-°°{X x Co;f}/j), constructed 
above, is the kernel of an operator P+{X) '■ C^{Cq) -^ C^°^'{X), which extends to 
an operator P+(A) : C-°°{C'o) -^ C-°°{X), and for g G C-°°{C'o) 



(12.7) 



(12. 



SUpp(P+(A).g) n Co, sing = 0, 

WF,,(P+(A)(?)c{(j/,VA,0): yesupp.g} 

U le e S \ (i?+(A) U i?_(A)) : 3C G WF(5), ^ 



CI, 



(12.9) 



WFsciiH - X)P+{\)g) 

C {C e S \ (i?+(A) U i?_(A)) : 3C e WF(g), ^ -'- C}- 
The actual Poisson operator is then given by 

(12.10) P+(A) = P+(A) - R{X + iO){H - A)F+(A), 
with a similar definition of P_(A): 

(12.11) P_(A) = P_(A)-i?(A-iO)(ff-A)P_(A), 

Indeed, if g e C^{C^) then {H-X)P+{X)g e C°°{X) and P+{X)g has an asymptotic 
expansion as in ( 12. 6D , so by Theorem 11.1, {H — X)P^{X)g = and P{X)g has the 
form ( 12. 1| ) (with changes as indicated in (12.2) if V is long-range). In addition, for 
g e C^°°{C'o), WFsc((g - A)F± ( A)ff) is disjoint from i?±(A) by Proposition 1121' . 
Hence, by Theorem 11.3 , ( 12.10| )-( 12.11 ) indeed make sense. We also immediately 
deduce from Theorem ll.S 



Proposition 12.2. Suppose that H is a many-body Hamiltonian satisfying (10.1). 
Then the Poisson operator P+{X) : C^(Cq) —^ C~°°{X) extends by continuity to an 
operator P+(A) : C-°°{Cl^) -^ C-°°{X). Moreover, for g e C-°°{Cl^) we have 

WFsc(P+(A)g)c{(2/,VA,0) : y G supp.9} U i?-(A) 

U le e t{X) \ R+{X) : 3C e WF(5), ^ -- C}- 

Our definition of the free-to-free S-matrix is based on asymptotic expansions of 



(12.12) 



generalized eigenf unctions. So let g G C^{C'q) and let u = P+{X)g. By (12.1) 
(modified as in (12.2) for long-range V) and Theorem 11.1, u has the form 



(12.13) 



-»Vax^(„-i)/2^_ + e'^^x'^^-''^/\+ 
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S-matrix by 



g. We then define the free-to- free 



(12.14) 



5(A) : C^{C',) ^ C°°(C^), SiX)g = v+\ci. 



We need a better description of the S-matrix to describe its structure. This can 
be done via a boundary pairing formula analogous to pl|, P roposition 13]. It gives 
the following alternative description of tlic S-matrix, seepq, Equation (5.7)] (or its 
analogue from [p5| in the non-Euclidean setting): 

Proposition 12.3. For A > the scattering matrix is given by 

1 



(12.15) 



S{\) 



2iVX 



((i/-A)P_(A))*P+(A). 



Proof. The following pairing formula was proved by Melrose |2l|, Proposition 13] 
for short-range V, but the same proof also applies when V is long-range. Also, the 
proof can be easily localized, see pq, Proposition 3.3]. Suppose that for j — 1,2, 



Uj e C 



(12.16) 

and fj - 
(12.17) 



Vj,± e ^phg(^ \ Co, sing), SUpp(w2,±) <^ X\ Co^sing, 

{H - X)uj e C°^{X). Let Oj- ± = Wj,±|ax- Then 

2i^/X {ai,+ a^- ai-a^)dh= / {ui f2 - fiu^) dg. 
Jdx Jx 



We apply this result with ui — P{X)g, U2 — P{—X)f. By the construction of P(— A) 
we conclude that a2,+ = /, 02,- — 0, while for ui we see directly from the definition 
of S{\) and P(A) that ai,_ = g, ai_+ = S{\)g. Substitution into ( 12.17 ) proves the 
proposition. D 



Propositions 12.1 and 12.2, when combined with (12.15), allow us to deduce the 
structure of the S-matrix. 

Theorem 12.4. Let (X^C) be a locally linearizable many-body space. Suppose that 
H is a many-body Hamiltonian satisfying (10.1). Then the scattering matrix, 5(A), 
extends to a continuous linear map C~°°{Cq) —f C~°°{C'q). The wave front relation 
of 5(A) is given by the generalized broken geodesic flow at time n. 

Proof. Let /, 5 G C~°°{C'q). Suppose also that there is no generalized broken ge- 
odesic of length TT starting at some C G WF{g) and ending at C' € WF(/). That 
means that for any ^ e S \ (i?+(A) U i?_(A)) we cannot have ^ ^-. C: C G WF(5), 
and ^ ^+ C\ C G WF(/), at the same time. Proposition 12.1 (with — signs instead 
of -|-) implies that 

(12.18) WFsc((i? - A)P_(A)/) C t \ (P+(A) U P_(A)); 

indeed, we also have WFsc((iif — A)P_ (A)/) C ^'^T^,X, so we can even replace WFsc 



byWFsc. Thus, by our assumption on WF(/) and WF(g), and by Propositions |12T 
12.2 , we have 

(12.19) WF,e((i/-A)P_(A)/)nWF,e(P+(A).g) =0. 
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But the complex pairing 

(12.20) {u,u')x = I ulUdg 

extends by continuity from u,u' G C^{X) to u,u' e C~°°{X) satisfying WFsc(m) H 
WFsc(u') = 0. To see this just let A e *°e°(^) with WF^^(A) n WFsc(w) = 0, 
WF^^(Id-yl*) n WFsc(m') = 0, and note that 

(12.21) {u,u')x = {Au,u')x + {u, {Id- A*)u')x 
extends as claimed. Hence, the pairing 

(12.22) (P+(A).9,(if-A)P_(A)/)x = (((i/ - A)P_(A))*P+(A)g, /)x 

defined first for f,g £ C^{Cq) extends by continuity to f,g £ C^°°{Cq) satisfying 
our wave front condition. In other words, g can be paired with every distribution 
whose wave front set has no elements related to WF(g) by the generalized broken 
geodesic flow at time tt. Thus, for any A € ^^(Co) with WF'(^) disjoint from 
the image of WF{g) under the generalized broken geodesic flow at time tt, and 
for any / e C-°°(C^), {AS{X)gJ)9x = {S{X)g,A* f)9x is defined by continuity 
from / € C^iCli), so AS{X)g £ C°°{Cl)). But this states exactly that WF(S'(A)5) 
is contained in the image of WF{g) under the generalized broken geodesic flow at 
time TT. n 



Appendix A. The proof of Proposition 3.3 



In this appendix we prove Proposition 3.3 under the assumption that C is totally 



geodesic, roughly following Lebeau's original proof in ||lj. As noted after the 
statement of the proposition we can proceed inductively, using the order on C. So 
assume that 7(^0) — Co £ S„(A) fl ^"^T^,, {Ca; X). The inductive hypothesis is that 



we have already proved the proposition for b with Ca C Ct,- Thus, by Deflnition 3.2, 
part (ii), there exists 6' > such that the conclusion of the proposition holds if we 
replace io replaced by i 7^ to, assuming \t — io| < <^'- Let C±{t) £ Sa-a, t ^ to, be 
the points given by the inductive hypothesis. We often write 

(A.l) ^±{t) = iy{t),z{t),T{t),^i±iy),,y{t)) 

in local coordinates, so e.g. T(^±(t)) — T{t). Note that 7r(^±(i)) — 7(i), hence the 
independence of the 7r-invariant coordinates, y, z, r and v, of the ± signs. 
Notice first that r is tt- invariant, so for t ^ to we have 

(A.2) 

d{Toj)/dt\t±^''HgT{i±{t)) = -2/i(e±(t)) - 2(r(e±(i))2 - A) = 2(T(7(i))2 - A) 

where we used that t^ + h = X in Sa-a- Thus, T{t) = T{'y{t)) is differentiable on 
(to — S' ,to + S') except possibly at to, it is continuous at to, and its derivative T'(t) 
extends to a continuous function on (to — 5' ,to + S'). Hence T(t) is differentiable at 
to and r'(to) = 2(T(to)^ — ^) = ^'^HgT{Co) for any ^0 £ ^A-\- Notice also that, with 
the notation of ( 10.5l| ) in the proof of Proposition 10.1 , T'(fo) — Wqt — {^^HgT){io)- 



In particular, 

(A.3) |r(t) -To\<Ci\t- to\ if \t - to\ < S'. 



MANY-BODY SCATTERING 



77 



In fact, the ODE T'{t) = 2{T{t)^ - A), satisfied for \t - to| < (J', has a unique C°° 
solution, so on (to — (5', to + 5'), T(i) is C°° and 

(A.4) \T(t) - {to + {W„T){t - io))| < C\t ~ iol'. 

From now on we only consider differentiability issues from the left at to; of course, 
the situation on the right is similar. We define the vr-invariant functions j] = y ■ jj,, 
ujo, oj and 4> — (j)^'^'^^ as in the proof of Proposition 10.1. It is shown there that 



there exist Co > and So > such that if e € (0, 1), (5 € (0,(5o), S < Cqc^ and 
^ = (y, z, r, n, v) e Sa-a satisfies To — r > —25 and 0(^) < 2(5 then ^^Hgcj) > cq > 0. 
So suppose that we fixed some 

(A.5) 0<T<min((5',Ci(5o) 

and let 

(A.6) 5 = C^T,e^2{5/CoY'\ 

Thus, for i G [to — T, to), \T{t) — To] < 25. As is a 7r-invariant function which van- 
ishes at ^0, we see that F = (j).nO^ satisfies F{t) < and dF/dt\t± — ^'^F[g(j){(,±{t)) > 



Co > for t g [to — r, to) (c f. th e proof of Proposition 7.1 



account the form of </> and ( [A.3| ) , we deduce that for t G 
satisfies 



after ([7.35D). Taking into 

to-r,to),w(t) = w(7(t)) 



(A.7) 

Applying this with t 

(A.8) 



w(t) < Cie^S^\t 
T we see that 

w(io -T)< C2T6 



io|. 



Since w is independent of e and 5, we have deduced that there exists (5i > such 
that 

(A.9) to-^i < t < to ^w(t) < Cjt-tol^ 

In particular, under the same assumption, 

(A.IO) uj^{t)<C'\t~to\\ 

so 

(A.11) -'- '- ■ '^°-^ 



K^)-(^o + (T7^Mt))l<c"|t-to|^/^ 



Since Wot ^ and T{t) is C°°, 
to and 
(A.12) 
|77(t)-(Wo77)(t-to)|<qt- 



'Woii 
this shows that rj{t) is differentiable from the left at 



|3/2 



W077 = '^Hgj^ii), I G 7r-i(eo) arbitrary. 



Using this and the definition of wo we also conclude that 
(A.13) \z,[t) - [Woz,){t - to)\ < C\t - to|3/^ 

(A.14) \i^,{t) - {Woi^,){t - to)| < C\t - to|3/^ 

This proves the proposition for the tt- invariant functions t, zj, vj and "q, and indeed 
it provides a better error estimate. However, we still need to estimate yj. 



To do so, we consider the second term in w, see (10.4 
(A.15) \\y{t)\^-ii-^r,{tf\<C\t^to\^ Mo - (A - 



Thus, from (A.E) 



h{zo,yo)Y'^- 
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Taking into account (A. 12), we deduce that 

(A.16) r{t) = \y{t)\ 
satisfies 

(A.17) \r{t)' - A^lit - to)'\ < C\t - to\'^\ 
Thus, 

(A.18) |r(i) + 2^oit - to)\ < C\t - to\^/\ 

Hence, r{t) is also dilFerentiable from the left at to, and in particular 

(A.19) \y{t)\^r{t)<C\t~to\. 
Now, 



(A.20) 



77W, 



vit)' 






i4-\i^±itW 



A^o 



By ( |10J4| ), ( lAlOl ) and <^J§, 

(A.21) ||M±(t)P - Mol < C^dyWI + Mty^^) < C'\t - to|. 



Thus, by dA.lSD , 

(A.22) 

In particular, for each j we have 

(A.23) 

Let 
(A.24) 6, = y,/r, 

so 9j is a TT-invariant function away from Cq, and we have |0j| < 1. Also let 

(A.25) 



|yW-^M±WI'<C|i-tor 

Mo 



Mo 



9,(i)-^, fo-^i<t<to. 



r(i) 

By the inductive hypothesis, 9j{t) is differentiable for f e (^o ~ '5i,to) from both 
the left and the right and 



(A.26) 
with 
(A.27) 
and 

(A.28) 

Thus, 

(A.29) 



^i-'-w-'t-^w-w-'l 



dyj/dt\t± = 2^ij,±{t) 



dr/dtU^ = ^\y{t)\-\d\y\ydtU^)^2^. 



'\,^2ritri,,Mt)-'4^), 



dt 



ity 



so by ( |A.23|) and ( |A.12|) , 



,d^. 



(A.30) l^l*± - 2r(t)-iy,(i)(^2^(i)-i - r7(t)r(t)-^| < C|f - to^/^ 
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But, by (A.l?) and (A. 12), this gives 



(A.31) \^\^^\<c\t-tor/'. 

Integrating from to — 6i to to gives that 0j__(io) = hmt^t„_ ^j(^) exists and 

(A.32) \0j,-{to) - 9jit)\ < C'\t ~ to\'^''. 

Returning to the original notation, 9j = Vj/r, we see that 

(A.33) \y,{t) + 2tioe,,-ito)it -to)\< C'\t - to\^'\ 

so yj{t) is difFerentiable at to from the left. We then let 

(A.34) |_(io) = (0,z(io),r(<o),i^(io),-Aioej,-(to)). 

Then the compositions of the 7r-invariant coordinate functions yj, Zj, t and Vj 
with 7 are all differentiable from the left at io and the derivative is given by ^'^Hg 
applied to the appropriate coordinate function, evaluated at i^-(O). Note also that 
from (A. 22) and (A.33D we have 



(A.35) |M±W-M-(io)| <C|i-ior/', ie(to-<5i,to). 

Since a general smooth 7r-invariant function / has the form 

(A.36) /(y, z, r, ^, v) = /o(z, t, j/) + ^ yj7j (z, t, A*, i') + X! Vjykfjkiy, z, t, ^, i/), 

/oj /j, fjk all C°°, this finishes the proof of the proposition. 
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